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Abstract 



One of the biggest challenges of present day theoretical physics is to find a consis- 
tent theory of quantum gravity. There are however obvious problems in constructing 
a quantum theory of general relativity. It has already been shown in the seventies by 
't Hooft and Veltman that perturbative quantum gravity is non-renormalizable in four 
dimensions. However, this does not mean that it is impossible to formulate quantum 
gravity by use of field theoretic techniques, on the contrary, it suggests that special, 
so-called non-perturbative methods have to be applied. In this thesis we investigate 
the importance of causality in such non-perturbative approaches to quantum gravity. 

Firstly, causal sets are introduced as a simple kinematical model for causal ge- 
ometry. It is shown how causal sets could account for the microscopic origin of the 
Bel^enstein entropy bound. Holography and finite entropy emerge naturally from the 
interplay between causality and discreteness. 

Going beyond causal set kinematics is problematic however. It is a hard problem to 
find the right amplitude to attach to each causal set that one needs to define the non- 
perturbative quantum dynamics of gravity. One approach which is ideally suited to de- 
fine the non-perturbative gravitational path integral is dynamical triangulation. Without 
causality this method leads to unappealing features of the quantum geometry though. 
It is shown how causality is instrumental in regulating this pathological behavior. In 
two dimensions this approach of causal dynamical triangulations has been analytically 
solved by transfer matrix methods. 

In this thesis considerable progress has been made in the development of more 
powerful techniques for this approach. The formulation through matrix models and 
a string field theory allow us to study interesting generalizations. Particularly, it has 
become possible to define the topological expansion. A surprising twist of the new 
matrix model is that it partially disentangles the large-N and continuum limit. This 



i 



makes our causal model much closer in spirit to the original idea by 't Hooft than the 
conventional matrix models of non-critical string theory. 

The power of the here introduced techniques greatly extends the scope of analytical 
calculations for causal quantum gravity. This ranges from simple matter coupling to a 
new causal perspective on non-critical string theory. 
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CHAPTER 



Introduction 



Finding a consistent tlieory of quantum gravity is a notoriously liard problem. Such 
a theory should give a fundamental quantum description of space-time geometry with 
general relativity as a classical limit. However, more than ninety years have passed 
since Einstein's discovery of general relativity, and still very little is known about the 
ultimate structure of space-time at very small scales. 

Quantum field theory one the other hand has proven to be a marvelously successful 
way to describe three of the four fundamental forces of nature: the electro-magnetic, 
the weak and the strong force. For the gravitational force however one does not have 
a well-defined predictive quantum field theoretic description yet, but one does have a 
very successful classical field theoretic description in the form of Einstein's theory of 
general relativity. Since the other forces can be described by quantum field theories, it 
seems natural that there also exists a quantum theory for the gravitational field. 

Another reason which supports the existence of such a theory of quantum gravity 
is the fact that gravity couples to all forms of energy. Hence, one expects the energy 
fluctuations at small distances due to Heisenberg's uncertainty relations to induce also 
quantum fluctuations in the gravitational field. This leads to the prediction that space- 
time geometry has a highly non-trivial microstructure at extremely small scales close 
to the Planck length, Ip = y^KG^^^ 1.616 x IQ-^^m. 

There are however obvious problems in constructing a quantum theory of general 
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relativity (see [2J for an overview). It has already been shown in the seventies by 
't Hooft and Veltman that perturbative quantum gravity is non-renormalizable in four 
dimensions [3J. There are several ways to confront this problem: 

One of the most prominent approaches to quantum gravity is string theory (see 
standard text bool^s l4l El IH [7]). In this unifying theory the worldlines of propagat- 
ing particles are replaced by worksheets of propagating strings. The issue of non- 
renormalizability is thought to be resolved from the start in string theory, since the 
point-li[<e interactions which cause the divergencies are replaced by extended interac- 
tions which cannot be localized. Whereas string theory has recently had some suc- 
cesses such as the calculation of strongly coupled gauge theories using the AdS/CFT 
correspondence (see for instance [8J), as a theory of quantum gravity it still relies on 
supersymmetry and extra dimensions. 

Another approach is loop quantum gravity m HQ] in which one uses a Dirac (con- 
straint) quantization [11] of Ashtel^ar's connection variables [12]. Recent advances in 
this approach have been made by use of covariant spin-foam models [Tol |13|. These 
approaches suggest that the ultraviolet divergences can be resolved by the existence 
of a minimal length scale, commonly expressed in terms of the Plancl^ length /p. Finite- 
ness of gravitational entropy can be seen as evidence for such fundamental discrete- 
ness of space-time. 

A more radical approach is causal set theory [HI [IS which postulates fundamental 
discreteness from the outset. In this approach discrete sets of space-time events are 
viewed as fundamental objects underlying continuum space-times. Even though some 
interesting phenomenological models using causal sets have been proposed, there 
has been little progress in formulating a quantum path integral for causal sets. 

There are other approaches which define quantum gravity non-perturbatively. Renor- 
malizability of the non-perturbative theory could be achieved by a nontrivial fixed point 
scenario as described by Weinberg [Tel. One possible example of such a scenario 
is the exact renormalization group flow method for Euclidean quantum gravity in the 
continuum fVT] . 

Another famous attempt to non-perturbatively define quantum gravity is dynami- 
cal triangulations (DT), a covariant path integral formulation, in which quantum gravity 
is obtained as a continuum limit of a superposition of simplicial space-time geome- 
tries (see for instance [18]). Two-dimensional DT has also been used as a worldsheet 
regularization of the bosonic string. A particular strength of this approach is that it ex- 
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presses the theory as a statistical sum over entirely geometric objects. Since the path 
integral assigns a probability to each contributing geometry, we also call it a model of 
random geometry. To perform the statistical sum one usually starts off in the Euclidean 
sector including many "acausal" geometries. This leads to unappealing features of the 
quantum geometries in dimensions higher than two. It can be shown within the ap- 
proach of causal dynamical triangulations (CDT) that causality is essential to regulate 
this pathological behaviour. In CDT one starts off with a sum over Lorentzian geome- 
tries and then analytically continues to the Euclidean sector to perform the summation 
prg] . Numerical simulations indicate that CDT in contrast to DT leads to a well-behaved 
continuum limit in four dimensions (201 [2T] - 

In this thesis we focus on the implementation of causality in models of random 
geometry. 

In Part I of this thesis we concentrate on the approach of causal sets. In Chap. 
[2] causal sets are introduced as a simple and instructive model of causal random ge- 
ometry. After defining causal sets we show how they can correspond to continuum 
space-times in the continuum approximation. Although we are not able to define a 
path integral for causal sets, we describe a simple phenomenological model where 
causal sets are randomly embedded in a fixed space-time. From this model we show 
in Chap. |3] how causal sets could account for the microscopic origin of the Bel^enstein 
or Susskind entropy bound. 

It seems that to define random geometry in terms of a non-perturbative path inte- 
gral over geometries more structure is needed. Such structure is provided by causal 
dynamical triangulations (CDT) on which we focus in Part II, III and IV of this thesis. 

Part II gives an introduction to CDT as a model of causal quantum gravity with em- 
phasis on analytical results in two dimensions. In Chap. |4] we give a brief introduction 
into path integrals in general with emphasis on the gravitational path integral. In Chap. 
[5] we introduce dynamical triangulations (DT) as a non-perturbative definition of the 
path integral for two-dimensional Euclidean quantum gravity. At the end of this chapter 
we comment on the problems one encounters in higher-dimensional DT. In the next 
chapter, Chap.jej CDT are introduced as a model of two-dimensional causal quantum 
gravity. A comparison to the Euclidean model (DT) is given as well as a brief outlool^ 
to results in higher dimensions. In Chap. [7] several relations between DT and CDT in 
two dimensions are discussed. In particular, we describe how CDT and DT can be 
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related by respectively introducing or "integrating out" baby universes, i.e. spatial topol- 
ogy changes. In Chap. [8] we discuss CDT for the case of non-compact geometries. 
We see how a semi-classical background emerges from quantum fluctuations yielding 
a simple model of two-dimensional quantum cosmology. 

In Part III we introduce a generalized CDT model. This is formulated as a third 
quantization of two-dimensional CDT, i.e. a model in which spatial universes or strings 
can be annihilated and created from the vacuum. In Chap. [9] we show how by intro- 
ducing a coupling constant to the process of off-splitting baby-universes one can non- 
perturbatively solve for the disc function with an arbitrary number of baby-universes. In 
Chap.[To]this result is embedded in the frameworl< of a string field theory (SFT) model 
of CDT. As an application of the resulting Dyson-Schwinger equations (DSE) of the 
SFT we calculate amplitudes of higher genus. 

The last Part IV deals with the development of matrix models for CDT. After a brief 
introduction to matrix models in the context of DT in Chap. [lT[ we formulate a matrix 
model description for the contributing triangulations in CDT in Chap. [12} In the contin- 
uum limit one recovers the results of the generalized CDT model introduced in Chap. [9] 
giving a combinatorial interpretation of this model. In Chap. [13] we show that also the 
continuum CDT model is described by a matrix model. This is particularly interesting, 
since contrary to DT the size of the matrices N does not have to scale with the cut-off 
in the so-called double scaling limit. Further, the expansion in l/A/^ reorganizes the 
power expansions in terms of topology, where terms with powers A/^^''^^ can be iden- 
tified with the continuum surfaces of genus h. This mal<es the new matrix model much 
closer to the original proposal by 't Hooft for QCD than the conventional matrix models 
of non-critical string theory. 

In Chap. [14]we conclude this thesis by summarizing the results. Further, we com- 
ment on possible applications of the newly discovered matrix model to simple matter 
coupling in CDT. Interestingly, it might be possible to give a simple derivation of the 
Onsager exponents from CDT. 

Finally App. [A] and [B] provide supplementary information to the chapters on causal 
sets and causal dynamical triangulations respectively. 



5 



Most of the novel results of this thesis have been published in research articles. In 
particular, the following chapters are based on the following articles: 

• Chapter|3]on D. Rideout and S. Zohren, Class. Quantum Grav. 23 (2006) 6195= 

• Chapter|8]on J. Ambjorn, R. Janik, W. Westra, and S. Zohren, Phys. Lett. 6 641 
(2006) 94-98; 

• Chapter H] on J. Ambjorn, R. Loll, W. Westra, and S. Zohren, JHEP 071 2 (2007) 

• Chapter [To] on J. Ambjorn, R. Loll, W. Westra, Y. Watabil<i and S. Zohren, JHEP 
0805 (2008) 032^ 

• Chapter [12] on J. Ambjorn, R. Loll, W. Westra, Y. Watabiki and S. Zohren, Phys. 
Lett. B 670 (2008) 224-230| 

idem, to appear in Acta Phys. Polon. B 39 (2008) 3355-3364. 

• Chapter [13] on J. Ambjorn, R. Loll, W. Westra, Y. Watabiki and S. Zohren, Phys. 
Lett. B 665 (2008) .252-256, 

The remaining chapters are mainly introductory parts which review existing litera- 
ture. These are included to make the thesis self-contained, but do not represent new 
work done by the author of this thesis. 
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Part I 

Motivating causality: Causal sets 
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CHAPTER 



Causal sets: Discrete causal geometry 



In this chapter we introduce causal sets as a simple model of causal random geom- 
etry. After giving a mathematical definition of causal sets we show how those can 
correspond to continuum space-times in the continuum approximation. Even though 
we are not able to define a path integral for causal sets we describe a simple phe- 
nomenological model where causal sets are randomly embedded in a fixed space-time 
(T4|. In the next chapter we will use this model to give a microscopic understanding of 
the Susskind or Bekenstein entropy bound. 

2.1 I Introducing causal sets 

Causal order suggests itself as a fundamental principle for quantum gravity because 
of the enormous amount of topological and geometrical information which it contains 

m- 

There are several reasons to support the assumption of causal structure as pri- 
mary in quantum gravity. It has been shown that, given merely the causal relations of 
events in a space-time manifold, one needs only the volume measure to recover the 
full geometry in the continuum [221123]. 

Being faced with the non-renormalizability of perturbative quantum gravity in four 
dimensions 't Hooft already proposed in 1 978 several radical ideas to possibly define 
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a theory of quantum gravity [24], of which one was to view quantum gravity as a theory 
of causal sets of space-time events (see also [25]). This Idea was later picked up again 
by Bombelli et al. |p3I and further developed by many others to what is sometimes 
referred to as causal set theory (see for instance 1261 [271 for some reviews). 

IVIathematlcally speaking, a causal set is defined to be a locally finite partially or- 
dered set C = {C, -<), namely a set C together with a relation -<, called "precedes", 
which satisfy the following axioms: 



Transitivity: If x -< y and y -< z then x -< z, Vx, y,z e C; 



Irreflexivity: x 7^ x; 



Local Finiteness: For any pair of fixed elements x and z of C, the set of elements 
lying between x and z is finite, card{y|x ^ y ^ z} < 00, where cardX means the cardi- 
nality of the set X. 



Fig.[2J]shows a graphical representation of a causal set as an example. Of these 
axioms, the first two say that C = (C, -<) is a partially ordered set or poset. The last 
expression, local finiteness. Is related to the "radical idea" of this approach, namely 
that space-time is fundamentally discrete. The motivation for introducing fundamental 
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discreteness is that it provides us with a way to account for the missing volume informa- 
tion. However, from a phenomenological point of view there are also other reasons, for 
example the finiteness of gravitational entropy which suggests that space-time might 
be fundamentally discrete. 

In the following we summarize some basic definitions for intrinsic causal set quanti- 
ties of which most are "discrete" versions of analogous concepts used to describe the 
causal structure of continuum space-times (cf. App. |A.1| |. 

The past of an element x g C is the subset past(x) = {y e C \ y -< x}. This 
corresponds to the causal past J (x) in the continuum approximation. The past of a 
subset of C is the union of the pasts of its elements. The future of an element x e C 
is the subset future(x) = {y e C\x -< y} which respectively corresponds to the causal 
future -/+(x) in the continuum approximation. 

An important concept for the formulation of the entropy bound from causal set the- 
ory is that of a maximal element in a causal set C = {C, -<). This is an element which 
has no successors, i.e. an element x for which $y e C such that x -< y. The set 
of all maximal elements in C is denoted by max(C) = {x e C\$y e Cs.t.x ~< y}. 
In analogy, a minimal element is one which has no ancestors, i.e. an element x for 
which $y e C such that y ^ x and the set of all minimal elements in C is denoted by 
mm (C) = {xeC\$y e CsX.y -< x}. 

Towards a continuum approximation 

Having given the precise definition of a causal set one might asW the questions: How 
could one actually be able to formulate a theory of quantum gravity using causal sets, 
and how do causal sets relate to the known classical notion of smooth Lorentzian 
manifolds? In the following we give a definition of what it means for a causal set C to 
be approximated by a Lorentzian space-time {M,g). 

Consider a strongly causal space-time {M,g). The map (p : C ^ M from a causal 
set € = {€,-<) into a space-time {M,g) is called a conformal embedding if x ^ y 
0(x) G J (0(y)), Vx,y G C. Consider the Alexandrov neighborhood J^{p) n J (g), for 
every p,q e M, which forms a basis for the manifold topology of {M,g) if {M,g) is 
strongly causal, which we assume throughout, spn'nkling The map is called a faithful 
embedding or sprinkling if it has the following property: The number of elements n 
mapped into an Alexandrov neighborhood is equal to its space-time volume V times 



12 



Chapter 2. Causal sets: Discrete causal geometry 



the space-time density, up to Poisson fluctuations. Thus, the probability of finding n 
elements in this region is given by the Poisson distribution 

P(n)= ^^ \ , (2.1) 

where p = If'' is the density set by the fundamental length scale If in d dimensions. In 
other words, this means that in a space-time region of volume V one finds on average 
N=Vp elements of the causal set embedded into this region, and the fluctuations are 
typically of the order v^. Further, the properties of the Poisson distribution lead to 
local Lorentz invariance in the continuum space-time [2811291 . 

Using the above definitions, we say that a space-time {M,g) approximates a causal 
set C if there exists a faithful embedding of C into {M,g). This notion gives a corre- 
spondence between causal sets and continuum space-times. In terms of this corre- 
spondence continuum causal structure arises from the microscopic order relations of 
the causal set elements and the continuum volume measure of a region arises from 
counting the number of elements comprising this region. Nevertheless, one still has 
to prove the uniqueness of this correspondence (up to small fluctuations). In general 
the precise formulation of such a uniqueness proof is a difficult mathematical prob- 
lem. However, it has been proven to hold for the limiting case p i-^ oo [30] and certain 
progress has been made in the generalization to large but finite p plj. 

Discussion and outlook 

The ultimate aim of causal set theory is to formulate a theory of quantum gravity using a 
path integral approach, where the single histories are causal sets. Continuum physics 
would then be supposed to emerge in a suitable continuum approximation. However, 
since there is basically no information on spatial distances of space-time events, there 
is no obvious way to find a discretized version of the gravitational action in this context. 
And even if such an action would be found it is not clear how the sum over different 
causal sets in the path integral could be performed. 

A different path in the search for causal set dynamics has been taken by analyzing 
certain sequential growth models of the causal sets 1321 [331 |34l [35). These growth 
models aimed to determine an effective classical dynamics of causal sets. One inter- 
esting result is that the "early universe" phase of such growth models closely resem- 
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bles de Sitter space, at least in the way that volumes of causal intervals scale with their 
length [36J. 

From the above discussion we can say that it is still very ambitious to view causal 
sets as a model for quantum gravity, but nevertheless, it still gives a very simple 
and instructive model of causal random geometry. In particular, we can take a fixed 
Lorentzian space time and view it as a continuum approximation of a causal set. Using 
the Poisson distribution, as mentioned in the previous section, we are then able to ran- 
domly draw causal sets which can be faithfully embedded in this space-time. From a 
phenomenological point of view we can now study certain properties of such an ensem- 
ble, as will be done in the next chapter with view on entropy bounds. Even though some 
of these phenomenological models give interesting predictions ||37l [281 [381 [391 [40] . 
one should still treat them with a bit of care for the following reason: Already in the 
case of a quantum mechanical particle we know that even though the classical path is 
smooth and nice, single histories of the path integral look nothing like this and are not 
even described by differentiable functions. Hence, it is not a priori clear, that classical 
space-times should approximate a single causal sets history. Keeping this in mind we 
use causal sets as a phenomenological model of fundamentally discrete space-time. 
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Entropy bounds from causal sets 



In this chapter we propose a measure for the maximal entropy of spherically symmetric 
spacelike space-time regions in terms of the ensemble of causal sets approximating 
this continuum space-time [391141]. A bound for the entropy contained in this region is 
obtained from a counting of potential "degrees of freedom" associated to the Cauchy 
horizon of its future domain of dependence. For different spherically symmetric space- 
WWe regions in IVIinl^owsl^i space-time of arbitrary dimension, we show that this pro- 
posal leads, in the continuum approximation, to Sussl^ind's or Bel<enstein's well-l<nown 
spherical entropy bound up to a numerical factor. 



Entropy bounds in gravity 

In the history of general relativity there has been a long discussion regarding the 
thermodynamics of gravitational systems. One of the most famous examples is the 
Bekenstein-Hawking formula for blacl^ hole entropy 14211131 HUBS], stating that the en- 
tropy of a blac[< hole is given by a quarter of the area of the event horizon in Plancl^ian 
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unitsQ 

Sbh = ^. (3.1) 

This bound is truly universal, meaning that it is independent of the characteristics of 
the matter system and can be derived without any l<nowledge of the actual microstates 
of the quantum statistical system. 

Along this line there have been several generalizations of this entropy bound. One is 
Susskind's spherical entropy bound [46], stating that the upper bound for the entropy 
of the matter content of an arbitrary spherically symmetric spacelil^e region (of finite 
volume) is, 

■Smatter ^ ~^ , (3.2) 

where A is the area of the boundary of the region]^ Even though this spacelil^e entropy 
bound cannot be generalized to arbitrary non-spherically symmetric spacelil^e regions, 
there exists a generalization in terms of light-sheets, namely the Bousso or covariant 
entropy bound [491 ISQ]. IVIore precisely, let A{B) be the area of any (c/-2)-dimensional 
surface B, then the {d - 1) dimensional hypersurface L is called the light-sheet of S if /. 
is generated by light rays which begin at B, extend orthogonally away from B and have 
everywhere non-negative expansion. The entropy flux through the light-sheet is then 
bounded by 

S(/.) ^ (3.3) 



''As a convention we will throughout the discussion regarding entropy bounds work in Planck units 
with 

c = ks — Ci\i = h — 1, 

where c,kB,GN,h denote the speed of light, Boltzmann's constant, Newton's constant and Planck's 
constant respectively. In these units the Planck length is given by 

l^^^^^^^^l, for 03. 

^There was an earlier proposed bound by Bekenstein [47] stating that the entropy of any weakly 
gravitating matter system obeys Scatter 2nER, where E is the energy of the matter system and R the 
circumferential radius of the smallest sphere that contains it. If one further assumes that this bound is 
valid for strongly gravitating matter systems, then gravitational stability in four dimensions implies that 
2E ^ R and hence Smatter ^ 2ttER ^ A/A. One sees that in four dimensions the Bekenstein bound is 
stronger then the Susskind bound, however, \n d > A gravitational stability and the Bekenstein bound 
only imply that Smatter ^ {d - 2)A/8 gS]. Hence, the geometrical Susskind bound is arguably more 
fundamental. 



3.2 An entropy bound from causal set theory 



17 



This entropy bound is also widely regarded as evidence for the holographic principle 
(STl I46l [50], stating that the maximum number of degrees of freedom carried by L is 
given by A{B)/A. 

These entropy bounds suggest that an underlying theory of quantum gravity should 
predict the bounds from a counting of microstates (see for example ||52l). This verifi- 
cation of the thermodynamic laws is an important consistency check for any approach 
to quantum gravity. Further, the finiteness of the entropy might already give some 
indications about the actual microstructure of space-time. There is a semi-classical 
argument that the description of a quantum theory of gravity by a local quantum field 
theory in the continuum, in the absence of a high frequency cut-off, leads to infinitely 
many degrees of freedom in a finite region, and therefore to a divergence in the en- 
tropy of this region |38l |53]. The entropy bounds therefore suggest that space-time 
might posses a fundamental discreteness at scales of order of the Planck scale. Con- 
tinuum physics would then have to emerge from this fundamental theory when making 
a continuum approximation at large scales. This suggests that to obtain a theory of 
quantum gravity one does not have to quantize the metric fields of the continuum ge- 
ometries, but should rather find a quantum theory of the discrete structure underlying 
those continuum geometries [54]. 

In the following we show how using causal sets as a phenomenological model of 
causal random geometry one can derive a notion of maximum entropy from a count- 
ing of potential horizon "degrees of freedom" of the fundamental theory reminiscent of 
former ideas in the context of black hole entropy [55', '38]. Using this measure we for- 
mulate an entropy bound for spherically symmetric spacelike regions within the causal 
set approach. We then show that in the continuum approximation, for different spheri- 
cally symmetric spacelike regions in Minkowski space-time of arbitrary dimension, this 
leads to Susskind's spherical entropy bound up to a numerical factor. 

3.2 1 An entropy bound from causal set theory 

In the previous chapter we introduced causal sets as a phenomenological model of 
causal random geometry. As explained earlier space-time volume arises from a count- 
ing of fundamental space-time elements. In the following we show how entropy bounds 
could arise from a counting of potential horizon "degrees of freedom" at the fundamen- 
tal level, giving a microscopic origin for Susskind's spherical entropy bound. 
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Figure 3.1: Illustration of a spacelike hypersurface r, its boundary B{T.), its future 
domain of dependence and its future Cauchy horizon H+{T). 

As already mentioned in the previous section, the spherical entropy bound states 
that the entropy of the matter content of a spherically symmetric spacelike region T (of 
finite volume) is bounded by a quarter of the area of the boundary of Z in Planck units 



in full units S ^ AkBC^/{^GNh), where A = Vol(i3(Z)) is the area of the boundary of this 
region. 

Surprisingly, the maximum entropy in ( |3.4| | can be determined without any knowl- 
edge of the microscopic properties of the thermodynamic system. A theory of quantum 
gravity however should be able to deduce ( |3.4| | purely from a counting of the fundamen- 
tal degrees of freedom at the microscopic level. (One may wonder in what manner a 
counting of degrees of freedom measures the entropy of a system. In a discrete con- 
text "degrees of freedom" are generally finite, and a state counting can be expected to 
yield something proportional to the exponential of the number of degrees of freedom. 
Thus measuring the entropy as the logarithm of the number of states can be seen to be 
equivalent to counting the number of degrees of freedom of the system.) In the follow- 
ing we want to give a notion of those fundamental degrees of freedom in the context of 
causal set theory leading to the formulation of an entropy bound within this approach. 

Consider a spherically symmetric spacelike region of space-time, potentially con- 
taining some matter system. We assume that the "back reaction" of the matter content 




(3.4) 



3.3 Evidence for the claim 
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upon the space-time geometry can be neglected, so that different states of the matter 
system lead to the same spherically symmetric space-time geometry. The entropy of 
that system must eventually "flow out" of the region by passing over the boundary of 
its future domain of dependence. But because space-time Is fundamentally discrete, 
the amount of such entropy flux Is bounded above by the number of discrete elements 
comprising this boundary. This Is a fundamental limit imposed by discreteness on the 
amount of information flux which can emerge from the region, and therefore on the 
amount of entropy which It can contain. Thus we have argued that an entropy bound 
arises from causal sets. 

Proposal Consider a spherically symmetric spacelike hypersurface T of finite volume 
in a strongly causal space-time M of dimension d ^ 3. Denote the future domain of 



dependence of this hypersurface by D^{T) (c.f App. \A1\ and Fig. [STTp. LetC = (C,^) 
be a causal set which can be faithfully embedded into D^{T). Then the maximum 
entropy contained in Z is given by the number of maximal elements ofC, 

5max = card{max(C)}. (3.5) 

Claim This proposal leads to Susskind's entropy bound in the continuum approxima- 
tion, 

Smax = (3.6) 

where A = Vol(S(r)) is the area of the boundary of this region T, if the fundamental 
discreteness scale is fixed at a dimension dependent value to be calculated. 



Evidence for the claim 

In this section we provide analytical and numerical evidence for the claim that the en- 
tropy bound ( |3.5| ) leads to Susskind's bound in the continuum. In all discussed exam- 



ples we consider {d - l)-dlmensional spherically symmetric spacelike hypersurfaces T 
In c/-dlmensional Minkowski space M'', and calculate the number of maximal elements 
in its domain of dependence C'+(Z). Since we assume that D+(Z) arises as a con- 
tinuum approximation of a causal set C, we know from the phenomenological model 
described above that the elements of C are faithfully embedded into D+(Z) according 
to the Poisson distribution ( |2.1| |. Hence the expected number of maximal elements in 
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D+(Z) is given by 

{n)=p [ c/x^exp{-pVol(j+(x)f|D+(Z))|, (3.7) 

where again p is the fundamental density of space-time. On the right-hand-side of 
pTZ) ) one integrates over all points x g D+(5I), where every point is first weighted by 
the probability of finding an element of C embedded at this point and further weighted 
by the probability of not finding any other element of C embedded in J+(x) f| 0^(51). 
Note that because of the fundamentally random nature of the discrete-continuum cor- 
respondence in causal sets we calculate the expected number of maximal elements in 
a sprinl^ling, even though the proposal ( |3.5| ) is phrased in terms of a fixed causal set. 



3.3.1 I The ball in Minkowski space-time 

In this section we want to calculate the expected number of maximal elements (n) in 
D+(E), where E is chosen to be a (d - l)-dimensional ball Sd_i(/?) with radius R in 
IVIinl^owsl^i space-time M''. Due to the spherical symmetry of the problem it is useful 
to introduce spherical coordinates x = [t, r, 0i, dd-i), where we choose the origin to 
be the futuremost event of D^{T). The volume element Vol ( J+(x) f| D+(Z)) is equal 
to the volume of the Alexandrov neighborhood of x and 0, Vold(r) = Vol(J+(x) (^-/"(O)), 
where we denote proper time by r = Vt^ - r^. Using the result for Volrf(r) as calculated 
in App. |A.2| and integrating out the spherical symmetry in ( |3.7| one obtains a general 



expression for the expected number of maximal elements in D+(Sd_i(/?)), 



(n) = p ^\^2V [ dt f dr r^-e-P°^(^^-a^ , (3.8) 

I (~2~/ "'O 

where the dimension dependent constant is defined in App. |A.2[ 

In the following we evaluate for various dimensions by analytical and numerical 
methods. 



2+1 dimensions 

\n d = 3 dimensions one can explicitly evaluate ( |3.8| . It is useful to express the result 
in terms of the expected total number of space-time elements faithfully embedded into 
D+(S2(/?)), N = pV, where V = fR^ is the volume of the domain of dependence 
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D+(S2(/?)). The expected number of maximal elements is 

(n) = 8 (e-^ - l) - 2A/Ex + 4r (3.9) 

where E„(x) is the exponential Integral defined by 

POO 

Enix)= / t-"e-'''dt. (3.10) 



We also measured the expected number of maximal elements numerically by "sprin- 
kling" a causal set into a 2+1 -dimensional 2x2x1 square box, which contains D+(S2(1)). 
By sprinkling we mean simply selecting Nc elements at random with uniform distribu- 
tion within the box, and computing the causal relation between each pair from the 
Minkowski metric. For each of 100 trials / = 1 ... 100, we deduce the set of elements 
which fall within D+(S2(1)), compute its cardinality A/,, and count the number of such 
elements n, which are maximal within that region. From these we compute the sample 
mean and its error, and repeat this computation for a range of values of Nc. These 
computations were greatly facilitated by utilizing causal set and Monte-Carlo toolkits 
within the Cactus computational framework [56]. 

The plot of the expected number of maximal elements (n) as a function of N for the 



unit disk Is shown in Fig. [3^ on a logarithmic scale. The agreement of analytical and 
numerical results justifies the numerical methods. 

For a large number of elements A/^oo we can use the asymptotic expansion of the 
exponential integrals 

En(x) oc — (l - - + ..) for Ixl ^ oo (3.11) 
X V X / 

yielding 

(n) =4r ( - 1 ... for A/^oo, (3.12) 



^3, 

where + ... are lower order terms in A/. The asymptotlcs are also displayed in Fig. 
|3.2| together with the full expression for the expected number of maximal elements. In 
terms of the density p and the radius R of S2{R) this result reads 
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Ball in 2+1 dimensions 



Analytic result 
Numerical data 
Asymptotics 




total number of elements N 

Figure 3.2: Shown is the expected number of maximal elements (n) as a function of 
the total number of elements In the domain of dependence of the unit disk In 2+1 di- 
mensional Minkowski space-time, on a logarithmic scale. Besides the analytical result 
and Its asymptotics, data points with error bars from Monte-Carlo simulations are also 
shown. 



up to lower order corrections. It Is Important to see that (n) oc/4/4, where A = 2nR is the 
length of the boundary of S2(/?). This is highly non-trivial, as one can see by looking at 
a snapshot of a numerical simulation (Fig. [3^. There one observes that the maximal 
elements do not align along the one-dimensional boundary B{S2{R)). Instead they are 
distributed along a hyperbola close to the two-dimensional Cauchy horizon /-/+(S2(/?)), 
with a density of maximal elements which decreases with distance from the center. 
Hence the fact that the expected number of maximal elements is Indeed proportional 
to the length A of the one-dimensional boundary B{S2{R)) for large A already gives 
very non-trlvlal evidence for the proposed entropy bound. 

To have the precise confirmation of the constant in (n) oc A/A one has to 
choose the fundamental length scale If = p^/^ to be If = 16v^/(35/^r(l/3)) ^ 4.41. 
This gives support to the belief that the fundamental discreteness scale of space-time 
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Figure 3.3: A snapshot from a simulation showing N = 7806 space-time elements 
forming the domain of dependence of the two-dimensional ball in 2+1 dimensional 
Minkowski space-time (illustrated in red) and n = 126 maximal elements therein (illus- 
trated in green). 



is of the order of the Planck length. 

Hence, using this value for the fundamental discreteness scale we have confirmed 
our claim, namely in the continuum approximation we have 



(3.14) 



or Srriax = kB2'K R / {A ^/ HGn / c^) in full units, when we set If = 16^/(35/^r(l/3)). One 
now sees that the assumption of large but finite N used to obtain the asymptotic be- 
havior is equivalent to saying that A is much larger than the Planck length. However, 
one can see that even for relatively small values of the length A, such as 10^ in Planck 
units, the approximation of (n) by its asymptotic expansion is already very accurate. 
For even smaller values of the length scale the Planck corrections become significant. 
However, as one can see in Fig. [3^ the corrections always decrease the expected 
number of maximal elements, such that (n) never exceeds the bound ( |3.14| |. 
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3+1 dimensions 



Clearly from a physical point of view the evaluation of ( |3.8| ) in 3+1 dimensions is the 



most interesting case. For d = 4 one can write ( |3.8| ) as follows, 



[n) = 2TTp dt dz^/W^e~P^'^' (3.15) 
JO Jo 



One can perform the first integration in ( |3.15[ ) by using the following integral relation 
15ZJ 



Jo 



X 



X nFn[-, ; , ;pu ,3.16) 

\n n n n n n J 

for 3?(/i) > 0, 3?(z/) > and n = 2, 3, where u) denotes Euler's beta function and 
pFq{3i, ap\ bi, bq\ z) is the generalized hypergeometric function defined through 

and (a)n = r(a + n)/V{a) are the usual Pochhammer polynomials. The second integra- 
tion, namely the one of the generalized hypergeometric function, can be obtained by 
use of the following relation [57] 

^""^ pFq{3i, api 6i, bq;z)dz = — p+iF„+i(ai, Bp, a; bi, bq, a + 1; z). (3.18) 

a 

Expressed in terms of the number of causal set elements sprinkled into D+(S3(/?)), 
N = pV, where the volume of D+(S3(/?)) is given by V = ^R^, the final result reads 

(n) = A/3F3Q,l,l;^,^,2;-^A/). (3.19) 

For a large number of elements N we can use the asymptotic expansion of the 
generalized hypergeometric functions (cf. App. |A.3[ ), yielding the asymptotic expression 
for the expected number of maximal elements, 

(n) = 3V2^A/^ + ... for/V-^oo. (3.20) 
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Figure 3.4: Shown is the expected number of maximal elements (n) as a function of the 
total number of elements in the domain of dependence of the 3-dimensional unit ball 
in 3+1 dimensional Minkowski space-time on a logarithmic scale. The plot shows the 
analytical result, its asymptotics, and numerical results from IVIonte-Carlo simulation. 



In terms of the fundamental density of space-time p and the radius R of S3(/?) this result 
translates into 

(„) = V6,i^^. (3.21) 

As in the previous case the number of maximal elements follows the right scaling ac- 
cording to the entropy bound, namely (n) oc A/A, where the area of the boundary of 
S3(/?) is given by A = AnR^. The factor of proportionality is a C(l) constant as in 
the previous case, supporting the assumption that the fundamental length scale If is 
proportional to the Planck length Ip. More precisely, to have an exact agreement with 
the Susskind bound the fundamental length scale in four dimensions will be given by 
//r = 1^6 in Planck units. Using this value for the fundamental length scale the asymptotic 
expansion for the expected number of maximal elements reads 

Smax = (n) = ^ (3.22) 
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or Smax = 47r/?2/(ec^/(4/iG) in full units. At first sight one might feel uncomfortable in ab- 
sorbing the order one constant into the fundamental length scale to exactly confirm the 
Sussl^ind bound. However, as discussed in the previous case, the scaling (n) oc/4/4 is 
nontrivial and already serves as a confirmation of the bound. Taking the phenomeno- 
logical law S^ax = (spherical entropy bound) as "data" gives the fundamental length 
scale in four dimensions to be 

If = ^ ^ 1.57 (in four dimensions). (3.23) 

In comparison to the previous case of 2+1 dimensions one observes that the funda- 



mental discreteness scale If depends on the dimension. In Sec. |3.3.1] we will derive an 
expression for the factor for arbitrary (even) dimensions, showing that this factor tends 
exactly to one as c/^oo. It is important to check that If is universal for all cases with the 



same space-time dimension. In Sec. |3.3.2| we will provide evidence that this is indeed 
the case. 

At this point it is interesting to note that a similar method of fixing the fundamental 
discreteness scale to obtain the right factor of proportionality in the entropy bound is 
followed in loop quantum gravity in the context of black hole entropy (cf. [58]). There 
one fixes the Immirzi parameter which can be regarded as a measure for the discrete- 
ness scale (through its relation to the lowest eigenvalue of the area operator) to obtain 
the right factor of a quarter in the black hole entropy. However, since there are several 
ambiguities in the relation between the Immirzi parameter and the fundamental dis- 
creteness scale in the sense one uses it in causal set theory, it is hard to compare the 
numerical values in any sense. 

4-(-1 dimensions 



In d = 5 dimensions one can also evaluate the integral in in a similar way to 
the calculation in 2+1 dimensions. As in the previous cases the result is expressed in 
terms of the number of causal set elements sprinkled into D+(S4(/?)), N = pV, with the 
volume of D+(S4(/?)) given by V = ^R^. The final result reads 

in) = U 32 (l-e-r.)-,3NE J 'L]-NE ^ 



3 V V / ^ Vl6/ 5 V16 



5(2A/)ir ( ^ ) -30(2A/)^r ( ^ ) ) . (3.24) 



3.3 Evidence for the claim 



27 



Ball in 4+1 dimensions 



Analytic result 
Numerical data 
Asymptotics 




total number of elements N 

Figure 3.5: Shown is the plot of the number of maximal elements (n) as a function of 
the total number of elements in the domain of dependence of the 4-dimensional ball 
in 4+1 dimensional IVIinl^owsl^i space-time on a logarithmic scale, its asymptotics, and 
numerical results. 



The plot of (n) as a function of N is shown in Fig.|33]on a logarithmic scale as well 
as the numerical results obtained from IVIonte-Carlo simulation and the asymptotic be- 
havior. As in the previous cases there is agreement between analytical and numerical 
results. 

For large N one can expand | |3.24| using the asymptotic expansion for the expo- 
nential integral ( |3.11| , yielding 

(^n) = ^r (1]{2N)I + ... for N ^ oo. (3.25) 



3 V5. 

Re-expressed in terms of the density p and the radius R of Sa{R) this result reads 

327ri/5 I 2 3 27r2/?3 

As for the lower dimensional cases this result gives the correct behavior (n) oc A/4, 
where A = In^R^ is the volume of the boundary of S^{R), providing further evidence 
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for the claim. Further, we can use the coefficient in ( |3.26| ) to fix the fundamental length 
scale in 4+1 dimensions If = p-^/^ to be If = 21^6^1/15/(58/15(5^^5)1/6^(8/5)1/3) ~ 1.22 
in Plancl< units. Using this fundamental length scale the number of maximal elements 
( |3^ reads 

Smax = (n) = , (3.27) 

or in full units S^ax = 27r^/?^/(g/(4(/iGA//c^)^/^), which confirms our claim. 



Generalizations to higlier dimensions 

In the previous sections we have seen that causal set theory can provide a fundamental 
explanation for Sussl<ind's entropy bound for the spacelil^e hypersurface Sd-i{R) in 2+1 , 
3+1 and 4+1 -dimensional IVIinl^owsl^i space-time. In this section we want to generalize 
these results to arbitrary dimensions. For the case of odd space-time dimension, it 
turns out that one can easily calculate the expected number of maximal elements, but 
one cannot write the results in a closed form. However, for even dimensions one can 
find a closed expression. 

As in the previous cases we will express the result of the expected number of maxi- 
mal elements in terms of A/ = p V, where V is the volume of the domain of dependence 
of Sd i(/?), given by 

d-l 

V ^ Vol (D+(S._i(/?))) = ^f(iry/^'- (3.28) 

Using the integration relations | |3.16| and ( |3.18| ) one can integrate for even dimen- 
sions, yielding 



{"> - ^ ( i °. 1; ^ ^.2: -2>-/v] . ,3.29) 



2 4 d d + 1 d + 3 2d -1 
5+1- 5+1 ■■■■ ■ —J- 

Note that this result is also valid for the case oi d = 2 dimensions, however it is not 
related to any entropy of the system, and is thus excluded from the proposal. 

For large N one can use the asymptotic expansion of the generalized hypergeomet- 
ric function (App. |A.3| to derive the asymptotics of the number of maximal elements for 



arbitrary even dimensions, yielding 



, ^ 2 7r(d — l) .,d-2 , ,„ 



3.3 Evidence for the claim 
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Figure 3.6: Shown is the plot of the fundamental length scale If in Planck units (/p = 
{HGn/c^Y'''^^/^) as a function of the dimensions d for even and odd dimensions and 
the analytic continuation. 



From this result one can see that the number of maximal elements scales like (n) ~/4/4, 
where A is the volume of the boundary of i(/?), i.e. 



A = \/o\{B{S,^,{R))) 



r(^) 



(3.31) 



Further, we obtain {n)=A/A precisely for the following value of the fundamental length 
scale, 



If 



2\ d 



\ 2 J 



d^sin(^)r(2^) 



1 

d-2 



(3.32) 



The analytic continuation of ( |3.32| as a function of the dimension is shown in Fig. 



3^ together with the explicit values for 2+1 and 4+1 dimensions as determined in the 
previous sections. One observes that the expression | |3.32| agrees with these values. 
This suggests that ( |3.32| also holds for arbitrary odd dimensions. For c/ = 2 the value 
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Figure 3.7: Illustration of the different hyperbolic spherically symmetric spatial hyper- 
surfaces Z parameterized by a, and the domain of dependence 



of ( |3.32| | diverges, since the Planck length /p = (/iG/v/c^)('' ^'/^ is not well defined in two 
dimensions. This also reflects the fact that the entropy bound only holds for dimensions 
d^3. For all values the fundamental length scale is of order of the Plancl^ length. 



.2 I Generalizations to different spatial hypersurfaces 

In the previous section we have derived the Sussl^ind bound for the case where the 
spacelil^e hypersurface was chosen to be a {d - l)-dimensional ball in d dimensional 
l\/lin[<ows[<i space-time. Further, from this we determined the fundamental discreteness 
scale of space-time. However, it is important to prove that the fundamental discrete- 
ness scale so determined yields the same entropy bound for all spacelike hypersur- 
faces of a certain dimension. 

In this section we show that the claim also holds for spacelike hypersurfaces in 
IVIinkowski space-time different from those discussed in the previous section. We in- 
vestigate hyperbolic spherically symmetric spacelike hypersurfaces given by 

t = -^r2 + a2(/?2-r2), ^ a ^ 1, (3.33) 



as shown in Fig. 



3.7 



together with its domain of dependence. For the special case of 



3.3 Evidence for the claim 
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a= 1 the spacelike hypersurface given by ( |3.33[ ) is equivalent to the {d- l)-dlmensional 
ball Sd-i{R) for which we have determined an analytic expression for the expected 
number of maximal elements (n) in the previous section. For other values 0<a< 1 one 
cannot determine the number of maximal elements analytically. In the following we will 
investigate this problem numerically for the physically most important cases of 2+1 and 
3+1 dimensions. 



2+1 dimensions 

For d = 3 dimensions we use Monte-Carlo methods to numerically obtain the number 
of maximal elements in the domain of dependence of the spacelike hypersurfaces T 
defined by ( |3.33| ) for different values of a as a function total number N of elements in 
the domain of dependence. Since all these spacelike hypersurfaces T have the same 
boundary B{T), it is useful to express (n) as a function of the length of the boundary 
A = 2txR. One can do this by using that N = pV, where the volume of the domain of 
dependence of T is given by 

V=^ — R\ (3.34) 
3 1 + a ^ ' 

Further, we use the value for the fundamental density of space-time as obtained in 



Sec. 3.3.1 , i.e. p=(35/2r(l/3)3)/(81927r). 



The results of the simulations are summarized in Fig. |3.8[ Shown is the expected 



number of maximal elements (n) as a function of the boundary length A. For the special 
case of a = 1 in ( |3.33| ) we know that the result is given by ( |3.9| ), where N is replaced 
by A/ = /4^/(2^^r(2/3)^). The asymptotics of this analytic result is given by (n) = A/A as 
shown earlier. For different values of < a < 1 the simulations show that even though the 
number of maximal elements as a function of A differs for small values A < 10^ in Planck 
units, in the expansion for large A, corresponding to the continuum approximation, all 
functions (n) for different a enter the same asymptotic expansion (n) = A/A, yielding 
Susskind's entropy bound. In addition, it shows that the prediction of the fundamental 
discreteness scale is universal in 2+1 dimensions, at least for all investigated cases 
of spacelike hypersurfaces in Minkowski space-time. A generalization to examples in 
curved space-time has still to be shown and will be investigated in future work. 
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Figure 3.8: Shown is the plot of the expected number of maximal elements (n) In 
the domain of dependence as a function of the length of the boundary A for differ- 
ent hyperbolic spherical symmetric spacelike hypersurfaces parameterized by a in 2+1 
dimensional Minkowski space-time. One sees that all functions approach the same 
asymptotic (n) = A/ A. 



3+1 dimensions 

As in the case of 2+1 dimensions we use Monte-Carlo methods to numerically obtain 
the number of maximal elements In the domain of dependence of the spacelike hy- 
persurfaces Z defined by | |3.33| for c/ = 4 and different values of a. Again, the result 
is expressed as a function of the area of the boundary of T, I.e. A = AtcR^. This can 
be done by using the relation N = pV and noticing that the volume of the domain of 
dependence of T is given by 



|_3 6(l-a2)i 

where we use p = 1/6 for the value of the fundamental density of four dimensional 
space-time. 



Vl^{2 - 5a2) + 3a' 




(3.35) 



3.3 Evidence for the claim 
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Figure 3.9: Shown is the plot of the expected number of maximal elements (n) as a 
function of the length of the boundary A for different hyperbolic spherical symmetric 
spacelil^e hypersurfaces in 3+1 dimensional IVIinkowski space-time parameterized by 
a. One sees that all functions approach the same asymptotics (n) = A/ A. 



The numerical and analytical results are shown in Fig. [3^ Displayed is the ex- 
pected number of maximal elements (n) as a function of the area A of the boundary. 
For the special case of a = 1 the analytical result was given by ( |3.19| ), with N replaced 
by A/ = A^l [2tx12^). The simulations show that for different values of < a < 1 all func- 
tions (n) enter the same asymptotics, (n) = A/ A, in the continuum limit, giving further 
evidence for the claim. For very small values of a the expected number of maximal ele- 
ments enters the asymptotic regime only very slowly. However, the upper value for A of 
the simulations is still very small (~10^^m). As in the lower dimensional case the sim- 
ulations show that the prediction for the value of the fundamental discreteness scale of 
four-dimensional space-time, namely If = a/6, is a universal quantity for this dimension, 
at least for all investigated cases of spherically symmetric spacelike hypersurfaces in 
IVIinkowski space-time. 
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Discussion and outlook 

Using the phenomenological model of causal set theory as introduced in the previous 
chapter we argued for a bound on the entropy in a spherically symmetric spacelil<e 
region from a counting of potential horizon "degrees of freedom" of the fundamental 
theory, namely the maximal elements of the future domain of dependence of the re- 
gion. It was then shown that, for different spherically symmetric spacelil<e regions in 
IVIinl^owsl^i space-time of arbitrary dimension, this leads to Suss[<ind's spherical en- 
tropy bound. The evidence was given in terms of analytical results for spatial {d - 1)- 
dimensional balls in c/-dimensional IVIinl^owsl^i space-time, and in terms of numerical 
results obtained from IVIonte-Carlo simulations for the case of hyperbolic spherically 
symmetric hypersurfaces in IVIinl^owsl^i space-time. 

So far results were given only in terms of examples in flat l\/linl<owsl<i space-time. 
Clearly this is a very restricted class and further evidence should be provided in differ- 
ent space-times. Spherically symmetric spacelil^e regions in curved space-time such 
as in Friedmann-Robertson-Wall^er cosmology will be investigated in future wiorW. 

Another step would be to formulate the proposal intrinsically in terms of order invari- 
ants, without mailing explicit reference to the continuum. Such a formulation may be 
fruitful in providing a fundamental understanding of Bousso's covariant entropy bound. 

Further worl< in progress in the direction of entropy bounds from causal sets is the 
implementation of entropy evaluation through linl< counting as proposed in ||38l|B3l for 
blacl< holes, to other situations such as dS space-time. 



Part II 

2D causal dynamical triangulations 
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CHAPTER 



4 



Path integrals and quantum gravity 



In Part I of this thesis we introduced causal sets as a simple model of causal random 
geometry. The remaining three parts of this thesis will be centered around the ap- 
proach of causal dynamical triangulations. As already discussed in the introduction, 
one way to confront the non-renormalizability of perturbative quantum gravity in four di- 
mensions is to find a proper non-perturbative definition of the gravitational path integral. 
Before describing dynamical triangulations as such, we first give a brief introduction to 
path integrals in general. Specifically, we comment on problems one encounters when 
formulating the gravitational path integral. 



Random paths and one-dimensional gravity 

Path integrals were first introduced by Dirac and Feynman ||59l 160] as a quantization 
scheme to first quantize physical systems. One of the most instructive examples is the 
free relativistic particle in D-dimensional IVIinl^owsl^i space-time M°. The amplitude of 
a particle moving from x to y can be expressed by the so-called propagator. 




(4.1) 
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Figure 4.1: Illustration of the path integral for a one-dimensional relativistic quantum 
mechanical problem, e.g. a propagating particle. One possible path of the configuration 
space (path space) is drawn. The "virtual" particle is propagating from x to y in a 
piecewise linear path of n steps each of size a. 



where the action is simply given by the mass m times the length of the path 



S[P{x,y)] = m[ dl. (4.2) 

Jp{x,y) 

The classical equations of motion are derived by choosing a specific parametrization 

x(0 : [0, 1] ^ M^, x(0) = x, x(l)=y. (4.3) 



The action now reads 



S[P(x,y)] = m f di^/WW^^ 
Jo 



(4.4) 



where we defined := dx^'/d^. The classical equations of motion can be easily ob- 
tained by varying the action 
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d k^" 







(5x^^(0 di |x| 

which are simply straight lines from x to y and reparametrizations of it. 



(4.5) 



4.1 Random paths and one-dimensional gravity 
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To be able to define the path integral ( |4.1| | one usually performs a Wick rotation by 
taking time imaginary 

t ^ -/V = -ix^. (4.6) 

One now deals with the corresponding Euclidean model in M^. 

One way to define the path integral (|4J} is to introduce a "geometrical" discretiza- 



tion scheme where one directly discretizes the length of the worldline (see for instance 
|18]). Therefore it is useful to work with the action ( |4.2| |, since its definition does not 
depend on parametrized quantities. Specifically, in the employed discretization each 
worldline consists of a piecewise linear path, where each edge is of size a (see Fig. 
4~T| |. Here a serves as a reparametrization invariant cut-off which will be taken to zero 



in the continuum limit. For a piecewise linear path of length n the action simply be- 
comes S = moan, where mo is the bare mass. Further, the path integral can be written 
as a sum over all possible paths connecting x and y. 

«y OD p n n 

/ vp{x. y)^J2 U ^(^ E -(y- ^))- (4.7) 

"^^ n=l ;=o /=0 

Here the ae, denote the vectors of the linear parts of the discretized worldline, where 
the e;'s are unit vectors in ]R°. The propagator can now be defined as follows 

Grr^oix. y) = E / n E ^' - - ^))- (4.8) 

f)=l /=0 (=0 

To perform the integration it is helpful to use the Fourier transform 

/oo /" " 

c/xe-'>(^-^)G.„(x,y) = ^e-"- / llde,e-'^^\ (4.9) 
n=l •' /=0 

The integral in the above expression can be expanded as 

j dee-'"P' = f{a\p\) = f{<d){l- c^{apf + 0{a'')), (4.10) 

where f (0) = VoI(Sd-i(1)) is given in ( |A.5| and c is an irrelevant constant. We can now 
perform the summation leading to 

oo ^ 

G..M=i:(^-""f(^ipi))°= i_,-...,(,ipi) - <^-^^> 
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One can define the continuum limit by talking the cut-off a to zero and at the same time 
tailing the number of steps in the discretized path to infinity. The latter is fulfilled when 
fi = moa is fallen to its critical value Hc = log f{0) at which Gmo(p) approaches its radius 
of convergence. Hence, we renormalize the bare mass mo such that 

^i- fic = m^c^a^ (4.12) 

where m is the physical mass. Using this renormalization we obtain 

G,(P) — (4-13) 

After a wave function renormalization the continuum propagator G^{p) becomes the 
well-l<nown expression for the Feynman propagator 

In terms of its Fourier transform this is precisely the solution of the field equation for 
the Green function of a Klein-Gordon field, 

{-dj + m^)Gm{x.y) = 5(°)(y - x). (4.15) 

In the above description we used a discretization scheme where we directly dis- 
cretized the length of the worldline. A different prescription is to use the parametrized 
action ( |4.4| | and to introduce a cut-off on the worldline coordinates (see for instance 



161]). However, since the square root in this expression is very difficult to handle, it is 
useful to \ooW at an alternative action for the free relativistic particle first introduced in 
I62j, 

S[x.g]i..y) = C d^./W)\{g-\^)x^ + m^), (4.16) 
Jo 

where g{Q is the internal metric of the one-dimensional worldline. The action ( |4.16| 
can also be viewed as the action of one-dimensional gravity coupled to D scalar fields 
x^. Solving for the equations of motion for g and plugging the result back into the 
equations of motion for x one finds | |4.5| ). This shows that the systems desribed by 
both actions ( [4^41 ) and ( |4.16| ) are equivalent on the classical level. Whether this is also 



true on the quantum level remains to be checl^ed. To obtain the propagator we now 
also have to integrate over the worldline metric g, 

G„(x,y) = / V[g] f'v.x e''5[^'^l(-), (4.17) 



4.2 Random surfaces and strings 
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where V[g] refers to the measure of the space of worldline metrics modulo diffeomor- 
phisms, i.e. coordinate transformations on the worldline. Let us in the following briefly 
sketch how to calculate the propagator ( |4.17[ ) (for details see |6T)). The first step is to 
perform the W\cW rotation. One proceeds by defining the integral over VgX for fixed g. 
Therefore one discretizes the worldline coordinates by forming A/ + 1 slices of arbitrary 
small parameter differences 

en=^n-^n-l, ^0 = 0, ^/V = 1- (4.18) 

The integration over VgX and the continuum limit ^ can now be performed in a 
similar manner as one usually does for the non-relativistic particle. The remaining in- 
tegral over V[g] can then be done after a suitable gauge fixing. The result is precisely 
given by | |4.14| |. This indicates that both actions ( |4^ and ( |4.16| not only lead to the 



same classical equations of motion, but also to the same dynamics on the quantum 
level. A comparison of the two discretization schemes further shows that the "geomet- 
rical" discretization employed above is computationally rather simple in comparison 
to the continuum method. Dynamical triangulations, as will be introduced in the next 
chapter, is an analogous "geometrical" discretization scheme for surfaces and higher- 
dimensional manifolds. 



Random surfaces and strings 

The natural generalization to the random paths, as considered in the previous section. 



are random surfaces or worksheets of propagating (closed) strings (Fig.|4.2||. In anal- 
ogy to the relativistic particle we expect those to describe the relativistic string. The 
simplest action for a surface F is given by its area 

S[F] = dA{F), (4.19) 



where A is a (cosmological) constant conjugate to the area. The action ( |4.19| ) is a direct 



analog of ( |4T2| ). The parametrized form of ( |4.19| ) is called the Nambu-Goto action 
and reads 

S[F] = \j d^^y/h, (4.20) 

where 

/, = det/,,, = det— — (4.21) 
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(a) (6) 

Figure 4.2: (a) A worldline parametrized by ^ is embedded into by the mapping 
xf". (b) Analogous to the worldline a two-dimensional worldsheet of a closed string 
parametrized by and ^ is embedded into M° by the mapping x^. 



is the determinant of the metric h^b, a, 6= 1, 2, of the embedded worldsheet induced by 
the mapping x. 

As for the relativistic particle we can also write an action in terms of the internal 
metric g^b, 3,b = l, 2, of the worldsheet. This is the so-called Polyal<ov action |i62J, the 
analog of | |4.16| ), 

S[g.x]= [ d'^^/\^\{g^'h,b + X). (4.22) 

J M 



where one should integrate over both g and x in the path integral [64]. Both ( |4.19| ) and 



( |4.22| | are equivalent on the classical level, however, it is not so obvious that this also 



holds at the quantum level (see for instance [18]). 

We will see in the following two sections how ( |4.22| | can also be viewed as the 
action for two-dimensional quantum gravity coupled to D scalar fields. An explicit "ge- 
ometrical" discretization to define the corresponding path integral will be discussed in 
the next chapter. 



4.3 A path integral for quantum gravity 
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4.3 I A path integral for quantum gravity 

in the previous section we discussed the Polyal^ov and Nambu-Goto action as simple 
generalizations of random paths to random surfaces. The natural action for a theory of 
quantum gravity in c/-dimensions is the so-called Einstein-Hilbert action, 

SEH{[g.,]) = / d'^,/MOm - R). (4.23) 

IdttG/v Jm 

where Gn is the Newton's constant, A the cosmological constant, R the scalar curvature 
of the metric g^t, and g = det g^t, the determinant of the metric. For manifolds with 
boundaries one also has to include the Gibbons-Hawking-Yorl< boundary term |65il66l 

([^^'^]) = «;;7- / ^'^'^v^^. (4.24) 

OTTb/v Jqm 

where h and K are the induced metric and curvature on the boundary. All physical 
degrees of freedom are encoded in the equivalence class [gsb] which is g^t, modulo 
diffeomorphisms, i.e. coordinate transformations. 

The expression of the gravitational path integral for the partition function is then 
formally written as 

Z(G/v,A)= [ P[g,fa]e'5^"(fc^]), (4.25) 

where the integration should be tal^en over the space of all closed geometries, i.e. 
the coset space Geom(M) = Metrics(M)/Diff(/W). Generalizations to geometries with 



'GHY 



boundaries are also possible. When trying to define the formal expression ( |4.25[ ) sev- 
eral problems arise: 

First of all, Vlg^t,] has to be defined in a covariant way to preserve the diffeomor- 
phism invariance. Unfortunately, to perform further calculations one would then have to 
gauge fix the field tensors which would give rise to Faddeev-Popov determinants L6ZJ 
whose non-perturbative evaluation is exceedingly difficult]^ 

The second problem is due to the complex nature of the integrand. In quantum field 
theory this problem is solved by doing the \N\ck rotation t t-^ -/V, where t is the time 
coordinate in l\/lin[<ows[<i space (x, t). Clearly, a prescription like this does not work out 



^ See (Ml [Ml IZQl for an evaluation in tine setting of two-dimensional Euclidean quantum gravity in the 
light-cone gauge. In (TTl a calculation for three- and four-dimensional Lorentzian quantum gravity in the 
proper-time gauge is presented. It is anticipated there that the Faddeev-Popov determinants cancel the 
divergences coming from the conformal modes of the metric non-perturbatively. 
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in the gravitational setting, since all components of the metric field tensor depend on 
time. Further t t-^ -/V is certainly not diffeomorphism invariant (as a simple example 
consider the coordinate transformation t for t > 0). Hence, the question arises: 
What is the natural generalization of the Wicl< rotation in the gravitational setting? This 
problem is often circumvented in the approach of Euclidean quantum gravity, where 
one does an ad hoc substitution 

g'SEwdg-il"'']) — > e^SEHilgifh)^ (4-26) 

In this substitution the path integral over Lorentzian manifolds M('°'' is replaced by an 
integral over Euclidean manifolds M^^"). Clearly, the path integral now also includes 
many "acausal" manifolds and we will see in the next chapter when discussing two- 
dimensional Euclidean quantum gravity what consequences this brings with it. The 
main interest of this thesis is however to analyze ways to define the gravitational path 
integral in a Lorentzian setting including a proper definition of the W\cW rotation. 

Another problem in the definition of ( |4.25| is the following. Since we are wor[<ing 
in a field theoretical context, some k\r\6 of regularization and renormalization will be 
necessary, and again, has to be formulated in a covariant way. In quantum field theory 
the use of lattice methods provides a powerful tool to perform non-perturbative cal- 
culations, where the lattice spacing a serves as a cut-off of the theory. An important 
question to asW at this point is whether or not the theory becomes independent of the 
cut-off. Consider for example QCD and QED on the lattice. All evidence suggests that 
QCD in four space-time dimensions is a genuine continuum quantum field theory. By 
genuine continuum quantum field theory we mean a theory which needs a cut-off at an 
intermediate step, but whose continuum observables will be independent of the cut-off 
at arbitrarily small scales. In QED the situation seems to be different, one is not able to 
define a non-trivial theory with the cut-off removed, unless the renormalized coupling 
Bren = (trivial QED). Therefore, QED is considered as a low-energy effective theory of 
a fundamental theory at the Plancl< scale. 

Finally, one could argue to also include a sum over topologies in the path integral 
( |4.25| . However, in higher dimensions such a sum over topologies is quite intractable, 
since for d > 3 there does not even exist an obvious classification of topologies in 
terms of a finite set of parameters. 

These are some of the reasons to first lool< at the simpler case of two dimensions. 



4.4 Two-dimensional Euclidean quantum gravity and Liouville theory 
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Two-dimensional Euclidean quantum gravity and Liouville the- 
ory 

One of the big advantages of the two-dimensional Einstein-Hilbert action is that the 
curvature term is a topological invariant. This follows from the Gauss-Bonnet theorem 
which states that 

/ «''ev¥(OI^ = 4vrx(M), (4.27) 

J M 

where x{M) = 2 - 2h - b, \s the so-called Euler-characteristic, h is the genus of the 
manifold M, i.e. the number of holes in the surface, and b the number of boundary 
components of M. 

Using the Gauss-Bonnet theorem we can write the partition function for two-dimensional 
Euclidean quantum gravity including a sum over topologies (using a rescaling of the 
couplings) 

oo 

Z{Gn. A) = J] e"^('^'')/^"Zft(A), (4.28) 
where the partition function for manifolds with fixed genus h is given by 

Z,(A) = / P[^]e-^(^'^) (4.29) 

JGeom(/V?,,) 

and 

S(g,A) = A\4, V,= [ d'^^MOl (4.30) 

Here Vg is the volume of the manifold M^. Eq. ( |4.28| ) is called the topological expansion. 
Note that since the sum is tal^en over closed surfaces, we have x{M) = 2-2/?. Often 
one is only interested in the genus zero contribution, i.e. the sum over surfaces with 
topology of S^. In this case we denote the partition function by Z(A) = Zo{^). 

If the path integral is tal^en over surfaces with b boundaries the Euler characteristic 
reads x{M) = 2 - 2h - b. Further, it is natural to include boundary terms in the action 

b 

S{g, A, Zi Z,) = AVg + Yl ^'^'^S). (4.31 ) 

1=1 

where Li{g) is the length of the /th boundary loop with respect to the metric g. Here 
the Z; take the role of boundary cosmological constants. One should notice that in two 
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dimensions the intrinsic boundary geometry is completely determined by its lengths 
Li{g). The partition function for the case of boundaries is given thus by 

l/l/(A,Zi Zb)= [ V[g]e-^^^''^-^' (4.32) 

JGeom(M) 

and is called the multi-loop amplitude. Instead of fixing the boundary cosmological 
constants Z, one could have alternatively also fixed the boundary lengths, 

r ^ 

l/l/(A,/.i U) = \ V\g\e-'^^^''^\{b[U-Lig)). (4.33) 

JGeom(M) ■_-y 



The length space versions of the multi-loop amplitudes are also called the Hartle- 
Hawking wave functions HZSJ. For the case of one boundary component we also call 
it the disc function. One observes that ( |4.33| ) and ( |4.32| | are related by Laplace trans- 



forms, 

/■oo ^ 

l/l/(A,Zi Z,)= / WdUe-^^'\N{KU U)- (4.34) 

So far we discussed only pure two-dimensional Euclidean quantum gravity, i.e. with- 
out any coupling to matter fields. If one would like to couple D scalar fields one 
should add the following term in the action 

Sm(^,x) = £^ c/^evW^^'^l^l^. (4.35) 

It is interesting to notice that SM^g, x) + S(g, A) is precisely the Polyakov action | |4.22| , 
discussed in the previous section. Hence we see that pure two-dimensional Euclidean 
quantum gravity corresponds to string theory in zero-dimensional target space, i.e. 
non-critical string theory. 

In the following section we explain how to define the path integral for two-dimensional 
Euclidean quantum gravity using a discretization scheme called dynamical triangula- 
tions. Before doing so let us briefly mention a different way of computing the gravita- 
tional path integral using continuum methods. This is done by fixing the metric to the 
conformal gauge 

g = e^g, (4.36) 

where g is the reference metric. One sees that the theory is now expressed in terms 
of a single field, the so-called Liouville field 0. There are interesting analytical results 
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for Liouville field theory on the quantized level (see for example lZ3l [Z4j). Further, 
whenever the continuum theory and the discretized theory (as described in the next 
chapter) can be compared the results agree [75]. However, contrary to what one might 
have expected the discretized theory allows one to compute many quantities which are 
not accessible from the continuum theory. 



48 Chapter 4. Path integrals and quantum gravity 



CHAPTER 



Two-dimensional dynamical triangulations 



In this chapter we introduce dynamical triangulations as a discretization of the Eu- 
clidean gravitational path integral. In the special case of two dimensions analytical 
results can be obtained. In particular, we show how both the Hartle Hawking wave func- 
tion as well as the gravitational propagator can be calculated explicitly within this ap- 
proach. The results agree with those of the continuum formulation of two-dimensional 
Euclidean quantum gravity through Louiville theory as briefly mentioned in the previous 
chapter. 



Geometry from simplices 

The use of discrete approaches to quantum gravity has a long history UZH]. In classical 
general relativity, the idea of approximating a space-time manifold by a triangulation of 
space-time goes back to the early work of Regge [77J and was first used in ||78l [79| 
to give a path integral formulation of gravity. By a triangulation we mean a piecewise 
linear space-time obtained by a gluing of simplicial building blocks. This one might 
think of as the natural analogue of the piecewise linear path we used to describe the 
relativistic particle. In two dimensions these simplicial building blocks are flat Euclidean 
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Figure 5.1 : Illustration of a positive (a) and negative (b) deficit angle at a vertex v. 



trianglesQ where flat means isomorphic to a piece of Euclidean space respectively. 
One could in principle assign a coordinate system to each triangle to recover the metric 
space {M,g^y), but the strength of this ansatz lies just in the fact that even without 
the use of coordinates, each geometry is completely described by the set of edges 
length squared {/?} of the simplicial building blocks. This provides us with a regularized 
parametrization of the space of all geometries Geom{M) in a diffeomorphism invariant 
way, and hence, is the first step towards defining the gravitational path integral. 

For the further discussion it is essential to understand how a geometry is encoded 
in the set of edges length squared {If} of the simplicial building blocl^s of the corre- 
sponding triangulation. In the case of two dimensions this can be easily visualized, 
since the triangulation consists just of triangles. Further, the Riemann scalar curvature 
R{x) coincides with the Gaussian curvature K{x) up to a factor of 1/2. There are sev- 
eral ways to reveal curvature of a simplicial geometry. The most convenient method is 
by parallel transporting a vector around a closed loop. Since all building blocks are flat, 
a vector parallel transported around a vertex v always comes back to its original orien- 
tation unless the angles 9; of the surrounding triangles do not add up to 2ti, but differ 



by the so-called deficit angle = 1tx - (Figure 5.1 . The Gaussian curvature 

located at the vertex y is then given by 

/<. = ^, (5.1) 



the next chapter, when discussion two-dimensional causal dynamical triangulations, we will con- 
sider a similar construction with Minkowskian triangles. 
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Figure 5.2: Shown is a triangulation and its dual triangulation (dashed line). Highlighted 
is a vertex v of the triangulation and its dual cell of volume V^. 



where \/^ is the volume associated to the vertex v; more precisely, the dual volume of 



the vertex as shown in Fig.|5.2[ Note that the curvature at each vertex takes the form 
of a conical singularity. Then one can write the simplicial discretization of the usual 
curvature and volume terms appearing in the two-dimensional Einstein-Hilbert action, 



\ d'x^\^\R[x) ^ J]e„ (5.2) 
I d^x^\Aetg\ ^ (5.3) 



where TZ = {If} is a triangulation of the manifold M described by the set of edge 
lengths squared. From this one can write down the simplicial discretization of the two- 
dimensional Einstein-Hilbert action, the so-called Regge action. 



where k is the inverse Newton's constant and A the cosmological constant. It is then 
an easy exercise of trigonometry to evaluate the right hand side of ( |5.4| | in terms of the 
squared edge lengths. 

From this notion of triangulations one can now write the gravitational path inte- 
gral as the integral over all possible edge lengths / Vl, where each configuration TZ is 
weighted by the corresponding Regge action ( |5.4| |. In this approach the triangulation 
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is fixed and only the length of the edges is varied. A potential problem one encounters 
with this ansatz is a possible overcounting of triangulations, due to the fact that one 
can continuously vary each edge length, as is clear from flat space-time. Further, one 
still has to introduce a suitable cut-off for the length variable /. Among other things, 
this motivated the approach of "rigid" Regge calculus or (Euclidean) dynamical trian- 
gulations (DT) where one considers a certain class T of simplicial space-times as an 
explicit, regularized version of Euclidean geometries, where each triangulation T e T 
only consists of simplicial building blocl^s whose edges have all the same edge length 
squared Ll = a^. Here, the geodesic distance a serves as the short-distance cut-off, 
which will be sent to zero in the continuum limit. 

Following early worl^ by Tutte HSOl [Ml [821 183) and Brezin et al. [84], DT have in a 
quantum gravity context first been introduced as an regularization of the bosonic string 
|851[86l|871[88l[89l. Further developments have been made in many contributions, for 
an overview the reader is referred to [18]. The motivation for using a discretization of 
Euclidean geometries instead of Lorentzian geometries is mainly to avoid the need to 
define a gravitational W\cW rotation. One simply starts off with a sum over Euclidean 
triangulations weighted by the Euclidean (real) Regge action. In the continuum limit this 
theory of (Euclidean) DT corresponds to two-dimensional Euclidean quantum gravity. 
It is important to notice that this theory is distinct from two-dimensional Lorentzian 
quantum gravity which we will introduce in the next chapter. 

Putting the above discussion into formulas we can give a definite meaning to the 
formal continuum path integral of two-dimensional Euclidean quantum gravity as a dis- 
crete sum over inequivalent triangulations. 



where 1/C(T) is the measure on the space of discrete geometries, with C(T) = 
|Aut(7)| the dimension of the automorphism group of the triangulation T. 

Now the problem of calculating the gravitational path integral reduces to a statistical 
physics problem. Here the weight e^R'=sge(7") jg ^j^g so-called Bolzmann factor and the 
sum together with the measure factor determines the entropy of configurations with 
equal weight. In the following section we solve this statistical problem analytically for 
the two-dimensional model. 




(5.5) 
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5.2 I The disc function 

in this section we calculate the disc function or Hartle-Hawking wave function of two- 
dimensinal Euclidean quantum gravity defined through DT. The disc function is defined 
as a sum over all triangulations with topology of a disc with fixed boundary of length /. 
In other words the disc function describes the amplitude for the creation of a universe 
from nothing [72]. Although rather "global" it is a simple observable which can be used 
to calculate other interesting observables. 



5.2.1 I Discrete solution 

As we are in two space-time dimensions we have seen that the curvature term in the 
Einstein-Hilbert action is a purely topological quantity. Hence, if we are interested in 
computing sums over triangulations of fixed space-time topology, the contribution from 
the curvature term reduces to an irrelevant normalization constant. 
The Regge action ( |5.4| is thus simply given by 



SRegge(r) = XJ2V. = Xa^N{T), (5.6) 

veT 

where A is the bare (dimensionfull) cosmological constant and A/(T) is the number of 
triangles in the triangulation T. The factor of comes from the volume of a single 
triangle where the order one constant has been absorbed into A. 
The expression for the disc function wi{g) now becomes 

rer(/) Ter(/) 

where we defined g = e as the fugacity of a triangle. The sum is fallen over all 
triangulations with boundary length / and a mark on one boundary edge. The mark 
is introduced for convenience to remove a symmetry factor of 1/1 in the sum which 
corresponds to the factor 1/C{T) in the above expression. We can further write ( |5.7| | 
as 

wi{g) = ^ wi^Ng'^, (5.8) 

N 

where m//,a/ is the number of triangulations with / boundary edges and N triangles. We 
now see that the problem of determine the disc function reduces to a purely combi- 
natorial problem, namely finding the number m//,/v of all possible triangulations with / 
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Figure 5.3: Shown is a regular triangulation with topology of the disc (a) and a unre- 
stricted triangulation with topology of a disc (b). 



boundary edges and N triangles. Fig. [5j3] (a) shows an example of one triangulation 
with / = 9 boundary edges and N = 16 triangles. This combinatorial problem was first 
studied by Tutte in 1 962 IM ISll El |83l. 

In the following we will not follow the procedure by Tutte but rather count the number 
of all possible triangulations of a wider class, so-called non-restricted triangulations. 
Non-restricted triangulations are more general than the ones considered above in the 
sense that they can include double links (see Fig.|5j3](b)). The double links are mainly 
included to facilitate the calculation and to be able to map the results to those of matrix 
models. The use of matrix models will be an important part of this thesis and we will 
discuss them in detail at a later stage. The justification for being able to include double 
links is that when taking the continuum limit the contribution coming from the double 
links will be "washed away" as they do not carry a factor of g. Further, using the results 
of Tutte it can be checked explicitly that both regular and non-restricted triangulations 
lead to the same continuum expression for the disc function (see for example [18J). 

When trying to solve a combinatorial problem it is often useful to introduce gener- 
ating functionals, 

oo 

w{g,x) = J2Mgy- (5-9) 

/=0 

Here x can be understood as the fugacity of a boundary edge or in other terms as 
being related to a boundary cosmological constant A, via x = e~''^'. 

Generating functionals are very useful tools especially since there are helpful to 
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(a) 





Figure 5.4: Elementary decomposition moves: If the marked edge on the boundary 
belongs to a triangle this triangle is removed (a), if the mark corresponds to a double 
link the double link is removed and the triangulation splits into two (b). 



find recursion equations for w{g,x), the so-called loop equations. Let us first make the 
observation that whenever one removes triangles from a triangulation for each volume 
decrement of one triangle the generating function is decreased by a factor of g and for 
each decrement in boundary length it is decreased by a factor of x. 

Starting with a triangulation one can now uniquely decompose it step-by-step to a 
single dot by the following two elementary moves: Firstly, if the marked edge on the 



boundary belongs to a triangle this triangle is removed (see Fig. [5^ (a)). Thus one 
boundary edge is removed and two new edges are created of which by convenience 
the one further counterclockwise is marked again. In this move the power of g in 
w{g, x) is decreased by one and the power of x is increased by one. The second move 



corresponds to the case where the mark belongs to a double link (see Fig. 5.4 (b)). In 



this case the double link is removed and the triangulation splits into two, each with a 
mark next to where the double link was attached. 

The step-by-step procedure described above leads to the following recursion rela- 
tion (see Fig.|5.5| 



w{g, x) = - {w{g, x) - wi{g)x - 1) + x^w^{x, g), 



(5.10) 



where the first term corresponds to the first move and the last term to the second 
move. The term wi{g)x + 1 in the bracket is subtracted to be in accordance with the 
initial conditions of the the recursion relation. One can now solve for w{g, x), however, 
the solution is implicit since wi{g) still has to be determined. 

To be in accordance with the notation adopted in matrix models and also in later 
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+ 



Figure 5.5: Graphical representation of the loop equation (5.10) for the disc function 



parts of this thesis one often rewrites this equation using a slightly different definition 
of the generating functional 

i^U,l/z) 



w{g,z) = wi{g)z- 



1-1 



(5.11) 



/=0 



Here z"^ takes the roles of the fugacity of a boundary edge. The formal solution to the 
recursion relation | |5.10| ) reads in terms of the newly defined generating function 

w{g,z) = \ {V\z) - ^{V'{z)Y-AQ[z)) , (5.12) 

where anticipating the discussion of matrix models in later chapters we introduced the 
notation 

V'[z) = z- gz\ Q{z) = 1 - w^{g) - gz. (5.13) 

The sign in front of the square root is determined by the initial condition wq^ = 1 which 
leads to the requirement w{g,z) = l/z + 0{l/z^). In principle one would now have 
to find a combinatorial solution for wi{g), however, there exists a way of determining 



an explicit solution for w{g, z) by making use of the analytic structure of ( |5.12| . There- 
fore we first solve | |5.12| | for the limiting case of g = 0, where only double links are 
present. In this case w{z) = w{g = 0, z) corresponds to the generating function of 
rooted branched polymers (see Fig.|5.6|). The implicit solution ( |5.12[ ) becomes explicit 
in this special case 

(5.14) 



w{z) = - i^z-Vz^-A' 
This is precisely the generating function of the Catalan numbers 

1/1/2/ (2/)! 



w{z) 



oo 



^2/+l ' 

/=0 



1/1/2/ 



(/ + !)!/! 



(5.15) 
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Figure 5.6: A rooted branched polymer created by gluing double links with one marked 
link. 

Hence we see that the number of rooted branched polymers with / links is given by the 
Catalan numbers 1/1/2/. 

Anticipating the general solution ior g ^ it is instructive to analyze the analytic 
structure of ( |5.14| |. The generating function w{z) is analytic in C with a branch cut 
on the real axis along the interval [c_(0), c+(0)] = [-2,2]. The endpoints of the cut 
precisely determine the radius of convergence of w{z) which is needed to find the 
critical points when taking the continuum limit later on. 

The key to finding the solution for the general case of g 7^ is to notice that the 
analytic structure of w{g,z) as a function of z cannot change discontinuously at g = 0. 
For g = the branch cut is determined by the following expression under the square 
root 

{V'{z)f - 4Q{z) = (z - c+(0))(z - c_(0)) = (z - 2)(z + 2) (5.16) 
For g 7^0 this expression becomes fourth order in z, 

{V'{z)f - AQ{z) = [z - (2 + 2g) + 0{g')] X 
[z + {2 - 2g) + 0{g^)] [gz - {1 - 2g') + 0{g')]\ (5.17) 

Since the analytic structure of w{g, z) should not change in the neighborhood of g = 0, 
this implies that the term under the square root must be of the general form 



{V'{z)f - 4(?(z) = (c(g) - gzf{z - c+(g))(z - c4g)), (5.18) 
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where we labeled the roots such that c-{g) < c+{g) < c{g). The solution for the disc 
function then reads 

w{g, z) = \[z- gz^ - {c{g) - gz)^{z-c^{g)){z-c^{g))) . (5.19) 

Matching the coefficients of ( |5.18[ ) in the expansion in z yields four equations which 
completely determine c^{g), c+{g), c{g) and wi{g). This gives the full combinatorial 
solution for the disc function in DT 



.2 I Continuum limit 

Having solved the gravitational path integral as a statistical sum over triangulations we 
now have to perform the continuum limit to extract the continuum physics. This is done 
by talking the lattice spacing a to zero while at the same time talking the number of 
triangles N and the discrete boundary length / to infinity, such that physical area A and 
length L stay finite. 

In terms of critical phenomena the continuum limit relates to a fine-tuning of the cou- 
pling constants, generally denoted by g, to the critical point gc of the phase transition. 
This means a divergence of the correlation length / in the scaling limit 

l{g) ~ ^ oo, (5.20) 

la 5c I 



while the cut-off a goes to zero as 

a(^)~ \g-gcr^O, (5.21) 

where z/ is a critical exponent. In terms of the cut-off a and the correlation length /, the 
physical correlation length is given by L = l{g) ■ a{g) and the continuum limit is defined 
as the simultaneous limit 

a{g)^0, /(g) ^oo for fixed L = l{g)-a{g). (5.22) 

From standard techniques of renormalization theory, we expect the couplings with 
positive mass dimension, i.e. the cosmological constant A and the boundary cosmo- 
logical constant A, to undergo an additive renormalization, 

A = ^ + A, A/ = ^ + Z, (5.23) 
a^ a 
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where A, Z denote the corresponding renormalized values. Introducing the critical 
values 

^, = e-^\ z, = e^^, (5.24) 

it follows that 

g = g,e-'"\ z = z,e'^. (5.25) 

The critical point for the coupling g can be extracted from the radius of convergence 
of wx{g) at which the number of triangles N diverges. Using the result for wi{g) which 
can be derived easily as described below ( |5.2.2| | one finds gc = 1/^/432. To obtain the 
critical value for z at which also the boundary length diverges we look at the radius of 
convergence of w{gc,z) as a function of z. From the analysis of the analytic structure 
we know that the radius of convergence is precisely given by the endpoints of the 
branch cut and hence we conclude that z^ = c+{gc) (in the following we omit the precise 
numerical values of the critical couplings, since they are not of essential importance for 
the discussion). 



Inserting the scaling relations ( |5.25[ ) and the critical values into the discrete solution 
of the disc function yields 



w{g, z) = wns + a3/2(Z - iv^)i/zWA + 0{3^), (5.26) 

where w^s = V'{z) /2 is the so-called non-scaling part of the disc function which consists 
of terms of lower order in the cut-off, more precisely an order one constant and a 
term proportional to aZ. Further, for convenience we rescaled A and Z by positive 
multiplicative constants to simplify the result. 

Introducing a wave function renormalization of the continuum disc function 

Wj,{Z) = a-^/^w{g,z) (5.27) 

we obtain 

1/1/a(Z) = (Z - lv^)y'7WA, (5.28) 

where we dropped the non-scaling part. This is precisely the result obtained from 
Liouville theory. 

We can now also go back to the physical length L The continuum limit of the 
defining equation for the generating function ( |5.11| is precisely the Laplace transform, 

POO 

1/1/a(Z)= / dLW;^{L)e-'-^. (5.29) 
Jo 
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Hence performing the inverse Laplace transform 



W^{L) 



2ni 



l/l/A(Z)e 



LZ 



of (5.281 we obtain 



1/1/a(/.) = /.-^/'(l + v^/.)e-^^ 



(5.30) 



(5.31) 



Further we can also calculate the resulting disc function without a mark on the bound- 
ary by simply dividing ( |5.32| by L, yielding 



where the (u) stands for "unmarked". 



(5.32) 



5.3 1 Geodesic distance and the two-loop amplitude 

5.3.1 I Defining geodesic distance and the two-loop function 

In the previous section we calculated the disc function of two-dimensional Euclidean 
quantum gravity using the framework of dynamical triangulations. This is in a sense 
a rather "global" observable and in the following we are interested in calculating more 
"local" observables. The quantity we calculate is the fixed geodesic distance two-loop 
amplitude which we also sometimes refer to as the gravitational propagator or simply 
propagator. 

To be able to calculate this quantity we first need to make some definitions of basic 
concepts. Firstly, we have to give a proper definition of geodesic distance. If we would 
reintroduce coordinates on the simplicial building blocks one could simply employ the 
standard continuum definition as used in general relativity. However, we have seen that 
the strength of this approach lies in the fact that we have a coordinate independent 
representation of geometries. Hence, we should better look for a geodesic distance 
definition intrinsic to dynamical triangulations. There are several such definitions. The 
one we employ here defines the geodesic distance d{e, e') between two links (edges) e 
and e' in a triangulation as the number of links of the shortest path between the corre- 
sponding links of the dual triangulation. Now we want to extend the above definition for 
the case of geodesic distance of boundaries of the triangulation, where each bound- 
ary consists of a collection of links. In particular, we define the geodesic distance of a 
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Figure 5.7: Fixed geodesic distance two-loop amplitude (S(/i, h'.g; t). 



boundary loop £ to a link e as mine'e£ d{e, e'). Finally, we say that a boundary loop £i 
has geodesic distance t to another boundary loop £2 if all linl^s / e £1 have a geodesic 
distance t to £2- One should notice that this definition is not symmetric in £1 and £2. 



As an example, Fig. [5^7] shows an amplitude with a boundary of length k and an- 
other boundary of length I2 at geodesic distance t to the first boundary. One observes 
that the geodesic distance induces a natural foliation or time-slicing for the triangula- 
tions by following the sequence of loops of geodesic distance r = 0, 1, 2, t to £ij^ 
Whereas the space-time topology of the triangulation is fixed to Si x [0, 1], the topology 
of the spatial sections changes in general as a function of geodesic distance r. We call 
these spatial topology changes. One should notice that spatial topology changes are 
naturally present in Euclidean dynamical triangulations. However, as we will argue in 
the next chapter, in a proper theory of Lorentzian quantum gravity the branching points 
of changing spatial topology relate essentially to points in the continuum geometry of 



^This notion of time essentially suggests the existence of a Hamiltonian for two-dimensional Eu- 
clidean quantum gravity. The formulation of such Hamiltonian dynamics lead to interesting results in the 

we will 



field of string field theory (SUl IHI We will not review this work, however, in Chap. [To 
propose an analogous string field theory formulation for Lorentzian quantum gravity. 
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anomalous causal structure. This is one motivation for studying triangulations without 
spatial topology change as we will do in the next chapter. 

We now have all ingredients at hand to define the fixed geodesic two-loop amplitude 

G{li,l2;g;t)= s'^''^' (5.33) 

TeT{k,i2.t) 

where the sum is taken over all triangulations with an initial boundary of length /i and 
a final boundary of length I2 at geodesic distance t to the initial boundary. For t=0 the 
initial and final loop coincide which leads to the following initial condition 

G{k,l2;g;t = 0) = 5i,,i,. (5.34) 

As for the disc function it will be useful to the introduce generating function 

00 

G{z,w;g;t)= J2 G{li,l2:g:t)z-''-^w-''-\ (5.35) 

/l,/2=0 



with the initial condition equivalent to ( |5.34| |, 



G(z,w;g;t = 0) = ^ — . (5.36) 

zw zw — 1 

Solving \5.33) is a purely combinatorial problem. In the following we describe two 
different methods to solve this problem. The first method of solving the geodesic two- 
loop amplitude is the so-called peeling procedure due to Watabiki [94]. It is very simple 
and instructive, since it essentially uses the above loop equations and only views them 
as a "time-dependent" process. The disadvantage however is that not all steps of the 
calculation are rigorous. The second method due to Kawai et al. [95J uses transfer 
matrices (see also p6l 197]). Though more complicated it is on the other hand better 
defined. Further, it introduces important techniques which will be essential to study 
problems arising in the next chapter. 

5.3.2 I The peeling procedure 

In Fig.|5T]we illustrate the fixed geodesic distance two-loop amplitude as decomposed 
in slices of thickness Ar = 1. In the spirit of the loop equation for the disc function, 
as discussed in the previous section, we can also decompose the triangulation by a 
peeling procedure, a sequel of elementary decomposition moves (see Fig.|5.8|). 
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Figure 5.8: Decomposition of the triangulation for the fixed geodesic distance two-loop 
amplitude by a peeling procedure . 



Similar to the loop equation for the disc function we can write down a recursion 
relation for the geodesic two-loop amplitude (see Fig.|5.9| 



G'ik. k; g: t) = gG{k + 1, k; ^; + 2 -1-2, k; g: t). (5.37) 

/=o 

Here the two decomposition moves are essentially the same as for the disc function. 
The factor of two comes from the two possibilities of having the final boundary on either 
of the two off-splitting triangulations. The prime on the expression on the left-hand-side 
means that the surfaces contributing to this quantity do not exactly have C2 at geodesic 
distance t from Ci. 

Rather than giving a precise definition of G'{li, l2\g\ t) one makes a proposal which 
approximates G'{k, h-.g; t) for large t. When comparing to the exact transfer matrix 
calculation in the next section we shall see that this approximation indeed leads to 
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+ 





+ 




Figure 5.9: The peeling decomposition moves are the same as already used in the 
loop equation for the disc function (recall Fig. [53) ): If the marked edge on the initial 
loop belongs to a triangle the triangle is removed. If it belongs to a double WnW then 
the double linl< is removed causing the triangulation to split. The off-splitting part of the 
triangulation with the final loop has topology of a cylinder, whereas the other part has 
topology of a disc. 



the correct result. The basic idea in approximating G'{li, h'.g; t) is that when applying 
the elementary decomposition moves /i times on h.g; t) we will on average get 
G{li. k',g', t - 1), since roughly speaking we peeled off a whole slice. Ignoring the fact 
that t is an integer we can then view G'{k, l2\g\ t)=G{k, h.g; t+ l//i). For large /i one 
can Taylor expand 

h-.g; t + l//i) = G{k, h-.g; t) + l/kd,G{k, h; g;t) + .... (5.38) 
Inserting this into ( |5.37D gives 

^G(/i,/2;^;t) = -/iG(/i,/2;^;t)+^/iG(/i + l,/2;^;t) + 

oo 

+ 2/i J2 Mg)G{li - / - 2, h; g; t), (5.39) 



or in terms of the generating functions 



-G{z, w; g-^t) = — [{z - gz^ - 2w{g, z))G{z, w; g; t)] . (5.40) 

Plugging in the result for the disc function | |5.2.2| | we see that the first term in the bracket 
cancels the non-scaling part of the disc function and we get 

^G(z, w; g; t) = -2^ [%, z)G{z, w; g; t)] , (5.41 ) 
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where we defined 

f{g,z) = w{g,z) - \{z - gz') = \{gz - c(g))v/(z-c+(g))(z-c_(g)). (5.42) 

Instead of trying to solve the differential equation and then take the continuum limit 
of the result, we will in the following take the continuum limit directly of the differential 
equation ( |5.41| | and its initial condition ( |5.36| . Therefore we extend the scaling relations 



(|5j25| to 

g = g^e-''^, z = z,e'^, w = z,e'^, (5.43) 

using the same critical values as described above. Inserting these relations into the 
initial conditions ( |5.36| l gives 



G,{Z,W-T = Q) = ^^, (5.44) 

where we introduced the following wave function renormalization for the geodesic two- 
loop amplitude 

G^{Z, W; T) = aG{z, w; g; t). (5.45) 

Further, inserting the scaling relations ( |5.43| and ( |5.45| into | |5.41| we observe that the 
geodesic distance should scale as 



T ~ tVa, (5.46) 

yielding the following continuum differential equation for the geodesic distance two-loop 
amplitude 

— Ga(Z, 1/1/; T) = -— [1/1/a(Z)Ga(Z, W; T)] . (5.47) 

One observes that time does not scale canonically as maybe expected. We will see in 
the next section that this is due to the fractal structure of the quantum geometry. 

This partial differential equation can be solved under the above initial condition, 
yielding 

where 7(7, Z) is the solution to the characteristic equation 

^z( r , z) = - i/i/a(Z( r , z)), z( r = o, z) = z. (5.49) 



66 



Chapter 5. Two-dimensional dynamical triangulations 



Equivalently we can write 



T 



dX 



Z{T.Z) 



(5.50) 



One can now insert the solution for the disc function ( |5.28[ ) into this equation and 
perform the integral, 



arctanli ( X + ^/K] 



H x = z 



(5.51) 



X=Z(T,Z) 



where a = ^/3/2. Inverting this formula gives 



Z(r,Z) = a^^Atanh ( ^^T + arctanh ( -^\/X + V^ 



A, 



(5.52) 



which together with ( |5.48[ ) yields the explicit solution for G^{Z, 1/1/; T). 

The geodesic two-loop amplitude is a useful quantity from which many other quan- 
tities can be calculated. For the discussion of the quantum geometry in the next sec- 
tion we will be especially interested in the dependence on the geodesic distance. For 
simplicity, we shall analyze this dependence by shrinking both boundaries to zero. The 
resulting quantity is called the geodesic distance two-point function and was introduced 
in [98]. It is much simpler than the full two-loop amplitude, but still contains all essential 
dependence on the geodesic distance. 

Let us first shrinl< the outer boundary to zero. Performing the inverse Laplace trans- 
form we get 

W^{Z{T,Z)) 



and hence 



Ga(Z,/.2;7-) 



G^{Z, L2 = 0; 7) 



-LZ 



1/1/a(Z) 
1/1/a(Z(T,Z)) 



(5.53) 



(5.54) 



1/1/a(Z) 

One could have also obtained this result directly from the Laplace transformed quantity 
( |5.48| ) by taking iimw^oo 1/1/'^Ga(Z, W; T), where the coefficient \s [3 = 1 which lead to a 
finite expression. The quantity Ga(Z, L2 = 0; T) is in a sense similar to W^{Z) with the 
difference that it has a marked point in the bulk which has a geodesic distance T to 
the entrance loop. This marked point comes from the boundary which was shrunken 
to zero. If we would integrate over T e [0, 00] this expression should correspond to the 
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disc function with a marked point anywhere in the bulk, given by -^^W^{Z). A quick 
check reveals that this is indeed true, i.e. that 

-^1/1/a(Z)~/ dTG^{Z,L2 = Q,T). (5.55) 







From ( |5.54| | we then get for the two-point function G^{T) ■= Ga(/-i = 0, /.2 = 0, 7), 



Ga(7) = j^j' ^'^^llzP^ = 1^a(Z(7, Z = oo)), (5.56) 



where we chose (3 = 3/2 to get a finite expression]^ Using ( |5.52[ ) one gets 



G.(r)^(a^)3^°^ (6.67) 
Sinn |v A / 

5.3.3 I The transfer matrix approach 

In this subsection we want to rederive the continuum result for the fixed geodesic dis- 
tance two-loop function as obtained by ( |5.48| ) in the previous subsection using the 
so-called transfer matrix approach. This exact method is more rigorous than the peel- 
ing method described in the previous subsection. Further, it employs useful techniques 
which will be of importance in the forthcoming chapters. 

It was argued in [95] that any triangulation T(/i, I2, t) can in principle be uniquely 
decomposed into triangulations of geodesic distance one, i.e. T(/i, I2, t = l), by slicing 
it in the way described above. This situation is shown in Fig. |5.7[ 



A nice feature of this construction is that it implies the following semi-group property 
or composition law for the propagator 

G(/i, h-.g; fi + t2) = Yl ^(/i' S' fi) h; g; t2), (5.58) 

/ 

or in terms of the generating functions 

/dx 

— G{z,x-^;g; ti)G(x, w;g; ta). (5.59) 

Fig. |5.10| gives an graphical interpretation of the composition law. Setting ti = 1 and 
t2 = f - 1 one obtains 

G(/i, h-.g; t) = Yl ^(/i' t = l) G(/, h-.g; t - 1). (5.60) 



/ 



^The different values for (3 when shrinking the initial and final loop to a point (i.e. 3/2 and 1) come 
from the fact that the initial loop was marked and the final loop not. 
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Here k-.g; t = 1) is the so-called one-step propagator which can be seen as the 
matrix element of the so-called transfer matrix T{g), 



G{li,l2;g;t = l) = {l2\T{g)\k). 



(5.61) 



Knowing the transfer matrix one can in principle iterate ( |5.60[ ) t-times to obtain the fixed 
geodesic distance two-loop amplitude 



G{k,l2;g;t) = {l2\r'{g)\k). 



(5.62) 



If one is however only interested in calculating the continuum expressions there ex- 



ists a shortcut to obtain (Sa(/-i, L2, T) directly from the continuum limit of ( |5.60[ ) without 
having to perform the iteration as we will see in the following. 

Our first task is now to calculate the one-step propagator G{li, k; g\t=l) by combi- 
natorial methods. So far we have been working with amplitudes which have a marl< on 
the initial boundary. In some cases it is however useful to work with amplitudes without 
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a mark on the boundary. The marked and unmarked propagators are related as follows 

G{k,h;g;t) = liG^"\li,l2;g;t), (5.63) 

where the (u) stands for unmarked. Remember that the marked amplitude has only 
the initial loop marked and not the final loop. In terms of generating functionals, \5.63} 

d 



reads 



G{z,w;g; t) = -^zG^'\z,w;g;t). 



For the unmarked propagator the composition law ( |5.60[ ) becomes 

G^"\li, l2;g; ti + t2) = J2' ^^"H/i. '■•g< fi) G("'(/. h'.g; t2), 



(5.64) 



(5.65) 



where the measure factor / in the sum is the direct manifestation of the symmetry factor 
C{T) in ( [5:5) . 

We will now follow the results of [95] to combinatorially determine G('')(z, w;g; t). 
It was shown there that every one-step propagator is made out of the following four 
elementary building blocks: 



'A 



g 




z^w 



g{zw{g,z) - 1) 



o 



w{g,z) 



(5.66) 
(5.67) 

(5.68) 

(5.69) 



Here the dotted lines refer to the initial boundary with weight 1/z per edge and the 
dashed lines refer to final boundary with weight 1/w per edge. Again each triangles 
carries a weight g. Further, the blobs represent the DT disc function w{g, z) as derived 
in SecEiMl 

Now we can write the one-step propagator as the composition of all these terms 



G(")(z,w;g;l) = 



log 1 



g gw{g,z) g{zw{g,z)-l) w{g,zy 



ZW 



Z^\N 



1(5.70) 
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Here the factor of l/k reflects the circular symmetry, i.e. 1/C{T). 

In the following we perform the continuum limit using the same scaling relations as 
in the previous subsection 



(5.71) 



with the same critical values. Inserting the scaling relations and the solution for (S(")(z, w; g; 1), 
i.e. ( |5.70| ), into ( |5.64| | we get to lowest order in a 



G{z, w;g; t = l) 



1 1 



1 d fW/,{Z) 



0{1). 



(5.72) 



aZ+W ,/adZ\Z+W^ 
Hence we introduce the following wave function renormalization for the geodesic two- 



loop amplitude to match the desired initial condition ( |5.44D , 

Ga(Z, W; T) = aG{z, w; g; t). (5.73) 
From the second term in the expansion of ( |5.72| | we see that time has to scale as 

T ~ tv^. (5.74) 



Inserting the full scaling relations and the solution (5.72| for the one step propagator 



into the composition law (5.59) we obtain 



:Ga(Z, W; T) = -— [1/1/a(Z)Ga(Z, 1/1/; T)] 



(5.75) 



dT ■ ' dZ 

This is precisely the differential equation for the fixed geodesic distance two-loop func- 
tion as derived in the previous subsection, i.e. ( |5.47| . Thus, we rederived the previous 
result obtained from the peeling method using the transfer matrix formalism. The solu- 
tion to the differential equation ( |5.75| is again given by ( |5.48[ ). 



5.4 1 Physical observables 

One of the simplest observables in two-dimensional Euclidean quantum gravity is the 
disc function or Hartle-Hawking wave function which describes the amplitude of cre- 
ation of a universe from nothing. The result was obtained to be 




(5.76) 
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or 

1/1/a(/-) = + \n\L)e-^^. (5.77) 

Having determined the dependence on the geodesic distance through the two-point 
function in Sec. |5.3.2| we can further calculate the Hausdorff dimension dn of the quan- 
tum geometry. This dimension estimator is formally defined as 

{V{T)) ~ 7^" for T ^oo. (5.78) 

Using the expression for the two-point function, i.e. | |5.57| , 

G^{T) ~ e-^^/" for r ^ cx). (5.79) 

we obtain 
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This equation reflects the fact that at large T the quantum geometry looks effectively 
one-dimensional, i.e. the volume is proportional to the length. Further one can ex- 
tract from this expression that the average spatial length of the quantum geometry at 
intermediate T behaves as 

(0 = ^^-^- (5.81) 
Hence for typical time scales of 7 ~ I/a/A we have 

{V{T))r^T' for7~^. (5.82) 



Thus we conclude that the Hausdorff dimension of two-dimensional Euclidean quantum 
gravity is given by dn = 4. This result is somehow unexpected since we started off 
with a theory of two-dimensional triangulations. In fact, this is an indication that the 



quantum geometry is fractal [991 |95l HOOi I101| . As we will see explicitly in Sec. [Tj 
the quantum geometry is dominated by a proliferation of baby universes at the cut-off 
scale, meaning that at every point in the geometry there splits off a geometry with disc 
topology. The infinite number of baby universes causes the increment in the Hausdorff 



dimension. This fractal structure can be nicely observed in Fig. 5.1 1 which shows a 



typical configuration of the path integral obtained from a IVIonte Carlo simulation [IJ- 



Discussion and outlook to higher dimensions 

In this chapter we introduced dynamical triangulations as a regularization scheme to 
define the Euclidean gravitational path integral. In particular, we considered the case 
of two dimensions in which analytical results can be obtained. We calculated the disc 
or Hartle-Hawking wave function, the fixed geodesic distance two-loop amplitude as 
well as the two-point function. Using these quantities we were able to analyze several 
properties of the quantum geometry. One very interesting observation is the fractal 
structure of the quantum geometry which manifests itself in a proliferation of baby uni- 
verses at the cut-off scale. This is related to the unexpected value of the Hausdorff 



dimension of c/h = 4. Fig. |5.1 1 [ illustrates this situation by showing a typical configura- 
tion of the path integral obtained from a IVIonte Carlo simulation [1]. 

It is interesting to note how this situation translates to higher dimensions. In an 
analogous manner to two dimensions one can define the gravitational path integral as 
a sum over d-dimensional simplicial complexes (see for example [18] for a detailed 
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exposition). The corresponding Regge action takes the following simple form 

ST{k,-2, k,) = k,N,[T) - k,_2Nd-2{T), (5.83) 

where Nd{T) denotes the number of cZ-simplices and A/^ 2(7") the number of [d - 2)- 
simpices in the triangulation T. The analog expression to ( |531 ) for the partition function 
is then written as the sum over closed c/-dimensional triangulations T^, 

Z[k,_2, k,) = V J^e-^-t''^--''''). (5.84) 

Since analytical methods are absent in c/ > 2, one has to perform IVIonte Carlo 
simulations of the quantum geometry (see for instance [18]). DT in three dimensions 
were considered in [1 02111 031 11 041 11 05| . The four-dimensional case was first analyzed 
in (TME0ZJL108J, followed by further progress in [iMIlIOlEIIlEIBEISlIlB]. The 



results for c/ = 4 are summarized by the following schematic phase diagram (Fig. 5.12 



"^Although we only consider the case of c/ = 4 here, the picture turns out to be essentially the same in 
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One observes that the model possesses an infinite-volume limit everywhere along the 
critical line k^{k2), which fixes the bare cosmological constant as a function of the 
inverse Newton constant k2 ~ G^y^ Along this critical line, there is a critical point 
k^. Below this critical point the geometries generically have a very large Hausdorff 
dimension dn ^ oo. This can be understood in terms of a situation where every point in 
the quantum geometry is effectively at Plancl^ distance to any other point. We call this 
the crumpled phase of the quantum geometry. Above k^ we find the opposite situation 
where a configuration contributing to the state sum lool^s lilies a branched polymer with 
Hausdorff dimension dn = 2 (see Fig. [5^. This is called the branched polymer phase 
of the quantum geometry. 

One could hope that right at the transition point k^ genuine extended triangualtions 
with a finite Hausdorff dimension might be dominating. For this to happen the phase 
transition at this point should be of second or higher order. Unfortunately, it was found 
to be of first l^inc^ which leaves us with little hope that a sensible continuum theory of 
Euclidean quantum gravity can be obtained from DT for d>2. 

This unsatisfactory situation is one of the reasons for incorporating causality in DT. 
We have already seen above that many of the unappealing features of the quantum 
geometry were related to the proliferation of baby universes at the cut-off scale. In 
a genuine Lorentzian frameworl^ we expect such "acausal" configurations to be sup- 
pressed and the path integral should be much better behaved. The development of 
such a genuine Lorentzian frameworl< by so-called causal dynamical triangulations is 
the subject of the next chapter. 



the case of ty = 3. 

^Evidence that the phase transition is of first order was given both in three dimensions I105M1151 as 
well as in four dimensions I116II117II1181 . 



CHAPTER 



Two-dimensional causal dynamical triangulations 



In the previous chapter we introduced the reader to the concept of dynamical triangu- 
lations (DT). As we discussed at the end of the chapter, DT seems to be incompatible 
with a well-behaved continuum limit in three and higher dimensions. This problem is re- 
lated to the failure of incorporating the Lorentzian signature into the gravitational path 
integral. In this chapter we see how this issue can be resolved in the frameworl< of 
Lorentzian quantum gravity defined through causal dynamical triangulations [19]. The 
presentation of this chapter closely follows | [T9ll119[[T20j . 

Incorporating causality in dynamical triangulations 

Inspired by early ideas of Teitelboim | |121[ 1122] . we want to construct a gravitational 
path integral, where we insist in starting from space-times with Lorentzian signature. 
However, to address the Lorentzian nature of the path integral it is essential to look 
at geometries that have an intrinsically defined notion of time. In the approach of 
causal dynamical triangulations (CDT), introduced by Ambjorn and Loll in 1998 [19] for 
the case of two dimensions^, this problem is addressed by studying piecewise linear 

^Further developments for two-dimensional CDT have been made in fT23] [T24l [T25l [126] (1271 [128] 
11291 11301 11311 11321 11331 which will be partically discussed in more detail in the following chapters. 
References to results in higher dimensions will be given in the last section of this chapter. 
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t + 2 



t+1 




Figure 6.1 : Section of a two-dimensional Lorentzian triangulation consisting of space- 
time strips of height At = 1. Each spatial slice is periodically identified, such that the 
simplicial manifold has topology [0, 1] x One sees that the lower strip with initial 
boundary length l{t) and final boundary length l{t + 1) consists exactly of /(t) up- 
pointing triangles and l{t + 1) down-pointing triangles. 



geometries that have an entirely time-sliced structure. By virtue of the time-slicing one 
can globally distinguish timelil^e and spacelike edges. We will see that this distinction 
is essential to define a consistent W\cW rotation. 

As in the case of DT one uses triangles at cut-off length to discretize the geometries. 
However, in contrary to DT, CDT uses IVIinkowsl^ian instead of Euclidean triangles. 
These have one spacelike edge of length squared = and two timelike edges of 
length squared /.^ = -a^. By construction one chooses triangulations with topology of 

X [0, 1]. Each triangulation consists of t time slices of height At=l. The geometry 
of an individual time slice is now determined by the composition of l{t) up-pointing 



triangles and l{t + 1) down-pointing triangles as illustrated in Fig. 6.1 



In this formulation the time t by definition takes the role of geodesic distance. Fur- 
ther, there are no spatial topology changes with respect to the time-slicing. Hence, 
the layered structure does not only introduce an intrinsic time-direction that enables us 
to define a Wick rotation, but it also excludes spatial topology changes which in the 
Lorentzian picture lead to singular points with ill-defined signature]^ 

Since the triangles are simply patches of flat IVIinkowski space, they are naturally 



^In the next chapter we will relax this constraint and we will show how one can make sense of spatial 
topology changes in a Lorentzian setting. 
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equipped with a local light cone structure. Furthermore, due to the sliced structure of 
the triangulations which glues spacelike to spacelike and timelike to timelike edges, 
the triangulation is naturally equipped with a global causal structure. By virtue of this 
global causal structure it is possible to define a Wick rotation for curved manifolds on 
the discretized level. In the triangulation context this is done by performing a Wick 
rotation on the timelike edges, i.e. U iU. The Wick rotation acting on the whole 
manifold is then defined as 

W : r'°^ = { T, Ll = a\L] = -a"} ^ T^" = { T, = a" , L] = a"}. (6.1 ) 

As in the continuum case this rotation should be treated with some care and one must 
show that the Lorentzian action defined with = -a^ can be analytically continued 



to the Euclidean action defined with Ll = a^. This is done in appendix |B.1[ where we 
discuss in more detail the Regge action for two-dimensional Lorentzian triangulations. 
It is shown there that under the above defined Wick rotation the Boltzmann weight for 
each triangulation indeed transforms as 

yy : e'^'^'^^seC'"'") g-SReggeC'""')^ (6.2) 

We shall see in the forthcoming section that the above defined kinematical structure 
reduces the computation of the path integral for the fixed geodesic distance correlator 
or finite-time propagator to a simple statistical mechanics problem. 



The discrete solution 

As for DT we fix the space-time topology of the triangulation to be of the form x [0, 1]. 
The Regge action is then simply proportional to the volume of the triangulation, 

SRegge(r) = Aa2A/(7), (6.3) 

where A is the bare cosmological constant and N{T) the number of triangles in the 
triangulation T. Note that an order one factor coming from the volume term has been 
absorbed into A. Since all slices have the same geodesic distance, the propagator or 
fixed geodesic distance two-loop function can be defined in a simple manner 

Gi"\k,k;t)= y e'-^^^'^{^), (6.4) 

Tcik.h.t) 
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where Tc denotes the causal triangulations with initial boundary length /i and final 
boundary length h constructed out of t time-slices. Again, Cr denotes the volume 
of the automorphism group of a triangulation. As for DT we add a superscript (u) to de- 
note the propagator with no marl<ed points on the boundary. As discussed above, after 
the Wicl< rotation the discrete sum over complex amplitudes is converted to a genuine 
statistical model with a real Boltzmann weight, 

Gl"Vi./2;0= E ^ e-^^^'^(^), (6.5) 
rc(/i,/2,t) 

where A and A differ by an order one constant due to the different volume of IVIinl<owsl<ian 
and Euclidean triangles (see App.jBJ}. As for the fixed geodesic distance two-loop am- 
plitude in the case of DT the propagator satisfies the following semi-group property or 
composition law (e.g. ( |5.65| )), 

g1"'(/i. /2; ti + t2) = 5^ / /; fi) gI^V. h; t2). (6.6) 

where again the measure factor / in the composition law comes from the periodic 
boundary conditions. Writing the composition law for ti = 1 we find the following 
relation for the one-step propagator, 

/2; t + 1) = J] / Gx{li. /; 1) Ga(/, t). (6.7) 

As for DT we introduce a generating function for Gx{li, h; t) 

Gl")(x,y;t) = $^x''y'Gl")(/c,/;t), (6.8) 

k,l 

where x and y can again be related to the boundary cosmological constants of individ- 
ual triangles, 

X = e-^'", y = e-^"'. (6.9) 
The analogous fugacity of a triangle is again given by, 

g = e-^'\ (6.10) 

The total Boltzmann weight of one strip can be determined by noting that for the 
one-step propagator in CDT a factor of gx is assigned for triangles that have the space- 
lil^e edge on the initial loop and a factor gy for triangles where the spacelil^e edge is 
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on the final loop. This relation is much simpler than in the case of DT. The one-step 
propagator is now easily computed by, 



oo ^ / oo oo \ 



(6.11) 



where the factor 1/k comes again from dividing by the volume C(T) of the automor- 
phism group for circular triangulations. Performing the summations in ( |6.1 1| ) we readily 
obtain 



f^^k\l-gxl-gy) 



- log 1 



gxgy 



{l-gx){l-gy) 



(6.12) 



From the inverse discrete Laplace transform of ( |6.12D it can be seen that the one-step 
propagator with fixed boundary lengths is given by 



(6.13) 



where the division by the volume of the automorphism group now maizes its appear- 
ance in the factor l/(/i + I2). Similar expressions for different boundary conditions can 
for instance be found in |1 1 9j . 

Instead of obtaining the continuum limit directly of the one-step propagator as 
done in the previous chapter we want to iteratively compute the finite-time propagator 
G{x, y; g; t) and then perform the continuum limit. To do so we rewrite the composition 
law \6.6) in terms of generating functions and obtain, 

f^-^- T. ^ G(x,z;g;ti)G(z',y;g;t2). 

j Ztti Ztti (1 - zz') 

Setting t2 = 1 and performing the contour integration over z one gets 



(6.14) 



G^''\x,y-g-t) = j 



dz' 



2tx iz' 



dz 



G^"\x,z-g-l) 



G^''\z',y-g-t-l). 



(6.15) 



z=l/z' 



Inserting the expression for the one-step propagator, i.e. ( |6.12[ ), yields 

G(")(^. y; t) = G^"' (r^' ^; t - 1) - G^"^ {g< y<g<t- 1) 



(6.16) 
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This equation can be iterated and the implicit solution written as 

G'"'(.,y;*; ^ log (^-^) - log {^^^) . (6.17, 

where Ft is defined iteratively by 

Ft{x)= f Fo(x)=x. (6.18) 

1 - gFt-i{x) 

The fixed point F of this equation as defined by Ft{x) = Ft_i(x) is given by 

By standard techniques one can use the fixed point to give the explicit solution to the 



iterative equation (6.181 



where we have defined 

At = l-F^'^\ Bt = F-F^'+\ Q = F^-F^\ (6.21) 

The explicit solution for the propagator is now obtained by substituting ( |6.20| in | |6.17| , 
yielding 

G^'\x, y; g-t) = - log (1 - Z(x, y; g; t)) , (6.22) 

where we have defined 

Z{x.y:g:t) = (l-^)- „ X - (6.23) 



(l - ftx) (l - |y 



The combined region of convergence of this result as an expansion in powers g'^x'^y'^ 
is 

\g\ < \ \A < 1. |y| < 1- (6-24) 
As in the case of DT we will use the radius of convergence to obtain the critical points. 
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The continuum limit 

To define the continuum limit of the theory we assume canonical scaling dimensions 
for the bulk and boundary cosmological constants as in the case of DT, 

g = g,e-''^'\ x = x,e-^^, y = yce-^^. (6.25) 

The critical values for these couplings are again determined by the region of conver- 
gence of the propagator, i.e. | |6.24| , yielding 

gc = 1/2, Xc = l, y, = l, (6.26) 

One can quickly convince oneself that these are the correct critical vales by looking at 
the average number of triangles in one strip 

(A/) = G(")(x,y; g; 1) = 1 - - + ^— — . (6.27) 

dg I- gx I- gy I- g[x + y) 

We see that the average number of triangles diverges at the critical point {gcXcYc) = 
(1/2, 1, 1). We hence arrive at the following scaling relations 

g = le-''''/\ x = e-^^, y = e-^^. (6.28) 
The continuum limit of the propagator can now be determined by inserting these scaling 



relations into ( |6.22[ ). To get a sensible continuum limit we require time to scale as 

T ~ 3t. (6.29) 

which yields the following result 

Gi:\x. Y.T) = - log (1 - Za(X, Y. T)) . (6.30) 

where 

Za(X, y, r) = (6.31) 

sinh^ VAr(X + VAcoth VA7)( V + VAcoth VA7) 

As in the case of DT the propagator undergoes a wave function renormalization 

G["\X,Y,T) = aG^"\x,y;g;t). (6.32) 
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Again the power of a is uniquely determined by the initial condition coming from the 
composition law, i.e. ( |6.14| |. 



The continuum expression for the propagator, where the boundary lengths are fixed 
instead of the boundary cosmological constants can now be obtained by an inverse 
Laplace transformation of ( |6.30| with respect to both X and Y, 



where /i(x) is a modified Bessel function of the first kind. The continuum version of the 
composition law ( |6.14| | for fixed lengths then reads 



POO 

G["\U,L2-T^+T2)= / dLLG['\u,L-T,)G["\L,L2,T2). 

Jo 



(6.34) 



An interesting observation is that in contrast to the Euclidean model the time or geodesic 
distance scales canonically. We will further comment on this when discussing the phys- 
ical observables of the model. 

6.4 1 Marked amplitudes 

To better compare to the results of the previous chapter for Euclidean DT it is important 
to introduce the same convention of a marl^ed edge on the initial boundary in CDT In 
terms of the propagator this convention reads, 

Ga(/i,/2;1) = /iG1"Vi./2;1)- (6.35) 

As before this marl^ing removes the measure factor in the composition law for the prop- 
agator (|6j6l 

Ga(/i, I2: h + t2) = Y, ^a(/i, /; ti) Ga(/, k; t2), (6.36) 

/ 

which corresponds to the following composition law for the generating function, 

Gx{x,z-'; h)Gx{z,y; ts). (6.37) 

Ztti z 

The one-step propagator with a marl< on the boundary can easily be computed from 



(6.12|| by taking a derivative with respect to the boundary cosmological constant. 



G(x, y ; g; 1) = x-f G^^\x, y ; g; 1) = (""^ (6.38) 

dx {1 - gx){l - g{x + y)) 



6.4 Marked amplitudes 
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Inserting this expression Into the composition law ( |6.36[ ) gives an iterative equation 
analogous to (|6.16[) , 



g 



l~ gx 



y;g;t-l 



(6.39) 



1-gx 

The explicit solution can be obtained from | |6.22| | by taking a derivative with respect to 
the boundary cosmological constant, 

G(x,y; g; t) = x^G^^\x,y; g; t) = „ T}'^ 'SP'''^ ^ (6-40) 

^ ^ dx ^ ' {At - Btx){At - Bt{x + y) + Qxy)' ^ ' 

where F, At, Bt and Q were defined In ( |6.19[ ) and ( |6.21[ ). 

To obtain the continuum marked propagator we could now apply the scaling rela- 
tions ( |6.28[ ) to this expression. Instead, one can also apply the scaling relations directly 
to ( |6.39| In the manner we proceeded for DT. Doing so one finds the following differen- 
tial equation for the continuum marked propagator. 



_^G,(X, Y. T) ^ - A 



W{X)Gf,{X, Y, T) 



(6.41) 



where 

W{X) = X^-A. (6.42) 

The notation W{X) is introduced In analogy to the Euclidean result, however. It should 
be noted that in the case of CDT W{X) is not the disc function. 
Solving the differential equation ( |6.41| | with Initial condition 

1 



Ga(X, Y; 7 = 0) 



gives 



G^{X, Y; T) 



X+Y' 
W{X{T,X)) 1 



W{X) X{T,X)+Y 
where X{T, X) is the solution to the characteristic equation 

dX{T,X) 



(6.43) 



(6.44) 



dT 



-1/1/(X(T,X)), X{T = 0,X) = X. 



Proceeding as In Sec. 5.3.2 we obtain. 



(VA + X) + e-2^^(VA-X) 



(6.45) 



(6.46) 
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Inserting this expression in | |6.44| gives the explicit solution for the continuum marked 
propagator. Further, one can check that if one marks the unmarked propagator ( |6.30| 
by taking a derivative with respect to the initial boundary cosmology constant X, the 
resulting expression precisely coincides with this result. Inverse Laplace transforming 
( |6.44| | with respect to X and Y yields the continuum marked propagator in the length 
representation, 

= sinhVAT Vm [^iiWKf) • ^'-"'^ 
which is simply ( |6.33| multiplied by U as expected. 

As in the case of DT we can also calculate the disc function from the propagator. 
In CDT it can be defined as follows from the finite-time propagator with one boundary 
shrunken to zero, 

1/1/a(/.) = / dTGA{L, L2 = 0, T). (6.48) 

Jo 

One notices that this definition differs slightly from the corresponding relation in the 
Euclidean setting, i.e. ( |5.55[ ). This is due to the fact that the mark coming from the 
zero-length boundary is always at latest T and cannot be anywhere in the bulk as in 
the case of DT. Therefore we do not have to differentiate with respect to -d/dA. 
Using ([633]) we get 

Ga(/-, L2 =0,T) = A^e-^Lcoth^T (g_4g) 

sinh vA7 



Inserting this into ( |6.48[ ) one readily arrives at the simple expression 

W^{L) = e"^^ (6.50) 

or in Laplace transform language 

^.(X) ^ (6.51) 
The corresponding result for the disc function without mark reads 

~L 

From | |6.49| ) we can also easily compute the two-point function for CDT, i.e. 

Ga(7-) = G^{U = 0,L2 = 0, T) = . (6.53) 

Sinn vA/ 



w!C\L) = (6.52) 



6.5 Hamiltonians in causal quantum gravity 
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Hamiltonians in causal quantum gravity 



As mentioned earlier the concept of time or geodesic distance directly relates to a no- 
tion of Hamiltonian dynamics. This can be seen explicitly when viewing the differential 
equation for the continuum marked propagator ( |6.41| as a W\cW rotated Schrodinger 
equation, 



_d_ 
dT 



G^{X,Y,T) = HxG^{X,Y,T). 



(6.54) 



Here Hx can be interpreted as the effective quantum Hamiltonian of the system and is 
given by 

Hx = (X2-A)J^ + 2X. (6.55) 



The effective quantum Hamiltonian in length space can be obtained from ( |6.54[ ) by 
inverse Laplace transformation, yielding 



The corresponding Hamiltonian for the case of unmarl^ed amplitudes reads 



(6.56) 



(6.57) 



These kind of Hamiltonians are analyzed in detail in App. |B.2[ where several properties 
such as the spectrum and eigenfunctions are computed. In particular, it is shown there 
that the Hamiltonian \<o.57) is self-adjoint with respect to the measure LdL This is also 
reflected in the corresponding semi-group property, where the integration is taken over 
the same measure. As we will see later for some comparisons it is useful to transform 
this Hamiltonian to one which is self-adjoint on a flat measure. This can be done by 
absorbing the measure in the wave functions corresponding to the initial and final loop, 
i.e. ipniL) = -^ipnlL). Commuting 



6.57 



with 1/VL gives 



(6.58) 



In the following we present a calculation by Nakayama [134] which reproduces the 
effective quantum Hamiltonian of CDT, i.e. ( |6.58| from a continuum calculation. IVIore 
precisely, Nakayama considered the non-local "induced" action of two-dimensional 



See App. B.2 for more details. 
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quantum gravity in the proper time gauge. The general form of this action was first 
introduced by Polyal<ov |,64J and reads 

S[g] = I dtdx^ Q^^^Z^/^ + ^) ■ (6.59) 



where R is the scalar curvature corresponding to the metric g, t denotes time and x 
the spatial coordinate. In the proper time gauge the metric has the form 




This is in contrast to the conformal gauge with g = e'f'g. It was shown by Nal^ayama 
that in this gauge the classical dynamics is described entirely by the following effective 
action 

where 

/(t) = i y dx^, (6.62) 

and where k is an integration constant coming from the solution for the energy-momentum 
tensor component 7oi = 0. As can be seen from Eq. Leant = 7r/(t) is precisely 



the length of a spatial universe at constant t, as calculated from the metric ( |6.60| . 

One can now quantize the action for k = (m+ 1)^, where m g N can be interpreted 
as a winding number. It was further argued in [134J that in the quantum theory one 
should shift m by -1/2 leading to k = {m+^y. The classical Hamiltonian corresponding 



^1 - "a V ' 21 ' 

to the effective action | |6.61 with k = (m + |)2 then reads 



H„ = n,/n/+ (m+|)'y + A/, (6.63) 

where □/ is the canonical momentum conjugate to /. To canonically quantize the sys- 
tem we mal^e the standard replacement □/ -/|^, yieldingj^ 

H. = -/^-| + (m+i)^i+A/. (6.64) 

One observes that for m = this Hamiltonian precisely coincides with the effective 
quantum Hamiltonian for CDT, i.e. ( |6.58| . An interpretation for m>0 in the context of 



CDT was given in | fT26| 



*For completeness the spectrum of this Hamiltonian is also derived in App. B.2 
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Quite remarkably, when computing the propagator, 

GjC\U.L2:T) = {L2\e-^^-\U). (6.65) 

one observes that 

V(-l)'"(2m + 1) / dTG["'\U, L2; T) = _v^e-^(^i+^=) (6.66) 

precisely coincides with the two-loop amplitude of Euclidean DT |134| . While a sum- 
mation over the winding number might seem natural, a complete understanding of this 
relation is still lacking. 

Before moving to the next section let us briefly mention that, since in CDT we are 
using a lattice regularization, the continuum variables are only defined up to a factor of 
proportionality which should be fixed by comparing to a continuum calculation. Recall 
that the length / appearing in Nakayama's Hamiltonian is not the physical length. Leant, 
but Lcont/T^- Hence, since both models are described by the same effective quantum 
Hamiltonian, it is natural to also define the physical length Leant = tt/. for CDT. We 
will use this assignment in Chap.[8]when discussing the emergence of a semiclassical 
background in the context of CDT. 



Physical observables: Comparing CDT and DT 

As in the case of DT one of the simplest observables is the disc function or Hartle- 
Hawking wave function which describes the amplitude of creation of a universe from 
nothing. The result is given by ^^5J}, i.e 

1/Va(X) = — (6.67) 

X + VK 



The corresponding Euclidean expression (5.28), here labeled with a superscript {eu) to 



distinguish from the analogous CDT expressions, is 



W^'"\X) = {X- ^^)^X + VK. (6.68) 

Another interesting quantity to calculate is the Hausdorff dimension dn of the quan- 
tum geometry formally defined as 

{V{T)) ~ 7^" for T ^00. (6.69) 
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We saw that in Euclidean DT we had a fractal dimension of c/h = 4 which was related 
to the non-canonical scaling of time. We already observed above that in CDT time 
scales canonically and we therefore expect the Hausdorff dimension to be dn = 2. In 
the following we will checl^ this by a simple calculation. 
From ( |6.53| we have 

Ga( 7") ~ e-'^^ for 7 ^ oo. (6.70) 

From this we obtain 

This equation reflects the fact that at large T the quantum geometry lool^s effectively 
one-dimensional. Further we see from this expression that the average spatial length 
of the quantum geometry at intermediate T behaves as 

(L) = - ^. (6.72) 

Hence for typical time scales of 7 ~ l/VK we have 

{V{T)) r^T^ for 7 ~ (6.73) 



which yields a Hausdorff dimension of dn = 2. In this sense the quantum geometry 
of CDT is much better behaved as a model of two-dimensional quantum gravity in 
contrary to DT which has a fractal dimension of dn = 4. 

To get a better picture of the quantum geometry of CDT it is instructive to \ooW at a 
typical space-time geometry (Fig. [6^. We see that in contrast to the Euclidean model 



there are no spatial topology changes. Further, as we have already seen above the 
average spatial length of the quantum geometry behaves as 

(L) ~ ± (6.74) 



A simple calculation as will be presented in Chap. |8] reveals that also the fluctuations 
scale as 

(AO = VOy^ - (L) - (6.75) 



Hence we see that the two-dimensional quantum geometry is purely governed by quan- 
tum fluctuations and there does not exist a sensible semiclassical background geome- 
try. We will see in Chap. [8] how this situation changes when one of the boundaries is 
taken to infinity and the geometry becomes non-compact. 
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Figure 6.2: A typical two-dimensional Lorentzian space-time. The compactified direc- 
tion shows the spatial hypersurfaces of length (L) and the vertical axis labels proper 
time T. Technically, the picture was generated by a IVIonte-Carlo simulation, where a 
total volume of A/ = 18816 triangles and a total proper time of t = 168 steps was used. 
Further, initial and final boundary has been identified. 



6.7 1 Summary and outlook to higher dimensions 

In this chapter we introduced two-dimensional Lorentzian quantum gravity defined 
through CDT 

The triangulations used in CDT have a global time-slicing without spatial topology 
changes. This is unlil^e DT where spatial topology changes are naturally present. The 
global notion of time enables us to define a gravitational Wick rotation. We solved the 
combinatorial problem and performed the continuum limit. As a result we saw that two- 
dimensional CDT and DT are distinct theories. Whereas the continuum dynamics of 
DT is completely dominated by spatial topology changes leading to a fractal dimension 
of dn = 4, the situation is different in CDT where we have dn = 2. In this sense two- 
dimensional CDT is arguably better suited as a model of two-dimensional quantum 
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gravity than DT. Before discussing several relations between both theories in the next 
chapter, let us first comment on the results for higher dimensions. 

In contrast to the disappointing situation for DT in higher dimensions, as described 
in the previous chapter, CDT lead to very interesting results for both d = 3 |135[ 11361 
[T371 [T381 [T391 [T40l MD as well as d = 4 HM HM [1451 EM MI\ (see also 
[211 [201 [T451 [T48| for a general overview). 

Let use briefly mention some of the numerical results obtained from IVIonte Carlo 
simulations in 3+1 dimensions. A very important non-trivial test for every non-perturbative 
formulation of quantum gravity is whether it can reproduce a sensible classical limit at 
macroscopic scales. The numerical results indicate that the scaling behavior of the 
spatial volume as a function of space-time volume is that of a four-dimensional uni- 
verse at large scales, a first indication of sensible classical behavior [142| . IVIoreover, 
after integrating out all dynamical variables apart from the spatial volume as a function 
of proper time, one can derive the scale factor whose dynamics is described by the 
simplest minisuperspace model used in quantum cosmology [143[ 11451 11461 II 47| . In 
Chap.|8]we will describe an analogous model in two dimensions where a semi-classical 
background emerges from quantum fluctuations when the geometries become non- 
compact. 

Having passed the first consistency checks regarding the macroscopical structure 
of space-time it is very interesting what predictions one can make for the quantum na- 
ture of the m/crastructure of space-time. One important observable which has been 
measured is the spectral dimension of space-time which is the dimension a diffusion 
process would feel on the space-time ensemble. Surprisingly, this quantity depends 
on the scale at which it is measured. IVIore precisely, one observes a dimensional re- 
duction from four at large scales to two at small scales within measurement accuracy 
| 144| . This gives an indication that non-perturbative quantum gravity defined through 
CDT provides an effective ultraviolet cut-off through a dynamical dimensional reduction 
of space-time and might therefore be non-perturbatively renormalizable. The origin of 
the renormalizability could be a nontrivial fixed point scenario as described by Wein- 
berg [16]. It is interesting to note that similar results were also obtained in the exact 
renormalization group flow method for Euclidean quantum gravity in the continuum 
ri49l [1501 [T5T1 [1521 II 53| . 

An interesting observation at this point is that such an effective ultraviolet cut-off 
through a dynamical dimensional reduction of space-time provides us with a non- 
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perturbative mechanism which regulates the theory. To come back to the discussion in 
Part I of this thesis, this means in particular that it is priori not necessary to introduce a 
fundamental cut-off such as the fundamental discreteness scale by hand. 
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CHAPTER 



Relating Euclidean and causal dynamical 
triangulations 



In the previous two chapters we introduced two-dimensional Euclidean quantum gravity 
defined through DT and two-dimensional Lorentzian quantum gravity defined through 
CDT. Even though both theories have some physically very distinct features such a 
the Hausdorff dimension it is nevertheless possible to establish certain relations be- 
tween them. In particular, following mi |T25J, we describe how CDT and DT can be 
related by respectively introducing or "integrating out" baby universes, i.e. spatial topol- 
ogy changes. 

Introducing spatial topology changes: From CDT to DT 

In this section we discuss the incorporation of spatial topology changes into the dis- 
cretized frameworl^ of CDT and show how this model relates to DT in the continuum 
limit. Although there exist several ways of implementing spatial topology changes on 
the discrete level, they all lead to the same continuum description fTg) . 

Before constructing the model let us first discuss some Lorentzian aspects of spatial 
topology changes. While in the Euclidean model the presence of baby universes was 
natural, their appearance in the Lorentzian model is far from obvious. Even though the 
metrics are W\cW rotated from Lorentzian to Euclidean signature, one in general still has 



93 



94 



Chapter 7. Relating Euclidean and causal dynamical triangulations 




Figure 7.1 : Illustration of the double light cone causal structure at a Morse point in an 
up-side down trousers geometry. 



to perform an inverse Wick rotation. However, this is not so easy for geometries with 
baby universes, since these geometries do not admit a Lorentzian metric everywhere, 
as we already mentioned above. This is related to the fact that it is not possible to 
find a non-vanishing vector field everywhere. However, since the geometry is compact 
there are only finitely many isolated points where this is not possible. These are so- 
called IVIorse points [1 54[ 11 55[ 11 56[ 11 57[ [T58]. If one imagines the splitting of a spatial 
universe into two, there is a IVIorse point precisely at the moment where the spatial 
topology is that of a "figure eight". Embedded in the two-dimensional geometry 
lool^s lil^e an up-side down trousers and the IVIorse point precisely corresponds to the 
saddle point. Contrary to a generic point on the manifold, such a IVIorse point does not 
have a unique timelike vector field perpendicular to its spatial slice. In fact it has two 
future, and two past light cones which is referred to as a double light cone structure 
[158J and the resulting geometry is called causally discontinuous. Such a double light 
cone structure is illustrated in Fig.[7T|for the case of the up-side down trousers. One 
observes that both legs of the trousers carry a future light cone belonging to the IVIorse 
point at the splitting, while the mother universe carries the two past light cones. 

Since there are indications that the Einstein-Hilbert action develops complex valued 
singularities at these Morse points [1 59], it is not a priori clear that the Wick rotation as 
employed in CDT is valid at these points. However, since we want to relate our model 
to Euclidean DT, we can confine ourselves to the Euclidean sector of CDT. In Chap. [9] 
we present a generalized model of CDT which includes spatial topology changes, and 
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Figure 7.2: Illustration of the one-step propagator with an off-splitting "baby universe". 



where one in principle would like to Wick rotate back to the Lorentzian sector. However, 
in this model there is a coupling constant associated to every splitting which might be 
used to absorb possible divergencies. 

Let us now, after this short digression, focus on the explicit construction of the 
discretized CDT model with spatial topology changes. The first step is to generalize the 
one-step propagator by allowing the initial loop to have topology of a "figure eight". One 
possibility to do so is to non-locally identify two points of an initial spatial universe with 
topology For a boundary of length /i this so-called pinching leads to a combinatorial 
factor of /i in the full one-step propagator. A baby universe is then created by assigning 
a disc function to one of the loops of the "figure eight" and the initial loop of a "bare" 
one-step propagator to the other one (see Fig. ^L2) . Putting the above relation into 
formulas, we see that the new, or "dressed", one-step propagator is related to the old, 
or "bare", one-step propagator as follows 

/i-i 

Ga(/i, h; t = l) = Gf (/i, h; t = l) + 2 5^ kw{k-l, g) (/, h; t = l), (7.1) 

1=1 

where w{l,g) denotes the so far undetermined disc function and the factor of two 
comes again from the two different ways of attaching the disc function. One observes 
the analogy to the peeling equation for the fixed geodesic distance two-loop amplitude 
in DT 

It is readily checked that the "dressed" one-step propagator obeys the usual semi- 
group property, i.e. ( |6.36| , 



Gxik, h; ti + t2) = J2 /; ti)Gx{i, h; fs). (7.2) 

/ 

IVIaking the standard choice ti = 1 and t2 = t - 1 to relate to the one-step propagator we 
obtain 

GxUi. h; t) = J2 ^^(/i' 1) ^^(/' '2; t-i), (7.3) 
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or in terms of generating functions 

G{x,y;g;t) = 

(x, z-\ l) + 2x|- [w{g. x)Gl'\x. z-\ 1) 



dz 



Ini z 



G{z,y;g;t-1). 



(7.4) 



Inserting the explicit solution for the "bare" one-step propagator G^''\x, z; g;l), i.e. 
( |6.38| ), as derived in the previous chapter, we get 



G{x,y;g; t) 



dw{g,x) d 
2x^— + 2xw{g.x)- 



g 



,y:g:t-l 



(7.5) 



l-gx \l-gx 

At this point neither the dressed disc amplitude w{g,x) nor the propagator G{x,y;g; t) 
are known. In particular, w{g,x) is not fixed to be the Euclidean disc function. In 
the following we will show how one can use scaling arguments to uniquely determine 
the continuum propagator G^{X, Y; T) and the continuum disc amplitude W/^{X) from 
this equation. Let us assume the standard scaling relations for the boundary and bulk 
cosmological constants as used in DT and the "bare" CDT model, 

1 



g 



-aX 



y 



-aY 



(7.6) 



Inserting the scaling relations into and setting t=0 we obtain the initial condition 

1 



G^[X,y■T = Q) 



X+Y 

where we introduced the following wave function renormalization 

G^{X,y■T) = aGx[x,y■t) 
or in terms of the lengths variables 

G^{U,L2■T) = a-^G^{k,l2,t). 



(7.7) 



(7.8) 



(7.9) 



Extracting the scaling of the disc function is more difficult, since it actually depends 
on the scaling of time. To do so Ambjorn and Loll used the following combinatorial 
identity 

'-^^Y.Y.'^^X'i-'S-.'V-M. (7.10) 



g 



which in terms of the generating functions reads 

dw{g,x) 



g 



dg 



E 



dz _i dw{g,z) 

G{x,z ;g;t) 



liri z 



dz 



(7.11) 
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X 



Figure 7.3: Decomposition of the disc function with a marked point in the bulk into a 
propagator with arbitrary time and a disc function with a mark on the boundary. 



Fig. |7T3| illustrates this relation: If we introduce a mark in the bulk by taking a derivative 
with respect to the fugacity g, the disc function can be decomposed into a propagator 
and a smaller disc function which has this mark on its boundary. Further, one sums 
over all t to ensure that the mark can be anywhere in the bulk. 

Recall from Chap.[5]and Chap.|6]that the disc function for DT scales as 

w^^^\g,x) = + a3/2l/l/^-)(X) (7.12) 

while the disc function for CDT scales as 

w^'\g.x) = a-'wi^'\X). (7.13) 

Hence, we make the following general ansatz for the scaling relation which includes 
both cases 

w{g,x) = Wns + 3'^W^{X) + less singular terms. (7.14) 
Further, the general scaling ansatz for the time variable reads 

T = a't, 6>0. (7.15) 

which describes both the non-canonical scaling of time for DT, i.e. e = 1/2 as well as 
the canonical scaling of time for CDT, i.e. e = 1. Below we show that by allowing 
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the branching into baby universes to contribute in the continuum limit one is led to a 
non-canonical scaling of time as for DT. 

Inserting the scaling relations ( |7.6| ), ( |7.14| | and ( |7.15| ) into the combinatorial identity 
( |7.11| we obtain 

dwn. _ _,91/1/a(X) 



dg 



23"-^- 



dA 



j dT j dZ Ga(X, -Z; T) 



dz 



Zr dZ 



(7.16) 



where we have (x, g) = (xc, gc) in the non-singular part. 

From the last equation and the requirement that e > it follows that there are only 
two consistent choices for r]-. 

• Scaling 1: r] <0 

In this range the non-scaling part does not survive the continuum limit and Eq. 
( |7.16| becomes 



2a^ 



j dT j dZ G^{X, -Z; T) j 



1 ^W^{Z) 



dZ 



(7.17) 



Hence, we conclude that e = 1 as in the case of CDT. We observe that whenever 
the non-scaling part of the disc function is not present in the continuum result 
time must scale canonically. 

Scaling 2: 1 <r] <2 

In this range Eq. ( |7.16| ) splits into two equations 

^^2 ^1/1/a(X) _ 1 dWr 



- a ' 



(9A 



2a^ dz 



and 



dWr, 



dg 



dT I dZ Ga(X,-Z; T), 
1 5M/a(Z) 



dT / dZ Ga(X, -Z; T) 



(7.18) 



(7.19) 



z, dZ 

From the first equation we see that a''^ = 1/a^ while from the second equation 
we have a'^^V^" = 1- Combining these requirements we get e = 1/2 and r] = 3/2 
which precisely correspond to the scaling relations of DT, i.e | |5.46| and | |5.27| . 
Further, we observe that using e = 1/2 and r] = 3/2 Eq. ( |7.18[ ) reads 

dW^{X) 



(9A 



dTG/\{X, L2 = 0, T). 



(7.20) 
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Figure 7.4: (a) At every point in the quantum geometry there is an infinitesimal baby 
universe, (b) When trying to draw the baby universes on a plane one is forced to draw 
smaller and smaller triangles reflecting the fractal structure of the geometry. 



which is nothing but the relation between the disc function and the propagator as 

1/2 and r] 



obtained in DT, i.e. (5.551. In addition fore 



j dT j dZ G^{X, -Z; T) 



dZ 



= 3/2 Eq. ( |7.19[ ) becomes 
const, (7.21) 



where the constant is related to the non-scaling part of the disc function and does 
not play any role in the continuum theory. 



The relation {7.20) has an intersting interpretation in terms of baby universes, namely, 
it shows that at any mark in the bulk there is a baby universe whose boundary is of 
length of the cut-off. Further, since the mark can be anywhere, we conclude that near 
every point of the quantum geometry there is a baby universe with cut-off size bound- 



ary This situation is illustrated in Fig. 7.4 (a). Mapping the picture into the plane, as 



shown in Fig. [Z4] (b), clearly visualizes the fractal structure of the quantum geome- 
try. Another indication for the importance of geometries with infinitesimal boundaries 
comes from the Laplace transform of (|7.21[), i.e. 



JdT JdL Ga(/., L'; T) U W^[U) 



6{L). 



(7.22) 



Essentially, this equation shows that the distribution of geometries is peaked around 
universes that have infinitesimal boundary length. 
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Using the above scaling relations one can now analyze the scaling limit of (L5^ to 
obtain a differential equation for the dressed propagator. In order for the scaling limit of 
the equation to exist, the critical values x^, gc and w^s must satisfy two relations which 
can be straightforwardly determined from {7.5} . The remaining continuum differential 
equation equation reads 

d 



dT 



Ga{X, Y; T) 



d_ 
dX 



-A)Ga(X, T) 
1/1/a(X)Ga(X, Y; T) 



(7.23) 



The first term on the right-hand-side of ( |7.23| ) is related to the bare CDT model, while 
the second term corresponds to the creation of baby universes as is also present in 
DT. In particular, we saw in the previous chapter that the first term was related to the 
effective quantum Hamiltonian of CDT and hence represents the kinetic term of the 
more general model ( |7.23[ ). 

In the case where rj < (more precisely 77 < 1) the last term in ( |7.23| ) diverges. 
This reflects the fact that the bare model is incompatible with spatial topology changes. 
However, dropping the divergent term, i.e. excluding spatial topology changes we have 
already seen that for r/ < we need to have that 6 = 1 leaving us with the continuum 
differential equation for the propagator of the bare CDT model, i.e. ( |6.41| |, as a unique 
solution 



(X2-A)Ga(X, V; T) 



(7.24) 



For the second possible scaling, where 1 < 77 < 2, the last term on the right- 



hand-side of ( |7.23[ ) will dominate over the first while still being compatible with the time 
derivative term. Hence, the kinetic term of the model does not survive the continuum 
limit and the dynamics is entirely governed by the off-spliting baby universes, or as 
Ambjorn and Loll express it, once we allow for the creation of baby universes, this 
process will completely dominate the continuum limit. We have seen above that the 
only scaling consistent with 1 < 1] < 2\s {i],6) = (3/2, 1/2). Inserting the corresponding 
scaling relation into ( |7.23[ ) we obtain the following differential equation 

1/1/a(X)Ga(X, V;7)l. (7.25) 



d d 
Ga(X, Y; T) = -2- 



This is precisely of the form of the differential equation ( |5.47| for the propagator of 
DT. However, at this stage W^{X) and with that Ga{X, Y; T) are still undetermined. In 
the following we want to show how this equation combined with Eq. ( |7.20| ), completely 
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determines the continuum disc function W^{X) to be that of DT. Integrating ( |7.25[ l with 
respect to T and using the initial condition 

G^{X,Y-T = 0) = ^^ (7.26) 

i.e. in Laplace transform language 

Ga(X,/.2 = 0; 7 = 0) = 1, (7.27) 

we obtain 



1/1/a(X)|^1/1/a(X) 



(7.28) 



-1 = A 

dx 

From dimensional analysis one can easily see that W^[X) must have mass dimension 
-3/2 and hence can be written as V\/l[X) = X^F{\//\/X). This condition implies the 
following general form for the solution of the disc function 



1/1/a(X) = a/-2AX + 62X3 + c2A3/2. (7.29) 

We can further fix one of the two constants by noting that W/\{X) Is not allowed to 
have any singularities or cuts for 5R(X)>0. This requirement Is essentially the same as 
saying that the inverse Laplace transform of | |7.29| l, I.e. 1/1/a(/-), should be bounded for 
oo. This completely fixes the analytic structure of the disc function, yielding 



Upon rescaling of the boundary and bulk cosmologlcal constants we can absorb the 
constant b and obtain the disc function of DT, I.e. ( |5.28| ), 



1/1/a(X) = {X- ^VK)\/x + VK. (7.31) 



Integrating out baby universes: From DT to CDT 

In the previous section we showed how to obtain DT from CDT by Introducing spatial 
topology changes, i.e. by allowing baby universes to appear. In the following we want 
to show the opposite, namely how to obtain CDT from DT by "integrating out" the 
baby universes [125J. Since the baby universes completely dominate the continuum 
dynamics of DT such a relationship has to be established on the discrete level. 
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We remind the reader of the following differential equation for the discrete fixed 
geodesic distance two-loop function in DT, i.e. ( |5.41| |, which we obtained from the peel- 
ing procedure, 

^^(^''Hz, w-g- 1) = -2^ [f(-)(g,z)G(^")(z, w-g: t)] , (7.32) 
where we defined 

f^^^\g.z) = -\{z-gz^) + w^^^\g,z). (7.33) 
The initial condition of this differential equation was given in ( |5.36| ), i.e. 

G{z,w-g-Q) = - (7.34) 

zw zw — 1 



Here the first term in ( |7.33| ) corresponded to the first decomposition move, e.g. re- 



moving a triangle, while the second term corresponded to the second move, where 
a double linl< was removed and a baby universe (disc function) splits off. A first idea 
for integrating out the baby universes would therefore be to drop the second term in 
( |7.33| . However, as we are using unrestricted triangulations (recall Fig. ( |5.3| (b)) also 



double linl<s are present. Hence, there could still be off-splitting branched polymers 
which correspond to zero volume disc functions, i.e. ( |5.14| |, 



(z) = M/(-)(^ = 0,z) = \[z- Vz^) . (7.35) 



Thus we integrate out the (non-zero volume) baby universes by simply replacing the 
term i/i/(^")(g^, z) in ( |7.33| ) by w'^^p\z), yielding 



^G(z, w- g- 1) = -2^ [%, z)G(z, w- g: t)] . (7.36) 

with 

f{g,z) = ~^{z-gz') + w^'p\z), (7.37) 

subject to the same initial conditions {7.34} . Fig. [73] illustrates the peeling equation 
(TM} . 

We now want to show how to obtain the continuum CDT propagator by talking a 
suitable continuum limit of ( |7.36| . It is natural to use the standard scaling relations for 
the boundary and bu\k cosmological constants which were applicable for both DT and 
CDT 

g = g^e-'''^, z = z,e'^, w = z,e'^, (7.38) 
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Figure 7.5: Illustration of the peeling relation | |7.36| ). In each step of the decomposition 
either a triangle is removed or a double link which results into an off-splitting of a 
branched polymer. 



As before the demanded continuum initial condition 

Ga(Z, 1/1/; 7 = 0) = (7.39) 
already fixed the wave function renormalization of the propagator, i.e. 

Ga(Z, 1/1/; T) = aG{z, w; g; t). (7.40) 

Further, we note that to obtain a non-trivial continuum dynamics we demand that 
f{g,z) ~ 0{3^). This requirement completely fixed the scaling relations: Firstly, form 
f{g, z) ~ 0{a^) and the scaling relations ( |7.38| it follows that 

7 -at (7.41) 

as it is the case for CDT. Secondly, the requirement that the order one and order a part 
of f{g,z) should vanish determines 

1 
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Inserting the complete scaling relations into ( |7.36| ) yields (under rescaling of the cou- 
plings) 

^G,(Z, 1/1/, T) = - A [(z2 _ a)Ga(Z, 1/1/, T)] (7.43) 

with the above initial conditions. This differential equation precisely coincides with the 
differential equation for the propagator of CDT, i.e. ( |6.41| |. Hence we showed how to 
obtain CDT from DT by integrating out baby universes. 



The renormalization relation 



As an interesting relation between DT and CDT, it was shown in [125 ] that it is pos- 
sible to "renormalize" the couplings and time of DT in such a way that the differential 
equation for the DT propagator is matched to the one of the CDT propagator. This 
"renormalization" relation can be established both in the discrete as well as in the con- 
tinuum [1 25j . We will in the following only derive the continuum relation. 

Let us for the remainder of this chapter label all Euclidean quantities with a tilde. 
The differential equation for the Euclidean propagator thus reads, 



df 



where 



dX 



W~^{X) = {X - -V X 



A 



(7.44) 



(7.45) 



was the Euclidean disc function. The corresponding equation in the framework of CDT 
reads 



m{X)G/,{X, Y; T) 



with 



W/^{X) = - A. 



(7.46) 



(7.47) 



Here the hat (not to be confused with the tilde) should remind ourself that 1/1/a(X) is not 
the disc function of CDT but rather related to the effective quantum Hamiltonian of the 
system. 

To relate ( |7.44| | and ( |7.46[ ) we make the following ansatz for the renormalization 
relations of the couplings and time 



X = X{X,A), f = bK^''^T, K = K{N), 



(7.48) 



7.3 The renormalization relation 
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where b Is a constant and the factor of A^/"^ in front of the T has been chosen for 
dimensional reasons. To map the differential equation ( |7.44| | for DT into the differential 
equation ( |7.46| ) for CDT using the relations ( |7.48| , we must have the following wave 
function renormalizations 



Ga(X, Y; T) 



9X 



~GKi^^{x{x,^),Y{Y,^)■,b^'/' t) 



and 



dX 



m{x) = 2b^'/'{^)w~,^^^{x{x,^)). 



(7.49) 



(7.50) 



Integrating the second relation gives 



2 6A1/4 



dX 



1/1/a(X) 



X 



i/i4(X) 



(7.51) 



Since we require that X ^ oo should imply X ^ oo the integration constant must be 
k = 0. Inserting ( |7.45| and ( |7.47[ ) into ( |7.51[ ) with k = and performing the integrals 
gives 



b , x + Va 



2 1 



log 



V 



§av^/x 



A 



|AV4-v/X 



(7.52) 



7 



To further simplify this expression we perform a rescaling of both X and vA to choose 

3 . 



A 



A. 



(7.53) 



Using this relation \7.53) becomes 



X 



A 



Together with 



T 



2\JX + 
3 AV4 



T 



A 



6AV4 



(7.54) 



(7.55) 



these equations determine the whole set of renormalization relations. 

These renormalization relations map the differential equation of the DT propagator 
( |7.44| to the corresponding differential equation of the CDT propagator ( |7.46[ ). Further, 
it shows that the process of "integrating out" baby universes which we discussed in the 
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previous section can be understood as a renormalization of the cosmological constants 
and the time variable combined with a dressing of the propagator. 

It is interesting to notice that the relation ( |7.54[ ) is similar to one encountered in 
regularized bosonic string theory in dimensions d >2 [1 6011851 II 61), where the world 
sheet degenerates into so-called branches polymer. The two-point function of these 
branched polymers is related to the ordinary two-point function of the free relativis- 
tic particle by chopping off (i.e. integrating out) the branches, just leaving for each 
branched polymer connecting two points in target space one path connecting the two 
points. The mass-parameter of the particle is then related to the correspondin g par am- 
eter in the partition function for the branched polymers as X/y/K to X/\/K in (7.54|. 



Discussion and outlook 

In this chapter we explained several relations between two-dimensional Euclidean 
quantum gravity defined through DT and two-dimensional Lorentzian quantum grav- 
ity defined through CDT 

It was shown that when introducing spatial topology changes into the frameworl< of 
CDT one arrives at a continuum differential equation for the gravitational propagator 
which only allows for two possible scalings. In the first scaling time scales canoni- 
cally and spatial topology changes, i.e. off-splitting baby universes are forbidden. This 
scaling leads back to the 'bare" model of CDT without spatial topology changes as 
described in the previous chapter. However, in the second scaling where time scales 
non-canonical spatial topology changes are possible and the model naturally leads to 
DT. Further, we observed how in this scaling the baby universes completely dominate 
the continuum limit, explaining the fractal structure of two-dimensional Euclidean quan- 
tum gravity. 



In the Sec.[7y2]of this chapter we showed the opposite relation, namely, when inte- 
grating out the baby universes in DT one recovers CDT. This was done by introducing 
a new peeling equation similar to the one used to calculated the fixed geodesic dis- 



tance two-loop function in Sec. 5.3.2[ where outgrowing baby universes are replaced 



by branched polymers, i.e. zero-volume baby universes. 

In the last section of this chapter we showed how the process of "integrating out" 
baby universes could be understood as a renormalization of the cosmological con- 
stants and the time variable combined with a dressing of the propagator. 
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In summary, we showed many relations between CDT and DT by respectively in- 
troducing or "integrating out" baby universes. An unsatisfactory situation at this point 
is that we were in a sense not able to "regularize" the spatial topology changes. Ei- 
ther baby universes are absent as in CDT or the average number of baby universes 
is infinite as in DT which completely dominates the continuum dynamics. In Chap. [9] 
we show how one can naturally introduce a coupling which regulates spatial topology 
changes, leading to a finite number of baby universes in the continuum. Further, this 
formulation lead to the discovery of many other relations between CDT and DT within 
the framework of matrix models as will be discussed in Part IV. 
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CHAPTER 



The emergence of background geometry in two 
dimensions 



In the previous chapters we studied both Euclidean as well as causal dynamical trian- 
gulations with compact geometries. In particular, we saw that the quantum geometry 
in both cases does not posses a semiclassical bacl<ground, but is rather dominated 
entirely by quantum fluctuations. In this chapter we show how this situation changes 
when one studies the case of non-compact space-times. Before moving to our prime 
interest of non-compact space-times in CDT [ |129[ 162] let us first review some results 
of the Euclidean model [163]. 

Two-dimensional Euclidean quantum gravity with non-compact 
space-time 

In this thesis we have so far only considered two-dimensional Euclidean and Lorentzian 
quantum gravity with compact space-time. The study of two-dimensional Euclidean 
quantum gravity with non-compact space-time was initiated by the Zamolodchil^ovs 
(ZZ) in the context of Liouville theory [164| . In particular, they showed how to use con- 
formal bootstrap and the cluster-decomposition properties to quantize Liouville theory 
on the pseudo-sphere, i.e. the Poincare disl<. 
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It was later shown by Martinec [75] and Seiberg et al. [1651 1166| how the work 
of the Zamolodchikovs fitted into the framework of non-critical string theory, where the 
ZZ-theory could be reinterpreted as special branes, now called ZZ-branes. Particularly, 
they found that the ZZ-brane of two-dimensional Euclidean gravity was associated with 
the zero of of the disc function. 



WIC''\X) = {X- IVA)\/xWa. (8.1) 



2 

In |163[ I167[ \^68 ] it was shown how this could be understood in terms of worldsheet 
geometry, i.e. from a two-dimensional quantum gravity point of view. Let us briefly 
describe the main idea. We can use the fixed geodesic distance two-loop function, i.e. 
( |5.48| , to compute the average length of the final boundary. 



Further, we see from this that the average length of the boundary of a disc with a 
marked point in the bulk at geodesic distance T to the boundary is given by 

lLiT))y-,= ^-^ (8.3) 

\ V IIY.T x{T)+ Y ^ ' 

with X{T) = X{T, X = oo). Viewing X{T) as a "running" boundary cosmological con- 
stant with scale T, we see that both geodesic distance and boundary length diverge 
simultaneously if we choose Y = -X{oo) = -VA/2 which was precisely determined by 
= l/l//\(X(oo)) (recall Eq. ( |5.49| )). Hence, the geodesic distance from a generic point 
on the disk to the boundary diverges and in this way effectively creates a non-compact 
space-time. 

In this chapter we show that the same phenomenon occurs in the framework of 
two-dimensional Lorentzian quantum gravity defined by CDT. 



8.2 1 The emergence of background geometry 



Recall from Sec. 6.4 that the propagator of CDT is given by 



Mx,.:r,.%^^^^, ,8.4) 
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where X{T, X) is the solution of the characteristic equation 

dX 



-{X'-^). X{O.X) = X. (8.5) 



The explicit solution to this equation was given by ( |6.46[ ) which can be rewritten as 



X(t,X) = VAcoth VA(t + to), X = VAcoth VAto- (8.6) 

As described above we can view X{ T) as a "running" boundary cosmological constant, 
with T being the scale. For X > -VA we have X{T) VA for T oo, i.e. VA can 
be seen as a "fixed point", or in other words a zero of the "/9-function" -{X^ - A) in Eq. 

Using the semi-group property of the marked propagator 

POO 

Ga(X, Y; 7i +T2)= / dL Ga(X, L; T,) G{L, Y, T^). (8.7) 



we can now calculate the expectation value of the length of the spatial slice at interme- 
diate proper time t e [0, 7]: 



(8.8) 



Since the classical action is trivial, there is no reason to expect {L{t)) to have a clas- 
sical limit. We have already seen in Sec. [63] and Sec. |6j6] that for a generic situation 



the system is purely governed by quantum fluctuations. Except for boundary effects 
these quantum fluctuations are determined by the ground state of the effective quan- 



tum Hamiltonian Hx as introduced in Sec. 6.5 More precisely, an explicit calculation 



using ( |8.8| ) shows that in the situation where X and Y are larger than vA and where 
T > 1/y/K the system has forgotten everything about the boundaries and the expecta- 
tion value of L{t) is, up to corrections of order e^^*^ or e ^v^cr 0^ determined by the 
ground state of the effective Hamiltonian Hx- 

We will now study the situation where a non-compact space-time is obtained as 
a limit of the compact space-time in the same manner as described in the previous 
subsection for the case of DT. Hence we want to take T ^ oo and at the same time 
also take the length of the outer boundary at time T to infinity. Following the same 
reasoning as above, we see that the only choice of the boundary cosmological constant 



' Note that here W^ = - Ms not the disc function of CDT, as it was for the analogous term in DT. 
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Y independent of 7 where the length (/.(7))x,y,r goes to infinity for 7 ^ oo is the "fixed 
point" Y =-y/]\, since 

(L(T)) - 1 dG,{X, Y-T) _ 1 

which differs from the analogous DT expression, i.e. ( |8.2[ ), only in the definition of 

X{T,X). 

Inserting Y=-\/K into (8^1 one obtains in the limit T ^ oo: 

{L{t))x = ^ sinh(2v^(t + to(X))), (8.10) 
V A 



where to{X) is inversely defined through Eq. p^ . 

At this point we want to remind the reader that starting from a lattice regularization 
and talking the continuum limit the length L is only determined up to a constant of 
proportionality which should be fixed by comparing with a continuum effective action. 
In Sec. [63] we made such a comparison with a continuum calculation by Nakayama 
and concluded that L has to be identified with Lcont/T^- We thus have 

U,„,{t) = n{L{t))x = ^ sinh(2v^(t+to(X))). (8.11) 

V A 



We now consider the line element of the classical surface where the intrinsic geometry 
is defined by proper time t and spatial length Lcont{t) at constant t 

^ + % dO^ = dt^ + -"^'(^^^i' + ^(Xm ,,2 ,8.12) 
Ati^ 4A 

Here t > and tQ{X) is a function of the boundary cosmological constant X at the 
boundary corresponding to t=0 (see (|8j6])). A remarkable observation about Eq. ( |8.12| | 
is that the surfaces for different boundary cosmological constants X can be viewed as 
sections of the same surface, namely the Poincare disk with curvature R = -8A, since 
t can be continued to t = -to. The whole Poincare disk is then obtained in the limit 
X ^ oo, where the initial boundary is shrunken to a point and to = 0. 



The classical effective action 



Even though the classical action of our theory is trivial, we discussed in Sec. [63] that 
the effective quantum Hamiltonian of CDT can be derived from the following classical 
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effective action 





A/(0 + T ■ (8-13) 



4/(0 

With respect to this it is interesting to see in how far the results of this chapter relate to 

the classical solutions of this effective action. Solving the equations of motion of | |8.13| | 
gives 

l(t)= ^ sinh2^t, K > elliptic case, (8.14) 

/(t) = ^ cosh2VAt, K < hyperbolic case, (8.15) 

/(t) = e2^^ = parabolic case. (8.16) 

One observes that all solutions correspond to cylinders with constant negative curva- 
ture -8A. In the elliptic case, there is a conical singularity at t = unless k = 1 for 
which the geometry is regular at t = and corresponds precisely to the Poincare disc 
as described by ( |8.12[ ). However, k = 1 is exactly the value for which Nakayama's 
Hamiltonian corresponds to the effective quantum Hamiltonian of CDT. Concluding, 
we see that upon the required identification Uontit) = 7r/(t) the classical solution of 
Nal^ayama's effective action coincides with ( |8.1 1[ ). 



Quantum fluctuations 

In the considerations above we fixed the outer boundary cosmological constant to a 
specific value, namely the "fixed point" of the running boundary cosmological con- 
stant, rather than fixing the outer length of the boundary to a specific value. Since 
the boundary length and the boundary cosmological constant are conjugate variables, 
one pays the price that the fluctuations of the boundary size are large. In particular, it 



can be checked from (8.9 1 that the fluctuations are of order of the average length of the 



boundary itself, as was also true for the compact case 

{L'{T))x.y:r - {L{m,y.,r = -^^^^^^ = {L{T))ly.,r. (8.17) 

The same relation also holds for DT as was shown by Ambjorn et al. [163] . These 
large fluctuations are not only present at the outer boundary as described by ( |8.17| |, 
but also att < T. From this point of view it is even more remarkable that the emergent 
semiclassical background has such a nice interpretation in terms of the Poincare disc. 
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In the following we want to analyze the situation, where we by hand fix the boundary 
lengths U and L2 instead of the boundary cosmological constants. This is done in the 
Hartle-Hawl<ing Euclidean path integral when the geometries [g] are fixed at the bound- 
aries [72]. Recall from the previous chapters that the boundary geometry is complete 
specified by its length and that amplitudes with fixed boundary cosmological constants 
can be related to amplitudes with fixed boundary lengths by a Laplace transformation. 
Let us in the following analyze the temporal correlations between two spatial slices at 
intermediate times t and t + A with 0<t<t + A< T. For simplicity we only consider 
the situation where the entrance loop is shrunl^en to a point, i.e. U = 0. The temporal 
correlations are determined by the so-called connected loop-loop correlator 



{L{t)L{t + A))g^ = {L{t + A)L{t)),,,r - {m{L{t + A)),,,^. 



(8.18) 



Using the semi-group property p77} and the length version of the propagator pAj one 
obtains 



A sinli vA7 

2L2 sinh^ v^t sinh v^(t+A) sinh v^( T-{t+A)) 
sinh^VAr 



(8.19) 



A 



In comparison the average length is given by 



_ 2 s'mh VKt s'mh y//\{T - 1) sinli^ ^At 

^ ^^'^'"^^7K sinh Vat ^ \\uh'^/KT' 



(8.20) 

For fixed L2 and 7 ^ 00, i.e. in the situation of compact geometries, we obtain 

{L{t)L{t + A))[l^ = ^ e-^^^ (1 - e-^^^)' (8.21) 

and 



1 



L2 



l-e 



(8.22) 



These equations precisely agree with the picture presented above, namely that except 
for small t, where we are close to the initial boundary we have {L{t))L2 = 1/VA- i-S. 
the dynamics is determined by the ground state of the effective quantum Hamiltonian. 
Further, the quantum fluctuations are 

1 

/A 



(8.23) 
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Thus the quantum geometry is entirely governed by quantum fluctuations as is illus- 



trated in Fig.|6.2[ The time correlation between two slices L{t) and L{t + A) is also dic- 
tated by the same scale l/VK. Particularly, we observe that temporal correlations be- 
tween "elementary" spatial elements of size l/VK, separated in time by A fall off expo- 
nentially as e^^^. This is precisely what one would expect from the two-dimensional 
Einstein-Hilbert action with boundary terms, i.e. 

S[g] = l^jj dxdty/g{x, t) + X ^ dk+Y ^ dh. (8.24) 

If we force T to be large and choose a generic Y at which (/-2(T)) is not large, the 
quantum geometry will look lil^e a thin tube which expect close to the boundaries is 
"classically" of zero spatial extension, but due to quantum fluctuations has average 
spatial length l/y/K. 

Let us now repeat the above analysis for the case of non-compact geometries. 
Recall that only when the boundary cosmological constant was chosen to Y = -y/K 
one obtains a non-compact geometry in the limit 7 ^ oo. To implement this "critical 
value" of y in a situation where L2 is fixed instead of Y we fix L2{T) to the average 
value js^g) for Y=-\^, i.e. 

Ur)H4r))„._^^,^-L_l_^ ,8.25) 



From ( |8.19[ ) and | |8.20| | one obtains in the limit r — > cx) the following expression for the 
average length at intermediate time t 

{L{t)) = ^ sinh2VAt (8.26) 
V A 



which agrees with (8.10) as expected. Further, the loop-loop correlator equates to 



{L{t + A)/.(t))W = I sinh^ = ((^(0) - - ""'^')) ■ (^-^^^ 

In contrast to the analogue expression for compact geometries, one observes that 
( |8.27| | is independent of A. In particular, using A = we see that the quantum fluctua- 



tions for t > 1/Va are given by 



{AL{t)f ^ {L'{t)) - {L{t)f ~ ^m)- (8.28) 
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As in the case of compact geometries we can give an interpretation of ( |8.28[ ) in terms of 



correlations between "elementary" spatial elements of size l/vA. IVIore precisely, we 
can view the curve of length L{t) as consisting of N{t) ^ \/AL{t) ^ e^^*^ independently 
fluctuating "elementary" spatial elements of size 1/VK, each with a fluctuation of size 
1/a/A, i.e. each corresponding to a section of compact geometry. Thus the fluctuation 
AL{t) of the total spatial slice L{t) are of order y/N{t)/VK, i.e. 



e-^'". (8.29) 



mt)) 



Thus the fluctuations of L{t) around {L{t)) are small for t > 1/vA. In the same way 
we can understand the independence of the loop-loop correlator of A as the combined 
result, where L{t + A) is growing exponentially in length with e^^^ and the correlation 
between "elementary" spatial elements of L{t) and L{t + A) is exponentially decreasing 
with e~^^^, causing a correlation constant in A. 



Discussion and outlook 

In this chapter we described the transition form compact to non-compact quantum 
geometries within the framework of two-dimensional CDT In particular, the resulting 
quantum geometry can be viewed as the Poincare disk dressed with quantum fluc- 
tuations. Remarkably, it was shown that the quantum fluctuations of L{t) are small 
compared to the average value of L{t). In contrast to the analogous construction in the 
context of DT, this enables us to view the average geometry as a true semiclassical 
background. 

The main idea underlying this construction is similar to the appearance of ZZ- 
branes when described from a worldsheet perspective, i.e. form a two-dimensional 



quantum gravity point of view [163] . In DT the non-compactness arose when the run- 
ning boundary cosmological constant X{T) approached its fixed point described by 



l/l/;S^"^(X(cx))) = 0, i.e. to X(oo) = VK/2 (see Jsj). In the case of CDT the transition 
occurs in exactly the same manner, namely when the running boundary cosmological 
constant X{T) goes to for 7 ^ oo which is determined by W/\{X{oo)) = 0. How- 
ever, in the case of CDT Wa{X) is not the disc function, but is related to the effective 
quantum Hamiltonian of the system. It is interesting to observe that both processes are 
essentially the same, since both fixed points can be related by the mapping between 
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the coupling constants of both theories as derived in the previous chapter, I.e. ( |7.54| , 



X 2VXWA 
A V 3 AV4 ^ ' 



where the Euclidean couplings are labeled with a tilde again. From this expression we 
can see that X ^ \/K corresponds precisely to X ^ Va/2. 

Another Interesting observation Is that the geometries described by ( |8.12| |, for dif- 
ferent initial boundary conditions, are all sections of the same surface, namely the 
Poincare disc. Hence, they obey the Euclidean Hartle-Hawking no-boundary condi- 
tion. This is particularly surprising, since before rotating to the Euclidean sector we 
initially started with a path integral over entirely Lorentzian geometries. It would be 
Interesting to investigate whether this also holds In higher dimensional CDT. So far, the 



computer simulations, as briefly described In Sec.|6.7[ seem to be compatible with this 
condition. 
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Part III 



Third quantization of 2D causal 
dynamical triangulations 



119 



CHAPTER 



A field theoretic perspective of spatial topology 
changes 



In the second part of this thesis we introduced causal dynamical triangulations from a 
two-dimensional quantum gravity point of view in which the two-dimensional geome- 
tries were interpreted as a toy model for four-dimensional space-time. Tailing the string 
theoretic point of view we see the surfaces as worksheets of propagating strings. Us- 
ing this interpretation one is naturally lead to the idea of formulating a third-quantization 
of the theory, a so-called string field theory, in which strings can be annihilated, and 
created from the vacuum. IVlodels of string field theory (SFT) have been intensively 
studied in the frameworl< of the Euclidean model I90l [1691 [921 [TTOl [941 [T7T| . In this 
part of the thesis we develop an analogous formulation also for the Lorentzian theory. 
While in this chapter we compute certain "Feynman diagrams" of the model [130], we 
will in the next chapter introduce the whole string field theoretic frameworl^. 

Taming spatial topology changes 

In Chap.|7]we described how to introduce spatial topology changes into CDT. We saw 
that when viewing them as a purely geometric process^ one obtains Euclidean quan- 

^By this we mean that each distinct geometry (distinct in the sense of Euclidean geometry) appears 
with equal weight in the sum over two-dimensional geometries. 
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Figure 9.1 : Feynman rules: (a) shows a vertex of (/j^-theory which is assigned a cou- 
pling g. {b) shows the analogous interaction term for a splitting string. Similar to 
(y9^-theory we assign a string coupling g to this interaction. 



tum gravity. In particular, we recall from Sec. [7J] that there are two possible scaling 
relations entering the differential equation for the propagator ( |7.23| ), i.e. 



df 



G^{X, Y; T) 



d 



-2a 



dX 

v-iJL 

dX 



(X2-A)Ga(X, T) 



1/1/a(X)Ga(X, Y; T) 



In the scaling described by 



w{g,x) 



3^0 



a" 1/1/a(X), 
7"/a^ e 



1. 



(9.1) 

(9.2) 
(9.3) 



the creation of baby universes is forbidden, leading bacl^ to the "bare" CDT model. 
Further, in the scaling with 



w{g,x) 



a^O 



w.s + a'^W^iX), rj 
T/a\ e = 1/2, 



3/2 



(9.4) 
(9.5) 



the baby universes completely dominate the continuum limit, leading to Euclidean 
quantum gravity described by DT. 

Tailing a string field theoretic perspective, however, one can view each spatial topol- 
ogy change as a vertex of a Feynman diagram of the corresponding string field theory 
(see Fig. |9J}. In this spirit it is natural to assign a string coupling g to this process. 
This is in analogy to the coupling g in front of the interaction term of zero-dimensional 
V9^-theory, described by the action 

1 



S 3 

3!^- 



(9.6) 
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Figure 9.2: In all four graphs, the geodesic distance from the final to the initial loop is 
given by T. Differentiating with respect to T leads to eq. ( |9.8| ). Shaded parts of graphs 
represent the full, g^s-dependent propagator and disc amplitude, and non-shaded parts 
the CDT propagator. 

While in the Euclidean model such a coupling is absenj^we will see in the following 
that in the case of CDT it can be used to regulate the number of baby universes. From 
renormalization arguments one expects the coupling g also to scale, i.e. to be a non- 
constant function g = g(a) of the cut-off a. A geometric interpretation of this assignment 
will be given in the discussion at the end of this chapter. Since we are interested in 
a theory which smoothly recovers CDT in the limit as g ^ 0, it is natural to assume 
that f] = -1, as in CDT. Consequently, the only way to obtain a non-trivial consistent 
equation is to assume that g scales to zero with the cut-off a according to 



g = gs3 



(9.7) 



where gs is a coupling constant of mass dimension three, which is kept constant when 
a ^ [130J. With this choice, ([9J} is turned into 



_d_ 
df 



_d 
dX L 



(9.8) 



The graphical representation of ( |9.8| ) is shown in Fig. [9^ Differentiating the integral 
equation corresponding to this figure with respect to the time T one obtains ( |9.8| ). The 
disc amplitude 1/1/a,^,(X) is at this stage unknown. 

Note that one could in principle have considered an a priori more general branching 
process, where more than one baby universe is allowed to sprout at any given time 
step T. However, one observes from the scaling relation ( |9.7| that the corresponding 

^This is also reflected in the Hamiltonian of the string field theory l90l . 
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Figure 9.3: Graphical illustration of eq. ( |9.10| ). Shaded parts represent the full disc 
amplitude, unshaded parts the CDT disc amplitude and the CDT propagator. 



extra terms in relation ( |9.1| ) would be suppressed by higher orders of a and therefore 
play no role in the continuum limit. 

In the next section we show that the quantum geometry, in the sense defined above, 
together with the requirement of recovering standard CDT in the limit as 0, 
uniquely determines the disc amplitude and thus Ga,^,(X, Y\ T). 



9.2 1 The disc amplitude 



The disc amplitude of CDT was calculated in Sec. [M] i.e. ( |6.50[ ). In ( |6.48| it was de- 
termined directly by integrating G^{Ll, L2 = 0; 7) over all times. This decomposition is 
unique, since by assumption 7 is a global time and no baby universes can be created. 



In Sec.[Z2]it was shown that it could also be obtained from Euclidean quantum gravity 
(DT) by peeling off baby universes in a systematic way. By either method we found 

1 



(9.9) 



x + VK 

for the disc amplitude as function of the boundary cosmological constant X. In the 
present, generalized case we allow for baby universes, leading to a graphical repre- 



sentation of the decomposition of the disc amplitude as shown in Fig.[0| It translates 
into the equation 



+ gsj djj dUdL2 {U + /-2)Gi''i(X, U + L2- T)Wf,,,XLl)W^,,XL2) (9-1 



0) 



9.2 The disc amplitude 
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for the full propagator W^^gX^)^ where we have introduced a superscript (b) to indicate 
the "bare" CDT amplitudes, that is, 



<i(^) = l^A,,..o(X) = 1/1/a(X), 



(9.11) 



and similarly for G^^"] , quantities which were defined in eqs. (9.9i and (8.4 respectively. 
The integrations in | |9.10| can be performed, yielding 

1 



W/a,,.(X) 
Solving for W^gX"^) we find 



+ 



gs 



(9.12) 



W/a,,.(X) 



-(^^-A)+i3/a,,.(X) 

2gs 



(9.13) 



where we have defined 



K,s{X) = \j {X^ - + 4^s [gs WlgX^) + X - v^) . (9.14) 

The sign of the square root is fixed by the requirement that W^^gXX) ^ ^^{X) for 
gs 0, and 1/1/A,gs('^) is determined up to the value l/l/A,g^(VA). We will now show that 
this value is also determined by consistency requirements of the quantum geometry. If 
we insert the solution ( |9.13[ ) into ( |9.8| ) we obtain 

^Ga,,,(X, Y- T) = -^[m.JX) G,,JX, Y; J 

In analogy with ( |6.44[ ) and ( |6.45| ), this is solved by 

'^A,g.(X(7,X)) 1 



(9.15) 



Ga,^,(X, T) = 



l^A,^.W X{T,X)+Y' 



(9.16) 



where X(7, X) is the solution of the characteristic equation for | |9.15 



dX 
df 



= -14/a,^^(X), X{T = 0,X) = X, 



such that 



T 



dY 



(9.17) 



(9.18) 



Physically, we require that T can tal<e values from to oo, as opposed to just in a finite 
interval. From expression ( |9.18| for T this is only possible if the polynomial under the 
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square root in the defining equation ( |9.13| has a double zero, which fixes the function 

m,,xx) to 

^A.gs{X) = {X- C)^/iX+Cy-2gJC, (9.19) 



where 



C = UVK, U^ - U + 



gs 
A3/2 



0. 



(9.20) 



In order to have a physically acceptable l/l/A,g,(X), one has to choose the solution to 
the third-order equation which is closest to 1 . Quite remarkably, one can also derive 
|9.19[ ) from |9.13[ ) by demanding that the inverse Laplace transform W^^gXL) falls off 
exponentially for large L. In this region l/l/A,g^(X) equals wl{^^^{X) plus a convergent 
power series in the dimensionless coupling constant g^s/A^/^. 

One can check the consistency of the quantum geometry by using ( |9.16| ) in 

d 



OA 



poo /»oo 

^^A^) = dT dLG^,,XX, L- T) L W^,,XL) (9.21) 

JQ Jo 



which comes form the decomposition of a disc function with a mark in the bulk (see 
Fig.|7.3|), i.e. the continuum version of ( |7.10| |. Integrating ( |9.21| gives 

91/1/a,,,(X) _ 1/1/a,,.(X)-1/1/a,^,(C) 



w^,,XX) 

which is indeed satisfied by the solution ( |9.13| . 



(9.22) 



9.3 1 The loop-loop amplitude and the consistency condition 

We mentioned earlier that the propagator can be regarded as a building block for other, 
more conventional "observables" in two-dimensional quantum gravity. One of the most 
beautiful illustrations of this and at the same time a non-trivial example of what we 
have called quantum geometry is the calculation in two-dimensional Euclidean quan- 
tum gravity of the loop-loop amplitude from the fixed geodesic distance two-loop am- 
plitude [971. The full loop-loop amplitude is obtained by summing over all Euclidean 
two-dimensional geometries with two boundaries, without any particular restriction on 
the boundaries' mutual position. This amplitude was first calculated using matrix model 
techniques (for cylinder topology) |172| . 

To appreciate the underlying construction, consider a given geometry of cylindrical 
topology. Its two boundaries will be separated by a geodesic distance t, in the sense 
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Figure 9.4: Two different ways of decomposing the loop-loop amplitude into proper- 
time propagators and a disc amplitude. Two points touch in the disc amplitude 
1/1/, pinching the boundary to a figure-8, which combinatorially implies a substitution 
I^A,g,(/-) L 1/1/a,^,(/-) in the formulas. The time variables are related by T"i + T2 = T. 



of minimal distance of any point on the final loop to the initial loop. It follows that we 
can consider the geometry as composed of a cylinder where the entire final loop (i.e. 
each of its points) has a distance t from the initial one and a "cap" related to the disc 
amplitude, as illustrated in Fig. |9j4] (a). One can now obtain the loop-loop amplitude 
by integrating over all T and all gluings of the cap (we refer to [97J for details). An 
intriguing aspect of the construction is that the decomposition of a given geometry into 
cylinders and caps is not unique. One can choose another decomposition consisting of 
two cylinders of length Ti and T2, with 7i+ ^2 = T, joined by a cap, as illustrated in Fig. 



9.4 (b). As shown in [97J, the end result is indeed independent of this decomposition. 



The whole construction can be repeated for our new, generalized CDT model, in 
this way defining a loop-loop amplitude. IVIore precisely, although an exact equality of 
amplitudes corresponding to different decompositions lil^e those depicted in Fig.|9.4|(a) 
and (b) is not immediately obvious at the level of the triangulations of the discretized 
theor>j^ the continuum ansatz |9.23| ) below is self-consistent, in the sense that it leads 
to a non-trivial symmetric expression for the amplitude with a well-defined gs limit. 
The algebra is similar to that of [97]. 

We will denote the loop-loop amplitude by 1/1/a,^^(X, Y), and its Laplace transform by 



^because of the different arrangements of the proper-time slicings 
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W^^gXU, L2). The integral equation corresponding to Fig.|9^(a) is given by 

noo POO 

W^,gXLl, L2) = dT dL Ga,^.(/-i, L; T)LW;,,gXI- + ^2). (9.23) 
Jo Jo 

Laplace-transforming eq. ( |9.23| ), the integrals can be performed using eqs. ( |9.16| )- 
( |9.19| ). After some non-trivial algebra one obtains 



1 f[{X + C) + {Y+C)f 



f{X)f{Y)Ags V {f{X) + f{YW 
where we are using the notation 

f{x) = y/{x + cy - 2gjc = - Q- 

In the limit gs one finds 

l/l/i'i^(X, Y) ^ 



(9.24) 



(9.25) 



(9.26) 



A.g^v' - ■ / 2 Va(x + Va)2( y + y/Ky ' 

a result which could of course also have been obtained directly from \9.23) using | |6.44| |, 
( |6.45| ) and ( |9.9| ). We note that the corresponding expression in the case of Euclidean 
two-dimensional quantum gravity is given by 



1/1//!"' (X, y) 



2h{X)h{Y){h{X) + h{Y)y 



h{X) 



X + VA 



(9.27) 



which can be obtained from expressions similar to ( |9.16| )- |9.19| ), only with W/^^g^X) 
replaced by the Euclidean disc function 



wjC"\X) = {X- VA/2) h{X). 



(9.28) 



We observe a structural similarity between ( |9.24| | and \9.28) , with the function f{X) 
having the same relation to W^^g,{X) as h{X) has to W^""\x). 

Let us now prove the consistency condition that the two composition of the two-loop 
as shown in Fig. [9^ lead to the same result. Decomposing as in Fig.[9j4](a) leads to 

POO 

1/1/a,^.(/-i, L2- T)= dL G^,g^{U, L; T) L W^,g^{L + L2). (9.29) 
Jo 

while the decompositions of Fig.|9j4](b) reads 

poo POO 

l/l/A,g.(/-i, L2- T)= dL dU Ga,^,(/.i, L- 7i) L W^,gXL + /-') /-2GA,g.(/-', /-2; T - T,). 
Jo Jo 

(9.30) 
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Thus we arrive at the consistency condition: 

o pco poo 

0=—J^ dLj^ dU G^,,XLi.L:T,)LW^,,XI-+L')l-2GA.gAL'.L2:T-T,). (9.31) 

This condition was checked for the Euclidean model In [90]. Remarkably, it Is also 
satisfied In the our case. After a Laplace transformation eq. ( |9.31| ) reads: 



d 



dXdY dT, J_,^,, 2ni 



/oo+c 



dZ 



Ga,^^(X,-Z; T,)W^,,XZ)G^,,Sy.-Z;T-T,). (9.32) 



Using the explicit form of GA,gXX, Y; T), eq. ( |9.16[ ), we can perform the Z-lntegratlon in 



eq. \9.32\ . Using |9.12[ | we can express W/^g^X) In terms of Wf^^^X) given by | |9.10 
and finally using |9.17| | we can differentiate after Ti. The result is0 







dXdV 



X^{T,,X)-Y^{T-T,,Y) 



(9.33) 



which is satisfied. 

The existence of well-defined, symmetric expressions for the unrestricted loop-loop 
amplitudes In our generalized CDT model (at genus 0) and thus in standard two- 
dimensional CDT, formulas ( |9.24| | and ( |9.26| ), gives strong support to the claims that 
(i) the proper-time propagator does Indeed encode the complete information on the 
quantum-gravitational system, and (il) following the arguments given in [971 concerning 
the decomposition Invarlance of the loop-loop amplitude (c.f. Fig.|9.4||, the continuum 
theory Is diffeomorphism-invariant. 



Discussion and outlook 



The generalized CDT model of two-dimensional quantum gravity we have defined in 
this chapter is a perturbative deformation of the original model in the sense that It has 



'^The corresponding equation in tine case of non-critical string tineory is 

92 



= 



dXdY 



X{Ti,X)-Y{T- Ti,Y) 



Again one obtains tine remarkable result that the amplitude 1A'|^"'(X, Y; T) is independent of the subdivi- 
sion of T = T1 + T2 as shown in Fig.|9.4[ In addition we have in the non-critical SFT setting the additional 



consistency test that dTWj^"'\x, Y; T) = wj^^"\x, Y), where wj^"'\x. Y) is the so-called universal 



loop-loop correlator calculated from matrix model [1721 II 731 , i.e. 9.27 . This was verified in [90 
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a convergent power expansion of the form 

oo 

1/1/a,,.(X) = ^c„(X,A)(^)" (9.34) 



n=0 



in the dimensionless coupling constant gs/^^^^ [130 |. This implies in particular that the 
average number (n) of "causality violations" in a two-dimensional universe described 
by this model is finite, a property already observed in previous two-dimensional models 
with topology change [1741 11751 11281 1176] . The expectation value of the number n of 
branchings can be computed according to 

(„> . (9.35, 

which is finite as long as we are in the range of convergence of W/^g^{X). As already 
mentioned, this coincides precisely with the range where the function Wx^g^X) behaves 
in a physically acceptable way, namely, W^^g^L) goes to zero exponentially in terms of 
the length L of the boundary loop. The same is true for the other functions considered, 
namely, G^^g^U, L2; T) and GA^g^U, L2). 

The behaviour | |9.35| should be contrasted with that in two-dimensional Euclidean 
quantum gravity, and is reflected in the different scaling behaviours ( |9.3| ) and ( [931 ) for 
the time t. These scaling relations show that the effective continuum "time unit" in 
Euclidean quantum gravity is much longer than in CDT, giving rise to infinitely many 
causality violations for a typical space-time history which appears in the path integral 
when the cut-off a is fallen to zero. This phenomenon was discovered in the seminal 
paper I95l. 

As we will see in the next chapter, the calculations presented here should be seen 
as pertaining to the genus-0 sector of a generalized CDT model, which also includes a 
sum over space-time topologies. Although we have not given a precise definition of the 
higher-genus amplitudes in this chapter, one would expect them to be finite order by 
order. If the handles are as scarce as are the baby universes in the genus-0 amplitudes, 
it might even be that the sum over all genera is uniquely defined. Whether or not this 
is so will clearly also depend on the combinatorics of allowed handle configurations. 

In the context of higher-genus amplitudes, it is natural to associate each handle 
with a "string coupling constant", because one may think of it as a process where (one- 
dimensional) space splits and joins again, albeit as a function of an intrinsic proper 
time, rather than the time of any embedding space. An explicit calculation reveals that 
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in the generalized CDT model this process is related with a coupling constant gl (see 
next chapter), which one may think of as two separate factors of gs, associated with the 
splitting and joining respectively. 

How does the disc amplitude fit into this picture? From a purely Euclidean point 



of view all graphs appearing in Fig. [9^ have the fixed topology of a disc. However, 
from a Lorentzian point of view, which comes with a notion of time, it is clear that the 
branching of a baby universe is associated with a change of the spatial topology, a 
singular process in a Lorentzian space-time [159]. One way of keeping track of this 
in a Wick-rotated, Euclidean picture is as follows. Since each time a baby universe 
branches off it also has to end somewhere, we may think of marking the resulting "tip" 
with a puncture. From a gravitational viewpoint, each new puncture corresponds to a 
topology change and receives a weight l/G/v, where Gn is Newton's constant, because 
it will lead to a change by precisely this amount in the two-dimensional (Euclidean) 
Einstein-Hilbert action 

Seh = / ^'^v^^- (9-36) 

ZuGn J 



Identifying the dimensionless coupling constant in eq. (9.1 1 with g{a) = e^/^"'''', one 



can introduce a renorma//zec/ gravitational coupling constant by 



^ ^ ^InAa^. (9.37) 



This implies that the foare gravitational coupling constant G/v(a) goes to zero like 1/| In 
when the cut-off vanishes, a ^ 0, in such a way that the product q^I^n" /K^l'^ is inde- 
pendent of the cut-off a. We can now identify 

e-i/Gr = g^/^^l^ (9.38) 

as the genuine coupling parameter in which we expand. 

This renormalization of the gravitational (or string) coupling constant is reminiscent 
of the famous double-scaling limit in non-critical string theor^j^ In that case one also 
has gs oc e"^/^™", the only difference being that relation ( |9.37[ ) is changed to 

= ^J^ + -\nNa', (9.39) 



^It is called the double-scaling limit since from the point of view of the discretized theory it involves a 
simultanous renormalization of the cosmological constant A and the gravitational coupling constant Gn- 
In this article we have already performed the renormalization of the cosmological constant. For details 
on this in the context of CDT we refer to [1 9J. 
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whence the partition function of non-critical string theory appears precisely as a func- 
tion of the dimensionless coupling constant gs/^^^"^- 



CHAPTER 



A causal string field theory 



In the previous chapter we introduced a generalization of CDT by incorporating spatial 
topology changes regularized by a coupling gs. Using certain Feynman rules for prop- 
agation and splitting of strings (spatial universes) we were able to solve the genus zero 
disc function to all orders in gs. In this chapter we embed these results in a broader 
framework of a causal string field theory [131| . 



10.1 I The string field theory framework 

In quantum field theory particles can be created and annihilated if the process does 
not violate any conservation law of the theory. In string field theories one operates in 
the same way with operators which can create and annihilate strings. From the two- 
dimensional quantum gravity point of view we thus have a third-quantization of gravity: 
one-dimensional universes can be created and destroyed. In [901 iT69l I92l IT701 [941 fT 71 1 
such a formalism was developed for non-critical strings (or two-dimensional Euclidean 
quantum gravity). We will follow the formalism developed there closely and develop 
a string field theory or third quantization for CDT which will allow us in principle to 
calculate any amplitude involving creation and annihilation of universes [131 J. 

The starting point is the assumption of a vacuum from which universes can be 
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created. We denote this state |0) and define creation and annihilation operators: 

[^{L),¥{L')] = L6{L- U), vl/(/.)|0) = (0|vl/t(/.) = 0. (10.1) 

This assignment corresponds to working with spatial universes where a point has 
been marked. This avoids putting in certain combinatorial factors by hand when glu- 
ing together universes. The operators ^{L) and U't(/.) will be assigned dimensions 
dim [v|/] = dim [v|/t] = 0. 

We could alternatively have chosen creation and annihilation operators which cre- 
ate and annihilate universes without such a mark. Then instead of ( |10.1| one would 
have 

[^H{L),¥{L')] = L-'6{L- U), vl/(/.)|0) = (0|vl/t(/.) = 0, (10.2) 

and the dimensional assignment would be dim = l and dim = 1. One could 
even let create marked universes and annihilate unmarked universes if one just 
compensates for missing combinatorial factors by hand. Here we will use the assign- 
ment ( |10.1| ). 

Let us write the differential equation for the finite time propagator derived in ( |6.41| | 
using the boundary length rather than the boundary cosmological constant as variable 

^Gir\U,L,; T) = u[^^-A)g[-\U,L,; T), (10.3) 

For convenience we have in ( |10.3| ) also marked the exit-loop L2 in order to have symme- 
try between loops at the initial time and the loop at the final time, i.e. Gj(^\Li, L2, T) = 
/.2Ga(/-i, L2; T), where G^{Li, L2, T) \Nas given by ( |6.47| |. 
We can now also write 

Gi'^\U, L2, T) = (L2|e-™°(^)Ki), Ho{L) = + M, (10.4) 

where Ho{L) was derived in ([6j56]lQ Associated with the spatial universe we have 
a Hilbert space on the positive half-line, and a corresponding scalar product making 
Ho{L) self-adjoint (see App. |B.2| for more details) 

(^i|^2) = I ^ri{L)ML)- (10.5) 



^ We changed the notation from Hl in | |6.56| to HoiL) in | |10.4| to indicate that it is the first-quantized 
Hamiltonian. 
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The introduction of the operators ^{L) and in |10.1| ) can be thought of as 
analogous to the standard second quantization in many-body theory. The single par- 
ticle Hamiltonian becomes in our case the "single universe" Hamiltonian Hq. It has 
eigenfunctions tpn{L) with corresponding eigenvalues E„ = 2ny/K, n = 1,2, 

,Pn{L) = Le-^'p,^,{L), Ho{L)iJn{L) = E,iJn{L), (10.6) 

where Pn-i{L) is a polynomial of order n-l|^ We now introduce creation and annihilation 
operators al and corresponding to these states, acting on the Focl^-vacuum |0) and 
satisfying [an, a]„] = 6n,nT We define 

vl/(/.) = J2^nML), M/t(/.) = J2^lrn{L), (10.7) 

n n 

and from the orthonormality of the eigenfunctions with respect to the measure dL/L we 
recover ( |10.1| ). The "second-quantized" Hamiltonian is 

/CO Jl 
-¥{L)Ho{L)^\J{L), (10.8) 

and the propagator Gj("\Li, L2; T) is now obtained as 

GjC\U,L2; T) = (0|vl/(/.2)e-^^°vl/t(/.J|0). (10.9) 

While this is trivial, the advantage of the formalism is that it automatically tal^es care 
of symmetry factors (lil<e in the many-body applications in statistical field theory) both 
when many spatial universes are at play and when they are joining and splitting. We 
can follow [90] and define the following Hamiltonian, describing the interaction between 
spatial universes: 

H = Ho - gs j dU j dL2^^{U)WJ^{L2)^{U + L2) 

-agsj dU j dL2¥{U + L2)^{L2)^{U)- j ^p{L)Wf{L), (10.10) 

where the different terms of the Hamiltonian are illustrated in Fig. |10.1[ Here gs is the 
coupling constant we have already encountered in Chap. |9] of mass dimension 3. The 
factor a is just inserted to be able the identify the action of the two g^^-terms in ( |10.10| 



2 



The precise form of Pn-i{L) is not of relevance here, but it can be readily obtained from the general 



formulas in App. 



B.2 
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(a) 




L1+L2 




P{L) 



Li + L2 Li L2 L 

(6) (c) id) 



Figure 10.1 : The elementary terms of the string field theory Hamiltonian ( |10.10| |: (a) 
the single spatial universe propagator, (b) the term corresponding to splitting into two 
spatial universes, (c) the term corresponding to the merging of two spatial universes 
and (d) the tadpole term. 



when expanding in powers of gs. We will thinly of a = 1 unless explicitly stated differ- 



ently. Note that the sign of all the interaction terms in ( |10.10[ ) is negative. This reflects 
that we want these terms to represent the insertion of new geometric structures com- 
pared to the "free" propagation generated by Ho- These structures should thus appear 
with positive weight when we expand e^^. H is hermitian except for the presence of 
the tadpole term. It tells us that universes can vanish, but not be created from noth- 
ing. The meaning of the two interaction terms is as follows: the first term replaces a 
universe of length Li + L2 with two universes of length and L2. This is precisely the 
process shown in Fig. |9.2[ The second term represents the opposite process where 
two spatial universes merge into one, i.e. the time-reversed picture. The coupling con- 
stant gs seems to play the role of a string coupling constant: one factor gs for splitting 
spatial universes, one factor gs for merging spatial universes and thus a factor gf when 
the space-time topology changes. 

In a certain way the appearance of a tadpole term is more natural in the CDT frame- 
worl^ than in the original Euclidean frameworl^ in [90]. Recall the discussion in Chap. [9} 
In a Lorentzian setting there is no regular disc geometry, it has to end in a "puncture". 
The tadpole term is a formal realization of this process. Also recall that we could as- 
sociate this process with a gravitational coupling constant, in this way linking it to gs. It 
was done by observing that to each splitting off of a baby universe we have a puncture 
where the baby universe ends. Thus the coupling constant gs related to the splitting 
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of spatial universes could be identified with the vanishing of universes. This shift can 
be made explicit in our string field Hamiltonian H in |l0.10 ). In |l0.10 ) the coupling 



constant gs is associated with splitting and joining of spatial universes. No coupling 
constant is associated with the tadpole term, i.e. the vanishing of a spatial universe. 
However we can redefine and MJ'': 

\\} = gs\\}, vi>t = lvl;t. (10.11) 

gs 

With this substitution the coupling constant gs is shifted from the splitting term to the 
tadpole term, i.e. precisely the shift mentioned above. In addition the term associated 
with the joining of spatial universes will have the coupling constant gs, which matches 
a change in topology. 

Finally, let us identify the real coupling constant appearing in ( |10.10| ): let us mea- 



sure everything in terms of 1/vA which is the natural length scale of our universe. 
Introducing the dimensionless length variable I = L\/K, the dimensionless boundary 
cosmological constant X = X/VA, the dimensionless time variable f = tVA, the 
dimensionless tadpole density p{L) = p{L)/\/K, the dimensionless coupling constant 



gs = gs/l^^^^ (already introduced in eq. ( |9.34| |) and finally the dimensionless Hamilto- 
nian H = H/y/K, we can write 

Ho = VaHo, Ho = J ^MJ\i)HoCL)^{i). (10.12) 

where U'(I) = ^{L) and U'^(I) = Uft(/.) satisfy the same commutation relation as 
^{L), M/1^(/.) when expressed in terms of L, and H = Vah, where 

-agsj dU j dl2¥{U + l2)^(L2mL{)- j ^p(L)^(L). (10.13) 

From this equation it is clear that the real coupling constant in the theory is the dimen- 
sionless gs, precisely the "double scaling" coupling constant which already appeared in 



the calculation of W^gX'^) ^"^cl G^^gX^' T) (i-e. |9.34[ )). From the discussion above 
we also observe that the expansion parameter for topology change of space-time is 

g^ = ^sVA^ 

In principle we can now calculate the process where we start out with m spatial 
universes at time and end with n universes at time T, represented as 

G,,,XLi. .., L^: L[. ... C; T) = (0|vl/(/.;) ... vl/(C) e-^V(/.0 ... vl/t(/. J|0). (1 0.14) 
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10.2 1 The genus zero limit 

10.2.1 I The disc amplitude 

Let us consider the simplest amplitude: a single spatial universe which disappears in 
the vacuum. This is precisely the disc-amplitude considered in the previous chapter. 
The topology of space-time was not allowed to change in the calculation in Chap. [9| 
We can incorporate that in the SFT-picture by choosing a = in ( |10.10| |, i.e. the genus 
zero limit. The disc-amplitude can then be expressed as: 

Wa,M= lim W/,,g^{L,T)= lim (0|e-^^("=°)vl/t(/.)|0). (10.15) 

It will describe all possible ways in which a spatial loop can develop in time and disap- 
pear in the vacuum without changing the topology of space-time. Note that the tadpole 
term in ( |10.10| ) is needed in order that the amplitude ( |10.15[ ) is different from zero since 
the state \ L) = V'^(/-)|0) is orthogonal to the vacuum state |0). We note that if a = we 
have 

(0| e-^"("=°)vl/t(/.i) ■ ■ ■ ^\U)\0) = (10.16) 
(0| e-^^("=°)vl;t(/.,)|0)(0| e-^^("=°)vl/t(/.2)|0)(0| • ■ ■ |0)(0| e-™("=°)vl/t(/.^)|0), 

this factorization being a consequence of the fact that if we start out with m spatial 
universes there is no way they can merge at any time if a = (it is easy to prove 
( |10.16| using the algebra of the vU's). 



Following [90] we obtain an equation for W/\_gXL) by differentiating ( |10.15[ ) with re- 
spect to T and using H|0) = 0: 

0= lim ^1/1/a,^^(/., 7) = lim (0|e^^"("=°)[H(a = 0),vl/t(/.)]|0). (10.17) 

The commutator can readily be calculated and after a Laplace transformation eq. 
( |10.17^ reads 

^{{X'-A)W;,,,XX) + gsWlJX))=p{X), (10.18) 

where the last term on the left-hand-side of ( |10.18| ) is a consequence of the factoriza- 
tion ( I10.16K 

Eq. ^10.18^ has the generalized CDT solution ([9T9j-([9j20) discussed in Chap. [9] if 

p{X) = l, i.e. p{L) = 6{L), (10.19) 
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which is a reasonable physical requirement: the spatial universe can only vanish in the 
vacuum when the length of the universe goes to zero. 

10.2.2 I Inclusive amplitudes 

Above we have understood how to reproduce the generalized CDT disc amplitude 
W^ gX^) the connected amplitude arising in SFT in the limit a = 0. We now want to 
understand how to reproduce the proper-time propagator G^y^^X, Y\ T) in the context 
of SFT. We have an "entrance" loop at time 7 = and an "exit" loop at time T. However, 
the propagator Y\ T) allows baby-universes to split off and propagate further 

than time T if they only vanish into the vacuum eventually (see Fig. |9.2[ ). 



We can reproduce this result in the a = limit of SFT by introducing the "inclusive" 
Hamiltonian [90]. Since we are working in the « = limit, we only have universes 
branching, not merging, during the time evolution, and all the branching universes 
except one have to vanish in the vacuum. The branching process is dictated by the 
term 

gsj dU j dL2^\U)¥{L2)^{U + L2) (10.20) 

in the Hamiltonian H, eq. ( lO.IO i. Once the branching has occurred, only one of the 
two universes can connect to the exit loop at time T, the other universe has to continue 
until it eventually vanishes in the vacuum, a process which is allowed to occur a time 
later than T. This scenario is captured by the replacement 

vl/t(/.i)vl;t(/.2) ^ W^,,XLi)'^KL2) + ^\U)W^,,XL2) (10.21) 

in eq. ( |10.20| ). Thus we arrive at the following "inclusive Hamiltonian" 

Hind = j ^-^^\L)H^[L)^[L) - 2gsj dU j dU 1/1/a,^,(/-i)v1/^(/-2)M/(/.i + L^), (10.22) 

and we obtain the following representation of G^y^X^i, L2; T) 

G[";1{U,L2- T) = (0|vl/(/.2)e-^^'-vl/t(/.,)|0). (10.23) 



Differentiating eq. |10.23| ) with respect to T, commuting Hind to the right to |0) (and 
using H,v,c/|0) = 0) one obtains after a Laplace transformation and removal of the mark 
on the final boundary ( |9.8| ). Thus the generalized CDT proper-time propagator also 
has a simple SFT description. 
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10.3 I Dyson-Schwinger equations 

The disc amplitude is one of a set of functions for which it is possible to derive Dyson- 
Schwinger equations (DSE). Here we consider a more general class of functions. De- 
fine the generating function: 

Z(J; T) = (0|e-^^ ciLj(L)^HL)^Qy (10.24) 

We have 



(10.25) 

j=o 



We have already seen that in the case where the coupling constant a = we had 
factorization: 

Z(J, T;a = 0) = e/^^-^(^)'^A..a,T) (^q ^S) 

where W/^g^L, T) denotes the disk amplitude where the universe decays into the vac- 
uum before or at time T, and where Wf, g^{L, T = oo) was the disc amplitude we already 
calculated. 

Following [90] we can obtain the DSE in the same way as for the disc amplitude, 
the only difference being that when the constant a is no longer zero these equations 
do not close but connect various amplitudes of more complicated topology. However, 
the equations can be solved iteratively. We denote 

Z{J) = lim Z(J; T), (10.27) 

T— >oo 

Z(J) being the generating functional for universes that disappear in the vacuum. We 
now have 

0= lim A(o|e-^"e/^^^(^)^'(^)|0)= lim (Ole"^^ H e/^^-'W^'^^^lO). (10.28) 

Commuting the U'(/.)'s in H past the source term effectively replaces these operators 
by LJ{L), after which they can be moved to the left of any U't(/.) and outside (0|. After 
that the remaining U't(/.)'s in H can be replaced by 6/6J{L) and also moved outside (0|, 
leaving us with a integro-differential operator acting on Z(J): 

0= r dLJ{L)o(L,J,jj)z{J) (10.29) 
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where 



LJ, 



6J 



HoiD- 



S{L) 



'5J{L) 

'^''-Jo 5J{L')5J{L-L') 



(10.30) 



Z{J, T) is a generating functional which also includes totally disconnected uni- 
verses which never "interact" with each other. It is of more interest to restrict ourselves 
to the study of "connected universes", i.e. universes where space-time is connected. 
The generating functional for connected universes is obtained in the standard way from 
field theory by taking the logarithm of Z(J, 7). Thus we write: 



F(J, 7) = logZ(J, 7), 



and we have 



TH 



vl/t(/.i)---vl/t(/.„)|0), 



5"F{J, T) 



(10.31) 



(10.32) 



and we can readily transfer the DSE ( |10.29| )-( [T"0.30| ) into an equation for the connected 
functional 

F[J) = lim F(J, T). (10.33) 

From ( |10.29| -( [Ta30| we obtain 

-gL t dU^I^^^^-agL r dUUm^^^\ (10 34) 



From ( |10.34[ | one obtains the DSE by differentiating | |10.34[ ) after J{l) a number of times 
and then taking J{L) = 0. 
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Let us introduce the notation 
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as well as the Laplace transform 

poo POO 

Xn)= dU-- dL„e-^'^'--^"'-"W{U Ln) (10.36) 

Jo Jo 

for the higher loop amplitudes. 

Let us differentiate eq. ( |10.34[ | after J{L) one, two and three times, then tal<e J{L) = 



and Laplace transform the obtained equations. We obtain the following three equa- 
tions (where Ho{X)f{X) = dx[{X^ - A)f(X)]): 

= Ho{X)W{X)-l + gsdx(^W{X,X)+W{X)W{X)y (10.37) 

= {Ho{X) + Ho{Y))W{X.Y)+gsdxW{X.X.Y) + gsdYW{X.Y.Y) 
+2g,{dx[W{X)W{X, Y)]+dY[W{Y)W{X, Y)]) 

+2...*a.(^(Mr)) (10.38) 

= {Ho{X) + Ho{Y) + Ho{Z))W{X.Y.Z) 

+gsdxW{X,X, Y,Z) + gsdYW{X, Y, Y,Z) + gsdzW{X, Y,Z,Z) 
+2gsdx[W{X)W{X, Y, Z)] + 2gsdY[W{Y)W{X, Y, Z)] 
+2g,dz[W{Z)W{X, Y,Z)] + 2g,dx[W{X, Y)W{X,Z)] 
+2gsdY[W{X, Y) W{ Y, Z)] + 2gsdz[W{X, Z) W{ Y, Z)] 

^^^^ mX.YyWjX.Z) -^ ,10.39) 

The structure of the DSE should now be clear. 

We can solve the DSE iteratively. For this purpose let us introduce the expansion 
of W{Xi Xn)\n terms of the coupling constants gs and a: 

OD OO 

W{X, X,)= 5^ 5^ g^W{X, X,;m,k). (10.40) 

k—n—l m—k—1 

The amplitude W{Xi Xn) starts with the power {agsY~^ since we have to perform 

n mergings during the time evolution in order to create a connected geometry if we 
begin with n separated spatial loops. Thus one can find the lowest order contribution 
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to W{Xi) from (|10.37[|, use that to find the lowest order contribution to W{Xi,X2) from 



( |10.38| ) and use this again in | |10.39| ) which involves X2, X3), etc. Returning to 
eq. | |10.37[ ) we can use the lowest order expression for W{Xi, X2) to find the next order 



correction to W{Xi), use this and the lowest order correction for W{Xi, X2, X^) to find 
the next order correction to W{Xi, X2), etc. 

Two remarl<s are in order: firstly, the integration constants which come by integrating 
( |10.37| |-( [T"0.39[ ) and the corresponding higher order equations are uniquely fixed by 



the requirement that the correlation functions fall off for the lengths /., ^ 00, i.e. the 

requirement that the Laplace transformed amplitude W{Xi X^) is analytic for X, > 0. 

Secondly, the expressions obtained for W{Xi X^) can of course be obtained directly 



from a diagrammatic expansion, using the interaction rules shown in Fig. |10.1[ the 
propagation being defined by Ho, and then integrating in a suitable way over the times 
Ti involved. 

Let us just list the first few orders: 

W(X;0,0) = (10.41) 

x + Va 

W(X;1,0)= ^ + 3^ . (10.42) 
W(X,Y;1,1) = ^ ^ (10.43) 

2Va{x + VKy{Y + VKy 

WiX. Y. Z; 2, 2) ^ ^^'^ + + Y ^ Z) ^ sVKjXY ^ XZ ^ YZ) ^ XYZ 

4Ai(VA + X)3(v/A+ V)3(v/A + Z)3 

These amplitudes involve no change in space-time topology. The genus one and genus 
two amplitudes to lowest order in gs, i.e. without any additional baby universes, become, 
using the results ( |10.41[ )-( [T"0.44| | in the iteration as described above. 



= 15Ai + llAX + 5^X^ + X3 

32Ai(v^ + X)5 

W(X;3,2) = fll319A^ + 19951A^X+21555A5X2+ 



2048At(VA+X)9 
16955A2x2+9765Aix^+3885AX^+945v^X^ + 105X^) . (10.46) 



In terms of diagrams this genus 2 amplitude corresponds to the following three dia- 
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grams (integrated suitably over the times 7,): 




W{X;3,2) 



As mentioned above the amplitude W{Xi X„) starts with the power {agsY^^ 

coming from merging the n disconnected spatial universes. The rest of the powers 
of ags will result in a topology change of the resulting, connected worldsheet. From a 
Euclidean point of view it is thus more appropriate to reorganize the series as follows 



W{X, X„) = {agsY^'Y.'^ag^yW.iX, X„) (10.47) 

h=Q 

oo 

l/l//,(Xi X„) = Y^d^i^i X,;n~l + 2h+j,n-l + h) (1 0.48) 

i=o 

and aim for a topological expansion in ag^, at each order solving for all possible baby- 
universe creations which at some point will vanish into the vacuum. Thus l/l//,(Xi, ... ,Xn) 
will be a function of gs although we do not write it explicitly. The DSE allow us to obtain 
the topological expansion iteratively, much the same way we already did as a power 
expansion in gs. 

Since we have W{X, Y) = 0{a) this term does not contribute to the lowest order 
and from the DSE | |10.37[ ) we obtain a closed equation for l/l/o(X): 



Ho(X) Wo{X; gs) + gsdx W^{X; gs) 



This equation is of course just eq. | |10.17| | and we have 

Wo{X) = 1/1/a,^,(X). 



(10.49) 



(10.50) 



Knowing Wo{X) allows us to obtain Wo{X, Y) from ( |10.38| ) since W{X, Y,Z) is of 
order 0{a^). Thus the 3-loop term does not contribute to the lowest a order of eq. 



(|10.38|l, which is 0{a), and we have to the lowest order: 

Ho{X)+2gdM{X) + Ho{Y)+2gsdYWo{Y))Wo{X,Y) 



-2d.d.('M§^) 



(10.51) 
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Thus Wo{X, Y) is entirely determined by the knowledge of Wq{X). We note that using 
the definition ( |9.12[ ) we can simplify ( |10.51| : 

J^(m/a,,,(X)1/1/o(X, V)) + ^(l^A,,.(V)l/l/o(X, Y) 



1 ( W^,,XX)-W,,,Xy) \ 

R.dXdYK x-Y ) ^^^■^'^> 



gsdXdY 

The solution Wo{X, Y) can readily be found from eq. ( |10.52[ ), yielding 

■ yry Y)- ^ ^ f [{X+C) + {Y+C)r \ 

where 

f{x) = ,/{x+Cy-2gs/C = m,gXX)/{X - C). (10.54) 

and 

C=UV]\, U^-U + ^ = 0. (10.55) 



as in Sec. [92 



In fact this solution was already found in Sec. |9.3| since we by definition have 

POO 

Wo{X,Y)= / dT W/,_gXX,Y;T), (10.56) 
Jo 

where W/^g^X, V; 7) is the Laplace transform of the loop-loop function W/^g^U, L2; T) 
defined in ( |9.30| l with Ti = 7/2. When expanded to lowest order in gs we reproduce 
( [1043) . 

It should now be clear that one can iterate the DSE in a systematic way as a power 
series in the number of handles of the world sheet exactly the same way as we iter- 
ated the DSE as a function of the coupling constant gs: As an instructive example we 
calculate the genus one amplitude Wi{X). We expand 

W{X) = Wo{X) + ag^W,{X) + a^g^W2{X) + --- (10.57) 
W{X,Y) = agsWo{X,Y) + a^g^W,{X,Y) + --- 

While ( [10.49[ ) was the 0th order in a of eq. ( [10.37| |, the 1st order reads 



d_ 
dX 



(i4/a,^,(X) 1/1/1 (X) + l/l/o(X,X)) = 0, (10.58) 
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where Wq{X,X) is given by ( |10.53| |. The integration constant is fixed by the require- 
ment that l/l/i(x) is analytic for X > 0, i.e. that Wi{L) falls of as /. ^ oo. We then 
obtain 



It is seen that if we expand Wi{X) in powers of gs, the first term will reproduce (10.45|) 



10.5 1 Discussion 



We have developed here a string-field theory based on the CDT quantization of two- 
dimensional quantum gravity. It shares many properties with the original non-critical 
SFT, and the whole formalism was borrowed from non-critical SFT. Yet, it is different 
and in some ways simpler. The tad-pole term is simpler. It is simply p{L) = 6{L), telling 
us that universes can only disappear in the vacuum if they have zero spatial "volume". 
This is in accordance with the interaction between spatial universes, which preserves 
the total length (the total spatial volume). In non-critical SFT the evolution in proper 
time results in a process where the original spatial universe at proper time 7=0 spawn 
an infinity of (infinitesimal) baby universe during the time evolution. This is linked to the 
fact that the proper time in non-critical SFT has the anomalous length dimension 1/2. In 
our new CDT-based SFT the situation is different. The proper time T has the canonical 
dimension 1 , and the number of baby universes created during the time evolution is 
finite (see | |935) ). 

As discussed in Chap. |9] it is not possible to connect the non-critical SFT and the 
CDT-based SFT by a simple analytic continuation in the coupling constant gs, not even 
in the a = limit. As already shown in Chap.|7| starting out with a discretized, regular- 
ized version of the theory, the Euclidean theory (quantum Liouville theory) is obtained 
if the "bare" dimensionless coupling constant g is of order one. However, the relation 
between the bare coupling constant and the dimensionful continuum coupling constant 



gs used here and in the previous chapter is, as mentioned in Sec. 9.1 



g = gs3' (10.60) 



As discussed in Sec. 9J]the generalized CDT continuum limit corresponds to gs fixed. 



a ^ 0, and thus to g{a) 0. The fact that g{a) goes to zero in CDT SFT is of course 
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related to the finite number of baby universes generated in this theory. In contrast we 
have an infinite number of baby universes generated in non-critical SFT where g is of 
order one. 

However, there is clearly a deep connection between the Euclidean theory and 
the CDT theory awaiting to be fully understood. In Sec. |7j2] it was shown that if one 
integrate out the "excessive spawn" of baby universes in Euclidean two-dimensional 
quantum gravity one recovers the CDT theory and the mapping between the dimen- 



sionless variables X/vA in the two theories was found, i.e. ( |7.54[ ). This mapping was 
later discovered by Seiberg et al. [165| as the uniformization map from the algebraic 
surface representing the "semiclassical" non-critical string to the complex plane. The 
singular points of the algebraic surface corresponds to ZZ-branes where there is a tran- 
sition from compact to non-compact topology [163[ I167[ 1168] . These singular points 
are mapped to points in the complex plane where one has a similar transition from 
compact to non-compact geometry in CDT context as discussed in Chap.|8} 

It is interesting to try to generalize the present CDT-based SFT to include the cou- 
pling to matter. In particular, the presence of a c = 1 barrier can then be addressed. 
Since the c = 1 barrier partly can be understood as the result of an excessive creation 
of baby universes, such that the two-dimensional worldsheet is torn apart i161[l177| it 
is clear that the CDT theory may behave differently at c = 1. Numerical simulations hint 
that there still might be a barrier for large c [1781 . but nothing is presently known with 
certainty and CDT-based SFT may provide a useful analytic tool. Work in this direction 
is in progress. 

Equally interesting is the possibility of performing a summation over wordsheets of 
all genera. Again, since the double-scaling limit in CDT-based SFT is different from the 
double-scaling limit in non-critical string theory and the penalty for creation of a higher 
genus surface is larger in the sense mentioned above, viewing the creation a higher 
genus worldsheet as a successive creation and annihilation of a baby universe, one 
could hope the result of such a summation is better behaved and less ambiguous than 
was the case in non-critical string theory. Work in this direction is also in progress. 
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CHAPTER 



Matrix models for two-dimensional Euclidean quantum 
gravity 



In the previous chapter we computed higher-loop and higher-genus amplitudes in the 
framework of a SFT for the continuum model of generalized CDT. An earlier formulation 
for DT was introduced by Kawai et al. [90J as discussed above. If one would like to 
perform similar calculations in the discrete framework of DT one could either formulate 
a discrete STF as has been done in [94] or make use of so-called matrix models and 
loop equations derived from them. 

In this chapter we want to introduce the use of matrix models to solve combinatorial 
counting problems of DT on the discretized level. Those techniques are very powerful, 
since they allow for many generalizations especially in the context of matter couplings. 

11 .1 I Counting triangulations with matrix models 

In the following we want to give a short introduction to how certain matrix integrals can 
be used as generating functionals for random triangulation. The basic idea underlying 
this construction goes back to the seminal work by 't Hooft HZU on the large-N limit 
of QCD, followed by work on the saddle point approximation JgH 1180] . For further 
references the reader is referred to several reviews | [74ll181[[73l[T8] . 

Before entering the discussion of matrix integrals let us first briefly recall the Feyn- 
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Figure 11.1: Building blocks for Feynman graphs: (a) the scalar propagator and (b) the 
scalar three-point vertex. 



man expansion of zero-dimensional ^p^ field theory whose action we already mentioned 
in ([9]6), 



with 



ln{g) 



'Itx 



oo 

E 

n=0 



— ; ( 



S 3 

3!^ 



(11.1) 



(11.2) 



The integral ( |1 1.1| | is formal since it is clearly not convergent until some analytic contin- 
uation is performed. In terms of Wick contractions the integral ( |1 1 .2) counts all possible 



ways of connecting n three-valent vertices (Fig. |1 1 .1 (b)) by pairings, represented by 



the propagator {(p(p) (Fig.|1 1 .1 (a)). These Feynman diagrams already look very much 



alike the dual triangulations in DT. However, to represent Riemannian geometry more 
structure is needed. This structure is provided when replacing the usual Feynman 



graphs by so-called "fat-graphs" as shown in Fig. |11.2| ln terms of the integral expres- 
sion ( |1 1 .2| ) this means that we replace the field ip by traces over N x N Hermitian 
matrices 4> = {(pafs), i-e. 



/oo 
c/0e-'^(^^^<^^-S*^^^) = 5^/„(^, N) 
n=0 



with 



ln{g, N) = j d4>e 



n\\ 3 ^ 



(11.3) 



(11.4) 
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(a) 



Figure 1 1 .2: Building blocks for fat-graphs: (a) the Hermitian matrix propagator and (b) 
the Hermitian matrix three-point vertex. 



Here the integration measure is defined using the matrix elements 

a<f3 a<l3 



(11.5) 



The propagators become now oriented links (see Fig. 1 1 .2 (a)) 



(11.6) 



and the three-valent vertices become also oriented (Fig. |1 1 .2| ;b)). From this we see that 
the integral ( |1 1 .4| accounts for all closed, possibly disconnected, unrestricted triangu- 
lations with n triangles (recall Fig. [53J. Since each vertex comes with a factor of gN 
and each propagator with a factor l/N one can convince oneself that the contribution 
of a given triangulation with n triangles carries a factor 



(11.7) 



where the Euler characteristic of the triangulation is given by x{T) = 2 - 2h, with h 
being the genus of the triangulation]^ Further, it is noted that the factor 1/n! in | |1 1 A) 
combines with the number of permutation of the vertices yielding the factor 1/C{T) 
which already appeared in ( |5.5| ). To make contact with dynamical triangulations we 



^The relation | |1 1 .7} is also valid for the case of triangulations with boundaries, where the Euler 
characteristic becomes x( 7") = 2 - 2/? - fa with fa being the number of boundary components. 
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see that we have to assign g with the bare dimensionless cosmological constant and 
N with the bare dimensionless Newton's constant Gn as follows|3 

g = e-\ A/ = e^/^". (11.8) 

Eq. ( |1 1 .71 ) now corresponds to the Bolzmann weight of the triangulation. We conclude 
that the partition function for dynamical triangulations is given by 

-(-")^^' <-^' 

where the logarithm is introduced to precisely account for the number of connected 
(non-restricted) triangulations. Here Z{g, N) as introduced in ( |1 1 .3[ ) is commonly writ- 
ten as 

Z{g.N) = ld<Pe-''^^^^\ = - f 0^ (11.10) 

where V(0) is called the potential of the matrix model. From ( |1 1 .7) we see that Z{g, N) 
has the expected topological expansion 

oo 

Z{g,N) = Y,N'-''Z,{g). (11.11) 

In the so-called planar or large-A/ limit one takes N ^ oo and obtains the genus zero 
partition function 

Zo{g) = N-^Z{g, N) + 0(l/A/2). (1 1 .12) 

We have seen how to represent the partition function of DT as an integral over 
Hermitian matrices. One of the strengths of using matrix models is that it allows for 
many generalizations. One possible generalization is to replace the potential \/(0) in 
fTTlO) by 

oo 

This integral counts for closed, possibly disconnected, general non-restricted triangu- 
lations which include j-gons, each weighted by the corresponding coupling gj. Gener- 
alizing dl 1 .7| to arbitrary potential gives 

oo 

C-i(7)A/^(^) (11.14) 



^The relation for g is tine same as used in previous chapters, only that we sometimes written it in 
terms of the bare dimensionfull cosmological constant, i.e g = e^^^^. 
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where N{J) is the number of j-gons in the triangulation. 

Having established the connection between matrix models and dynamical triangu- 
lations we now want to show how to actually calculate certain quantities, such as the 
disc function, within this framework. 



11 .2 1 Loop equations from matrix models 

Let us for the following consider the matrix model with general potential as described 
in ( |1 1 .13[ ). We abbreviate the set of couplings by g = {gi, g2, gs, ■■■). 
Denoting the expectation value by 

we see that 

wi{g,N) = ^{tr<p') (11.16) 

accounts for the disc function for general non-restriucted triangulations with boundary 
length / and arbitrary genusj^ We can now also introduce the generating function 

w{g, ^■^) = f2 ^'(^' = ^ ('^ ^) ^ ^) 

It is sometimes useful to further express the disc function in terms of the density of 
eigenvalues p(A) defined as 

P(A) = (^XJ5(A-A,)^, (11.18) 

where A,-, / = 1, 2, .., A/ are the eigenvalues of the matrix 0. We can now write the disc 
function as ^ oo \ 

wi{g,N)= I dXp{X)X', w{g,N,z)= I dX-^^. (11.19) 

In the limit A/ ^ oo the density of eigenvalues has finite support [c , c+] on the real axis. 
This implies that the analytic continuation of w{g, N, z) is well defined on C/[c_, c+] and 
we can invert ( |1 1.19[ ) as follows 

27r/p(A) = lim {w{g, A/, A - ie) - w{g, N,X + ie)) . (1 1 .20) 



^The factor of 1/N is introduced to cancel tine factor A/^ ^, b=l coming from tine Eulerclnaracteristic. 
TInis is done to make tine leading order in the large- A/ limit of order one. 
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To derive the loop equations for the disc function from the matrix model [182111721 
I173| we will make an infinitesimal change in the field variable of the functional integral. 
This is a standard method employed in field theories. Let us consider the following 
infinitesimal field redefinition 

5:0^0 + 6^-, (11.21) 
z-cj) 

where e > is an infinitesimal parameter. In the limit N ^ oo this field redefinition 
is well defined for z e C/[c_,c+]. Under this field redefinition the measure and the 
potential transform to first order in e as 

c/0 ( 1 + ) (11.22) 

\ z-cj) z-cj) J 

and 

tr \/(0) ^ tr V{(i)) + e tr ^'(0)) ■ C 1 -23) 

Using that log Z(g, N) as given in ( |1 1.13[ ) should be invariant under the transformation 
8 we obtain the following identity 

5\ogZ{g,N) = Z-\g,N) j d<p\^[t^^^ -A/tr (^-1^ \/'(0)^ | e"'^*^ = 0. 

(1 1 .24) 



The first term in the brackets is by definition (see ([lT40)) 

1 




N'w\g,N,z) + w{g,N,z,z), (11.25) 

where w{g^, N, zi, zs) is the two-loop amplitude, i.e. the sum over all triangulations with 
two-boundaries with boundary cosmological constants zi and Z2. The second term can 
be evaluated using the density of eigenvalues, 



N \ z — (f) / z — A Jd; 2ni z — w 

where the first equality follows from ( |1 1.19[ ) and the second equality from ( |1 1 .20} . The 



contour C encloses counterclockwise the cut [c , c+] of w{g, N,z) but not z as shown 
in Fig. |1 1 .3[ As discussed above ([TTTT9]| for this relation to hold we must have N oo 



which implies finiteness of [c , c+]. Inserting | |1 1 .25| ) and ( |1 1 .26| into ( |1 1 .24) yields the 
loop equation 

dw V'(w) . 1 

^^->w{g, N, z) = w'{g, N, z) + —w{g, N, z, z). (1 1 .27) 

Zni z — w — — /V^ — 
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Figure 1 1 .3: The contour C encloses counterclockwise the cut [c_, c+] but not z. 



Since we are only interested in the genus zero disc function we can take the large- A/ 



limit which allows us to drop the last term from ( |1 1 .271 , yielding 

dw V'{w) 



Q 7.11] z 



w 



-w[g,z) = w {g,z). 



(1 1 .28) 



The solution to | |1 1.28[ ) reads 

w{g,z) = \ {V'{z) - M{g,z)^{z-c^{g)){z-c_{g))) , (11.29) 

where M(z) is a polynomial whose degree is one less than V"(z). Further the density 
of eigenvalues reads 

1 



P(A) 



27r 



M{g,\)J[c+[g)-\)[X-c^[g)). 



(11.30) 



As for the loop equations of DT, the solution for the disc function ( |1 1 .29| is completely 
determined by the requirement that w{g, z) = 1/z + 0{l/z^). In particular, one can de- 
rive the following equations by imposing this requirement for ( |1 1 .29) and using ( |1 1 .28| 
(see for example [18] for more details): 

dw V'{w) 1 
'Co. 27r/ 1/1/ - z ^(m/- c+)(m/- c_) 



M{g. z)= (f 



(11.31) 
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and 



where we defined 



M_i(|;,c_,c+) = 2, Mo(|;,c_,c+) = 0, 

V'{w) 



M,{g,c_,c+)= (f ~ 



(11.32) 



(11.33) 



_)/c+l/2(^_ C_)l/2- 

Here the first equation ( |1 1.31| determines the polynomial M(g, z) as a function of c , c+ 
and the couplings g, while the second equation ( |1 1.32[ ) yields two algebraic equations 
which determine c (g) and c+(g) as functions of the couplings g. 
Let us now \ooW at two specific cases of the potential. For 

1 ,o 



V{ct>) 



(1 1 .34) 



we expect to obtain the generating function for rooted branched polymers. Inserting 
V'{z) into dl 1 .29| and using the conditions ([T]T31}-([lT32]| for the requirement that 

w{z) = 1/z + C(l/z2) we obtain 



(11.35) 



^z) = -(^z- v/(z-2)(z + 2) 
which is precisely ( |5.14| . Further, using ( |1 1 .30) one gets for the density of eigenvalues 

p{X) = ^^{z-2){z + 2) (11.36) 

which is Wigner's semicircle law for the eigenvalue distribution of Hermitian matrices 
| T83irT84| . 

For the potential 



v{cf>) = -4>' 



(11.37) 



2^ 3^ 

we expect, according to the considerations made above, to obtain DT. Inserting ( |1 1 .37} 
into dl 1.291 we get using | |1 1 .31h ( flT32) 



w{g, z) = ^[z- gz^ - {gz - c{g))^{z - c4g)){z - c_{g)) 



(11.38) 



which coincides with ( |5.2.2| |, where c{g), c (g) and c+(g) are determined by the same 
algebraic equations as in ( |5.2.2| |. 

We have seen how matrix models can be used to count generalized triangulations. 
Further, we showed how we can obtain DT from matrix models for the specific potential 
V(0) = 0^/2 - g(j)^/3. In particular, we computed the disc function. In the following 
section we want to extend our results by calculating higher-loop amplitudes. 
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11 .3 I The loop insertion operator and higher loop amplitudes 

Even though in Chap. |5] we only calculated the one-loop amplitude for DT (disc func- 
tion), we already encountered higher loop amplitudes in this thesis, namely the two- 
loop amplitude for DT in the continuum, i.e. ( |9.27| . 

IVIatrix model higher-loop amplitudes, i.e. amplitudes for generalized triangulations 
with b boundaries of lengths 4, 12, //,, can be defined as follows 

dg) = N'-mtr4>'^)-.-{tr4>'^))^, (11.39) 

where (...)^ denotes the connected expectation value defined as 

(Oi((/.) • • • O„(0)), = (Oi(0) ■ ■ ■ O„(0)) - (Oi(0)) ■ ■ ■ (O„(0)) . (1 1 .40) 

Equivalently one can also use the generating functional for the higher-loop amplitudes 

w{g,z^ z,) = N'-'({tr^—)...{tr^—)) . (11.41) 

- \Zi-0 Zb- (p / ^ 

To calculate the higher loop amplitudes we introduce the so-called loop insertion 
operator | fT72l [T85| 



dV{z) j^^zJ+^ dgj 

For a given potential with fixed coupling constants gf one uses these relations in the 
following way. Assume that gj can vary, one now acts with the loop insertion operator 
6/dV{z2) say on the disc function w{g, zi), and then set gj = gf again. Each term in the 
summation over j has the effect that it removes a >sided polygon^ leaving a marked 
boundary of length j to which the new boundary cosmological constant is assigned. 

We now obtain the higher-loop amplitudes by successive action of the loop insertion 
operator on the disc function 



"^We showed above that each triangulation comes with a weight proportional to H/^i g'"^''. where N{i) 
is the number of /-gons in the triangulation. Acting with 6/6gj therefore reduces the number of j-gons 
by one. Further it brings down a factor of N(J) which precisely accounts for the number of possibilities 
to choose a specific j-gon. 
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Inserting the solution for the disc function ( |1 1 .291 into ( |1 1 .431 one obtains for 6 = 2 
the two-loop amplitude for general potential || 182[l172[pr73J . 



w{g,Zi,Z2) 



Z1Z2 - |(c_ + C+)(zi + Z2) + C_C_| 



2 (Zi - Z2)2 I ^(Zi - C_)(Zi - C+)v/(z2 - C_)(Z2 - C+) 



(1 1 .44) 



The dependence on the potential is implicit through c (g) and c+{g). For the case of 
DT this expression can be used to derive the continuum two-loop amplitude ( |9.27| . 

We will make further use of the loop insertion operator in Chap. [Tsl where we will 
look at higher genus contributions. It is clear from the loop equation ( |1 1 .27| and the 
genus expansion 

00 

w{g. N.z) = J2 N'^'MS, z) (1 1 .45) 



that to calculate the higher genus disc function higher-loop amplitudes are needed. 
This is similar to the situation when calculating higher genus amplitudes from the SFT 



in Chap.QO} We will see in Chap.QSjthat this is not a coincidence, since the CDT-based 
SFT can be described by a matrix model. 



11 .4 1 Discussion and outlook 

In this chapter we introduced matrix models to facilitate combinatorial counting prob- 
lems in DT. Matrix model techniques are very powerful to compute many quantities on 
the discrete level and to study further generalizations of those. 

We explicitly saw how the disc function and higher-loop functions can be computed 
within this context. This construction can also be extended to higher genus contribu- 
tions of the disc function and higher-loop amplitudes as we also commented on above 
(see, for instance, [185] or the more recent papers |186[|187| ). 

Further, we saw that matrix model techniques generalize easily to arbitrary poten- 
tials. This interesting fact allows one to study simple matter models, so-called minimal 
models, coupled to the quantum geometry (see [188] for the simplest example). In 
the continuum limit of the pure DT matrix model the first moment of the potential was 
tuned to zero while the second moment was non-zero. This situation corresponds to 
the so-called 2nd multi-critical hypersurface in the space of coupling constants. Being 
able to introduce higher terms in the potential we can generalize this situation by going 
to the mth multi-critical hypersurface, where the first k moments, k < m oi the potential 
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are zero. It was seen that those minimal models correspond to so-called (2,2m - 1)- 
conformal field theories coupled to two-dimensional Euclidean quantum gravity. 

The so-called two-matrix model allows for further applications with respect to mat- 
ter coupling [TM HM [T9T1 [TM [TM [TM [Tig. This model is used to study the Ising 
model coupled to two-dimensional Euclidean quantum gravity, where the two matrices 
are used to describe the triangles with up-pointing and down-pointing spins respec- 
tively. We will further comment on this in the final discussion of this thesis. 
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CHAPTER 



A loop equation and matrix model for CDT 



Matrix model techniques have been a useful tool in the field of DT. However, until re- 
cently 11331 [T96l [T971 no matrix model for two-dimensional CDT has been found. In this 
chapter we introduce a newly discovered matrix model for discretized CDT |196[|197| . 
After motivating geometrical loop equations for the "bare" CDT model in the next sec- 
tion, we see that to make the WnW to matrix models it is important to include (weighted) 
spatial topology changes. In particular, we present a matrix model derivation of the 
generalized CDT disc function as derived in Chap. [9] purely from continuum methods. 



12.1 I Geometrical loop equations for CDT 



In this section we compute the generating function for a set m//,a/ of triangulations which 
are somewhat similar to the original set of causal triangulations and which lead to the 
same continuum physics: Let N denote the number of triangles and / the number of 
linl^s at the boundary (which has one marl^ed WnW), and assume the topology is that of 
the 6\sW and denote the generating function w{g, x): 

oo oo 

w{g, x)=Y, ^i.Ng'^^' = Yl ^'(^)^'- 2.1 ) 

/,A/=0 /=0 

Here we adopted the definition of the generating function as used in ( |5.9| , where g is 



the fugacity of a triangle and x is the fugacity of a boundary edge. The triangulations 
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-O 



(a) 



Figure 12.1 : Illustration of the two composition moves to add a triangle. The white dot 
on the right-hand-side shows the position of the mark before the triangle was added 
whereas the blacl< dot shows the mark after the triangle was added. 



considered here do not have the strict time-sliced structure as those considered in 
Chap.|6](e.g. see Fig.|6.1|), but have a spiral structure, since they are derived using loop 
equations, i.e. through a peeling procedure (see Fig. [731 . Nevertheless, we will see 
that this construction leads to the same continuum expressions as obtained in Chap. 
[6} This is similar to the derivation of the fixed geodesic distance two-loop function 



for DT considered in Sec. |5.3[ where both the transfer matrix method and the peeling 
procedure lead to the same continuum results. 

The triangulations can be generated by recursively adding triangles. In our model 
there are two possible moves. Firstly, one can glue two edges of the additional triangle 
to the triangulation, one to the marl^ed edge and the other one next to it in the clocl^wise 



direction (Fig. |12.1| (a)). Secondly, one can add a triangle by simply gluing one of its 
edges to the marl^ed edge of the triangulation and assigning the new mark to the 



new edge further clocl^wise (Fig. 12.1 (b)). Together, the two moves give the following 



generating equation for large n and / (see Fig. |12.2| |, 



w{g,x) 



g 



{w{g, x) - wi{g)x - 1) + gxw{g, x). 



(12.2) 



The one in front and the subtraction of wi{g)x + 1 are related to the initial conditions. 
Equation ( |12.2| | is a simple linear equation and the solution is given by 



w{g,x) = g 



(1/g - l/1/i(g)) X 

gx^ -x + g 



(12.3) 



where as in previous situations wi{g) can be determined by demanding that the sin- 
gularity structure of (|12.3[) does not change discontinuously near g = 0. The poles of 
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Figure 12.2: Graphical representation of the loop equation ( |12.2| | 



( [123) are located at 



1 ± y/l - Ag^ 



(12.4) 



Since the expansion of wi{g) needs to be a power series, we have that ivi =x_. Hence 
the disc function is given by the following simple expression 

1 



w{g,x) 



Wi{g) = Wi{g) 



1 - wi{g)x' 

Using the same scaling relations as in the transfer matrix formalism of CDT, 

1 



g 



X 



-aX 



we reproduce the continuum disc amplitude of CDT, e.g. ( |6.51 

1 



X 



A 



1/1/a(/-) 



,-VA/. 



(12.5) 



(12.6) 



(12.7) 



12.2 1 A matrix model for generalized 2D causal quantum gravity 

To include weighted spatial topology changes, as has been done in Chap. |9] in the 
continuum, we reintroduce the quadratic term in the loop equation | |12.2| 

wp{g, x) = 1 + - {wf3{g, x) - wi{g, (3)x - 1) + gxwp{g, x) + (3x^wp{g, xf, (1 2.8) 

where /5 is a coupling constant that determines the rate of the spatial topology fluctua- 
tions (see Fig. |12.3[ ). 

To conform with matrix model conventions it is useful to introduce the following 
notation 

Wp[g,X= 1/z) 



wp[g, z) 



(12.9) 
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Figure 12.3: Graphical representation of the loop equation (12.8). 



which has already been employed in ( |5.1 1| 

With these conventions the loop equation is given by 



where 



f3w0{g, zf - V\z)wp{g, z) + Qp{g, z) = 0, 



V[z) = -gz + \z^-\gz\ V'[z) = -g + z-gz' 



and 



Qp{g,z) = l-g{w^{g,(3) + z). 
The solution of the loop equation | |12.10[ ) is of the following form 



(12.10) 

(12.11) 
(12.12) 



MS,^) = ^ (v'{z) - ^V'{zy-APQ,{g,z)^ . (12.13) 

At this stage the solution of the disc function is still implicit since it depends on wi{g, (3) 
through | |12.12[ ). Demanding the solution to have only one cut in the complex z plane, 



as discussed in Sec. |1 1 .2[ gives the explicit solution 

Ms, ^) = ^ (-^ + ^ - gz^ + {gz - c)^J{z-c^){z-c^)^ 
where c is the solution of a third order polynomial, 

2c^ - 3c2 + (2g' + 1) c = g\l - 2/5), 

and 



(12.14) 



l-c±V2^{l-c)c-g^ 



g 



(12.15) 



(12.16) 
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Written in the form of ( |12.10[ ) it is easily seen from the considerations in Sec. |1 1 .2 
that this equation corresponds to the loop equation of a simple matrix model, 

Zf,{g,N) = J c/0exp(^-|tr\/((/))y (12.17) 

where is a A/ x A/ Hermitian matrix and the functional form of the potential V(0) is 
given by 

\/(0) = + 1^2 _ 1^^3 (12.18) 



We see that in comparison with the matrix model of DT, i.e. ( |1 1 .37| |, we have an 



additional linear term in the potential. For later convenience we also took the factor 1//3 
out of the potential, but alternatively one could have included it in the potential to be in 
accordance with the conventions used in | |1 1 .13| |. 

12.3 I Taking the continuum limit 



In the continuum limit of DT, as described in Sec. |5.2.2[ or equivalently of the matrix 
model for DT, the critical value for the boundary cosmological constant z coincides with 
the critical value of c+ only, i.e. = c+{gc)- As a result, this standard limit has the 
peculiar feature that it leaves a non-scaling term as a memory of the discrete theory, 
since 

^) = ^ns + 3'/'WlC''\Z) + 0{3^). (12.19) 

where Wns = {z - gz^)/2 is non-scaling part and l/l/^^"^(X) the continuum disc function 
of DT 

Wl-\Z) = - ^) \fzWK (1 2.20) 

In our specific model the matrix potential is such that the critical points of c+{gc) 
and c (gc) coincide. This leads us naturally to the universality class of two-dimensional 
CDT implying the same scaling relations as before ( |12.6[ ), provided one also scales the 
coupling constant (3 in exactly the same way as described in ( |9.7| |, 

(^-\gsa\ c = \e^^. (12.21) 

Contrary to the standard continuum limit of the one matrix model as employed by DT 
our new continuum limit is free from leading non-scaling contributions. Inserting the 
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scaling relations of the bare model (|12.6[), i.e. 



g=\e-'"^'\ z = e 



aZ 



(12.22) 



and ( |12.21^ into | |12.14D yields 

3 

where the continuum disc function W^gXZ) is given by 

= ^ {-{Z' - A) + (Z - C)V{Z+Cr-2gs/C 



(12.23) 



(12.24) 



Inserting the scaling relations ( |12.22| and ( |12.21| | into the algebraic equation for c 
( |12.15| ) further gives 

C^-AC + gs = 0. (12.25) 

We see that ( |13.1| ) and ( |13.1[ ) agree precisely with the solution for the generalized CDT 
model as derived in Chap. [9] purely in a continuum frameworl^, i.e. ( |9.19| )-( |9^20] |. 



12.4 1 The finite time propagator 



Having derived the disc function of generalized CDT from a matrix model we can now 
also implement the decomposition moves of the above loop equation into an corre- 
sponding peeling equation for the finite time propagator. This is in analogy to the 
calculation of the fixed geodesic distance two-loop function in DT which has been de- 
scribed in Sec. 15.3.21 



Recall from Sec. |5.3.2| that in the peeling procedure it is natural to interpret the loop 
equation as a time dependent process [|g?ll125[[T98| . where each decomposition move 
on a boundary of length / is seen as a "V/-th" part of a time step. Following 1197] we 
can obtain an equation for the finite time propagator similar to \5.39) and ( |5.40| for DT. 
For large /i we have 

d 

^G^(/i, l2;g; t) = gkGp{k - 1, l2\g\ t) - liGp{k, h-.g; t) + gkGp{k + 1, l2\g\ t) + 

oo 

+ 213k J2 ''lis. P)Gp{li - / - 2, h: g: t). (1 2.26) 

/=o 
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Here the first term on the right-hand-side is new in comparison with ( |5.39| and it pre- 
cisely corresponds to the new move in the loop equation. In terms of the generating 
function we have 



d d 
-Gp[z,w;g;t) = — [{-g + z - gz^ - 2pWf}{g, z))Gp{z, w; g; t)] 



(12.27) 



Using the scaling relations defined above, i.e. ( |12.22| | and ( |12.21| |, and letting time 
scale as 

T^ta (12.28) 



we obtain 

d 



d_ 
dX 



(X2-A + 2g,l/l/A,,^(X)) Gj,,JX,Y- T) 



(12.29) 



This is precisely the result obtained in Chap.|9} However, whereas in Chap.[9]we were 
generalizing CDT within a continuum frameworl<, we now found an underlying discrete 
theory reproducing those results. 



12.5 I Discussion and outlook 

In this chapter we developed a loop equation and a matrix model for generalized CDT 
|196[ [T97J. For the continuum limit of this matrix model it is important that both end- 
points of the cut are fallen to the same critical point, i.e. c+{gc) = c_{gc) = z^. To 
accomplish this it is essential that the matrix model has a linear term in the potential 
with the same coupling as the cubic term. We gave an interpretation of this linear term 
in terms of a new move in the loop equation. 

In the next chapter we will see an intriguing aspect of this new matrix model con- 
tinuum limit, namely, that it does not require N to scale with the cutoff in the so-called 
double scaling limit | |199[ [2001 120 1| . IVIoreover the continuum results are also described 
by a matrix model. This matrix model has of course both an expansion in the coupling 
constants gs and A, as well as a large- A/ expansion in powers of l/A/^, which reorga- 
nizes the power expansions in gs and A in convergent "subseries". By comparing with 
the generalized causal dynamical triangulation model the powers of A/^''+2 in the large 
N expansion can be identified with the continuum causal dynamical surfaces of genus 
h [133| . In this sense our new continuum limit leads to a picture that is much closer in 
spirit to the original idea by 't Hooft [1 79] for QCD. 
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An advantage of having a matrix model description is its computational strength. 
Until now, only the transfer matrix formalism has been available for analytical computa- 
tions in CDT. Here we extended the available tools and presented new, more powerful, 
matrix model and loop equation methods. Importantly, these new methods allow one 
to analytically study simple matter models coupled to two-dimensional causal quan- 
tum gravity, such as minimal models or the Ising model. Several of these models are 
currently under investigation. 

It is interesting that two-dimensional CDT without spatial topology change is not 
described by a matrix model. Pure CDT is only recovered when taking the limit gs ^ 
of the loop equations for generalized CDT. There exists a rather complicated matrix 
model description for the one-step propagator of three-dimensional CDT [137[|139| . It 
seems quite plausible that if one allows for branching baby universes weighted by a 
coupling gs the matrix model could greatly simplify. 



CHAPTER 



A continuum matrix model for CDT 



In the previous chapter we found a matrix model for two-dimensional causal quantum 
gravity defined through CDT. In this chapter we show that not only the discrete, but also 
the continuum model of CDT can be described by a matrix model [133]. This is a very 
interesting result, since in the DT context matrix models have only been used on the 
discretized level and not in the continuum. 



13.1 I Mapping the DSE 

One interesting observation is that the continuum expression of the disc function for 
generalized CDT ( |9.19| , as derived in Chap.[9| is very similar to the expression of the 
disc function obtained from matrix models, i.e. | |1 1 .29[ ). Writing ( |9.19| ) in the following 
way it becomes more obvious: 

= ^ (-(x^ -^) + {x- c)v/(x-c_)(x-c+)) . (13.1) 

where the constants C, C± are determined by 

C^-AC + gs = 0. C± = -C±v^i^. (13.2) 
This analogy can be pushed much further. Introducing the notation 

V'{X) = -{A-X^), V{X) = - Ux -^xA, (13.3) 
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we can rewrite the DSE from the CDT-based SFT, i.e. ( |10.37| )-( [Ta!39l l, in the following 
way 

= dx( -V'{X)W{X)+W^{X) + aW{X,X)) - —, (13.4) 

= dx([-V'{X) + 2W{X)]w{X,Y) + aW{X,X,Y)^ + 

dY([-V'{Y) + 2W{Y)]W{X, Y) + aW{X, Y, Y)j + (13.5) 

= dx{[-V'{X) + 2W{X)]W{X,Y,Z) + aW{X,X,Y,Z)^ + 
dY{[-V'{Y) + 2W{Y)]W{X, Y,Z) + aW{X, Y, V.Z)) + 

dz([-V'{Z) + 2W{Z)]W{X, Y,Z) + aW{X, Y,Z,Z)j + (13.6) 

2dx[W{X, Y)W{X, Z)] + 2dY[W{X, Y)W{Y, Z)] + 

W{X,Z)-W{Y,Z) 



2dz[W{X,Z)W{Y,Z)] + 2( dxd 



X-Y 



W{X, Y)-W{Y,Z) W{X,Y)-W{X,Z) 
OxOz X - Z Y - Z 

Let us introduce the expansion^ 



l/l/(Xi X„) = 5^a''l/l/,(Xi (13.7) 



h=0 



As shown in Chap.QOl h can be interpreted as the number of handles of the worldsheet, 
and the equations above can be solved iteratively in h. More precisely, the equations 
at order a° allow us to determine Wo{x), Wo{x,y), and similarly the equations at 
general order determine l/l//,(x), l/14(x,y), etc. For example, one finds 

Wo{X) = 1(V'{X) + -{X - C)v/(X-C_)(X-C+)), (13.8) 
2\ gs / 

1 1 / xy- |(C_ + C+)(X + V) + c_c+ \ 
Wo(X. Y) = -^^^ I _ _ ^^^^^^ _ - II , (13.9) 



^ Note that both W and Wt, are still g-s-dependent, although we do not write the dependence explicitly 
here. 
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where the constants C and C± are defined as in ( |13.2[ ) above. 

Writing the amplitudes in this fashion leads one to the surprising realization that 
Wo{X) and Wo{X, Y) indeed coincide with the large-A/ limit of the disc function and 
two-loop function of a Hermitian matrix model 

Zgs{^ = j d<^e-^''^^'^^ (13.10) 

with potential 

\/(cD) = AcD _ J^cD3 (13.11) 
gs 3gs 

This is a potentially exciting result, because so far no standard formulation in terms of 
matrix models has been found for a continuum model of quantum gravity. We will in 
the following section show a more general result which identifies the Dyson-Schwinger 
equations derived above with the loop equations of a Hermitian matrix model with the 
cubic potential ( |13.1 1| . 



13.2 I Matrix loop equations 



In Sec. |1 1 .3| we explained how in a matrix model context loop equations for higher-loop 
amplitudes can be derived from the loop equation of the disc function by applying the 
loop insertion operator sufficiently many times, as summarized in ( |1 1 .43| . In order to 
compare with the Dyson-Schwinger equations of our string field theory, we differentiate 
the equations obtained with respect to X. For the potential ( |13.1 1| ) we arrive at the 
following equations 

= dx(^-V'{X)W{X)+W^{X) +^W{X,X)^ -J, (13.12) 
= d,([-V'{X) + 2W{X)]W{X,Y) + ^W{X,X,Y)^ + 

^a.a.(m_Mr)), ,3.13) 

= dx{[-V'{X) + 2W{X)]W{X,Y,Z) + ^W{X,X,Y,Z)j + 

2dx(w{X,Z)W{X,Y)^ + (13.14) 



X-Y J " ^\ X-Z 
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Using that W{Xi Xn) is a symmetric function of its arguments, we see that eqs. 

( |13.12| -( pr3.14| lead to exactly the same coupled equations for 1/1/ as do ( |13.4| |-( pi"3^ 
//we identify 

a = i,. (13.15) 



In this case the discussion surrounding the expansion ^31^ is nothing but the standard 
discussion of the large- A/ expansion 

oo 

l/l/(Xi Xn) = ^N-^^ yVh{Xi Xn) (13.16) 

of the multi-loop correlators (see [185J or for more recent developments [186[ I187| ). 
The iterative solution of these loop equations is uniquely determined by Wa{X) provided 

that W{Xi Xn) is analytic in those X, that do not belong to the cut of the matrix 

model. 



13.3 1 Relation to the discrete matrix model 

In the previous chapter we saw how the discretized model of CDT could be described 
by a matrix model with the following partition function and action 

Z^,'\g, N) = j c/^e-^'^'I^l S(^)[0] = |tr (-^0+^0^ - 3^0^) ■ (13.17) 

In comparison, in this chapter, we obtained a matrix model for the continuum model of 
CDT, described by 

Zj;)(A, N) = j dcDe-^*''!*], S(")[cD] = ^tr (^AcD - ^$3^ . (13.18) 

It is a very surprising result that both matrix models can be related by taking the con- 
tinuum limit on the level of the action. Inserting the standard scaling relations 

(^=\gs3\ g = \e-'"^'^ (13.19) 

together with a scaling of the matrices 

00 ^ 

= e^* = 5^-(act))", cDO = i„,„ (13.20) 
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into the action ( |13.17| | of the matrix model for the discrete CDT model gives 

S('')[(/.] = ^tr (^AcD - ^cD^^ + const. + C(a). (13.21) 

Hence we see that the matrix model for continuum CDT follows from the one for the 
discrete model by simply taking the continuum limit on the level of the action. 

This situation is very different to the one encountered in the matrix model for DT. 
Due to the absence of the coupling f3 one has to perform a double scaling limit in which 
N scales as 

/V(a) = e^/^"(^' = (^322) 
such that the dimensionless Newton's constant 

^ ^ ^ log(Aa') = log(A/(a)a5/2) (13.23) 



G/v(a) 4 

remains finite as a ^ and N ^ oo. In contrary, in our continuum CDT matrix model 
N does not have to scale with the cut-off. 



13.4 1 Discussion and Outlook 

Let us consider the matrix model corresponding to the potential ( |13.1 1| ). We can per- 
form a simple change of variables (t> ^ -(t> - y/K \n the matrix integral to obtain a 
standard matrix integral 

Zg^{m, N) = j d<^ Q-mrV(<^)^ (13.24) 
where the new potential (up to an irrelevant constant term) is given by 

\/(cD) = 1 (^^mct)2 + 1(1)3^ , m = 2Va. (13.25) 

It is amusing to note that the matrix integral ( |13.24| | is precisely the kind of matrix 
integral considered in the so-called Dijkgraaf-Vafa correspondence [I202[ I203| . where 
in this case V{<i>) is the tree-level superpotential of the adjoint chiral field (t>, which 
breaks the supersymmetry of the unitary gauge theory from A/" = 2 to A/" = 1. If one 
demands that this tree-level potential correspond to a renormalizable theory, its form 
is essentially unique, and precisely of the form ( |13.25| ) originally used by Dijkgraaf 
and Vafa, with gs a dimension-three coupling constant coming from topological string 
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theory and in the DV-correspondence related to the glueball superfield condensate in 
the gauge theory. 



In the "old" matrix model representation of non-critical strings and two-dimensional 
gravity one had to perform a fine-tuning of the coupling constants in order to obtain a 
continuum string or quantum gravity theory. This implemented the gluing of triangles 
(or, more generally, of squares, pentagons, etc.) which served as a regularization of 
the worldsheet. The fine-tuning of the coupling constants reflected the fact that the link 
length of the triangles (the lattice spacing of the dynamical lattice) was fallen to zero in 
the continuum limit. The situation here is different. Although CDT can be constructively 
defined as the continuum limit of a dynamical lattice, we have in the present \NorW 
been dealing only with the associated continuum theory. Thus in our case the matrix 
model with the potential ( |13.1 1| (or ( |13.24| ) already describes a continuum theory of 
two-dimensional quantum gravity. Its coupling constants can be viewed as continuum 
coupling constants and the role of N is exactly as in the original context of QCD, namely, 
to reorganize the expansion in the coupling constant gs. 't Hooft's large-A/ expansion 
of QCD is a reorganization of the perturbative series in the Yang-IVIills coupling gvM, 
with 1/N talking the role of a new expansion parameter. In this frameworl^, after the 
coefficient of the term l/N^^ of some observable has been calculated as function of the 
't Hooft coupling gf^ = g^^N, one must tal^e A/ = 3 for SL/(3), say. The situation in CDT 
string field theory is entirely analogous: starting from a perturbative expansion in the 
"string coupling constant" gs (in fact, in the dimensionless coupling constant g^s/A^/^, as 
described in previous chapters), we can reorganize it as a topological expansion in the 
genus of the worldsheet by introducing the expansion parameter a. For the multi-loop 
correlators this expansion is exactly the large-A/ expansion of the matrix model ( |13.1 1| ) 

and the coefficients, the functions l/l//,(Xi Xn), are exactly the multi-loop correlators 

for genus-/? worldsheets of the CDT string field theory with a = 1. 



As a "bonus" for our treatment of generalized (and therefore slightly causality- 
violating) geometries, we also obtain a matrix formulation of the original two-dimensional 
CDT model as discussed in Chap, [ej where the spatial universe was not allowed to 
split. Worl^ing out the limit as ^ of the various expressions derived above, we see 
that this model corresponds to the large-A/ limit of the matrix model where the coupling 
constants go to infinity, but at the same time the cut shrinl^s to a point in such a way 



13.4 Discussion and Outlook 



177 



that the resolvent (or disk amplitude) survives, that is, 

Wo(X) ^ = W!:''"\X). (13.26) 

The existence of a matrix model describing the algebraic structure of the DSE leads 
automatically to the existence of VIrasoro-like operators /.„, n > -1 [182[ l172| . which 
can be related to redefinitions of the time variable T In the string field theory. This line 
of reasoning has already been pursued by Ishlbashi, Kawai and collaborators in the 
context of non-crltlcal string field theory. It would be Interesting to perform the same 
analysis In the CDT model and show that reparametrization under the change of time- 
variable will reappear In a natural way in the model via the operators The results 
should be simpler and more transparent than the corresponding results In non-critical 
string field theory since we have a non-trivial free Hamiltonlan Hq In the CDT model. 
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CHAPTER 



Summary and conclusions 



In this thesis we discussed the incorporation of causality in models of random geome- 
try. 

In the first part of this thesis we introduced causal sets as a simple and instructive 
model of causal random geometry. Using a phenomenological model where causal 
sets are faithfully embedded in a fixed continuum space-time, we proposed a measure 
for the maximal entropy of spherically symmetric spacelil^e space-time regions. Using 
this entropy measure, a bound for the entropy contained in such a region was derived 
from a counting of potential "degrees of freedom" associated with the Cauchy horizon of 
its future domain of dependence. For different spherically symmetric spacelike regions 
in IVIinkowski spacetime of arbitrary dimension, we showed that this proposal leads, 
in the continuum approximation, to Susskind's and Bekenstein's well-known spherical 
entropy bound up to a numerical factor. It was interesting to observe how the inter- 
play between causality and discreteness leads to holography and finite entropy in the 
continuum approximation. 

In the remaining parts of this thesis we focused on the causal dynamical triangula- 
tions (CDT) approach to quantum gravity. 

In Chap. |4] we gave a brief introduction to path integrals. After a short exposition 
of the path integral for the relativistic particle and the relativistic string, we discussed 
problems one encounters when defining a path integral for quantum gravity. 
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In Chap. |5]we described dynamical triangulations (DT) as a non-perturbative defi- 
nition of the path integral for two-dimensional Euclidean quantum gravity. We showed 
an analytical solution and discussed the resulting quantum geometries. An intrigu- 
ing aspect is the fractal structure of the two-dimensional geometries and the related 
Hausdorff dimension a'H = 4. This behavior can be explained in terms of an excessive 
off-splitting of baby universes. We commented on how this proliferation of spatial topol- 
ogy changes leads to problems when trying to obtain a sensible continuum limit of DT 
in dimensions higher than two. 

In Chap. |6] we introduced two-dimensional Lorentzian quantum gravity which can 
be defined through causal dynamical triangulations (CDT). The triangulations used 
in CDT have a global time-slicing without spatial topology changes, i.e. off-splitting 
baby universes. This is in contrast to DT where spatial topology changes are naturally 
present. The global notion of time enables us to define a gravitational Wick rotation. 
We described how to solve the combinatorial problem analytically using transfer matrix 
techniques and performed the continuum limit. As a result one observes that two- 
dimensional CDT and DT are distinct theories. Whereas the continuum dynamics of DT 
is completely dominated by spatial topology changes leading to a fractal dimension of 
c/h = 4, in CDT one has c/h = 2. We briefly commented on results in higher dimensions, 
where the causality constraint is essential to obtain a sensible continuum limit. 

In Chap. [7] we explained several relations between both models. In particular, we 
showed how CDT and DT can be related by respectively introducing or "integrating out" 
baby universes, i.e. spatial topology changes. 

In Chap. [8] we analyzed two-dimensional CDT with non-compact space-times. In 
this situation a semi-classical background emerges from quantum fluctuations yielding 
a simple model of two-dimensional quantum cosmology. 

In the following part of this thesis we introduced a generalized CDT model which 
was formulated as a third quantization of two-dimensional CDT, i.e. a model in which 
spatial universes or strings can be annihilated, and created from the vacuum. In Chap. 
[9] we showed how by introducing a coupling constant gs to the process of off-splitting 
baby-universes one can non-perturbatively solve for the disc function with an arbitrary 
number of baby-universes. In Chap. [TO]this result was embedded in the framework 
of a string field theory (SFT) model of CDT. As an application of the resulting Dyson- 
Schwinger equations (DSE) of the SFT we calculated amplitudes of higher genus. 

One of the biggest achievements of this thesis is the development of a matrix model 
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for two-dimensional causal quantum gravity defined through CDT. After a short intro- 
duction to matrix model techniques employed in DT in Chap.[lT| we formulated a matrix 
model to solve combinatorial problems of the discretized CDT model in Chap. [T2| In 
Chap. [T3] we showed that also the generalized CDT model in the continuum is de- 
scribed by a matrix model. Interestingly, both matrix models can be related by talking a 
continuum limit already on the level of the action. Subsequently, the continuum theory 
also has a large- A/ expansion in which terms with powers A/^''^^ can be identified with 
continuum causal surfaces of genus h. In this sense, the role of N is exactly as in 
the original context of 't Hooft's large-A/ expansion of QCD, namely, to reorganize the 
expansion in the coupling constant gs. 

One advantage of having a matrix model is the computational strength. Until now, 
only transfer matrix techniques have been available for analytical computations in CDT. 
Here we extend the available tools by introducing powerful matrix model techniques for 
CDT 

One possible application is to lift the analytical methods to higher-dimensional CDT. 
In fact, there already exists a matrix model for three-dimensional CDT [137[|139| . How- 
ever, this model is rather complicated which is partially related to the fact that it does 
not allow for off-splitting baby universes. From the discussion above, it seems quite 
plausible that if one allows for the outgrowth of baby universes weighted by a coupling 
gs the matrix model could greatly simplify. 

Another important application comes from the fact that the new matrix model allows 
one in principle to analytically study simple matter models coupled to two-dimensional 
causal quantum gravity, such as minimal models or the Ising model. Several of these 
models are currently under investigation. 

With respect to applying the matrix model to couple the Ising model to CDT there is 
a very interesting feature, namely, that the triangulations appearing in CDT are much 
more regular than the triangulations appearing in DT. In particular, we observed that a 
generic triangulation appearing in DT has a fractal dimension dn = 4, while a generic 
triangulation appearing in CDT has a Hausdorff dimension dn = 2. If one couples the 
Ising model to DT the critical exponents will change from the Onsager-values (see for 
instance |20?1 12051 i206J) to the so-called KPZ-values | fT90l [T9T| . This change can 
be traced to the very fractal structure of generic triangulations in DT. However, we 
have seen that generic triangulations in CDT do not have this fractal structure and we 
expect the critical exponents to be much more similar to those of the flat case. Notably, 
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it has been shown In numerical simulations that the critical exponents for the Ising 
model coupled to CDT are precisely given by the Onsager-values [207l l178[[208J . By 
generalizing the one-matrix model of CDT, as presented in this thesis, to a two-matrix 
model one should be able to prove analytically that the critical exponents are indeed 
given by the Onsager-values. In fact, if possible, this would be by far the simplest 
derivation of the Onsager exponents from an Ising-lattice model. \NorW in this direction 
is in progress. 

We hope to have convinced the reader that the incorporation of causality in mod- 
els of random geometry leads to interesting models of causal quantum gravity and to 
fascinating new results. As described above there are still many interesting challenges 
and we will hopefully report on them soon. 



APPENDIX 




Annexes: Causal sets 



Basic definitions regarding the causal structure 

In this appendix we state some of the basic concepts regarding the causal structure 
of continuum spacetimes. For further references the reader is referred to [2091 1210] . 



where we follow the conventions of [209| . 

Let {M,g) be a time orientable spacetime. We define a differentiate curve X{t) 
to be a future directed timelike curve if at each point p g A the tangent is a future 
directed timelil^e vector. Further, X{t) is called a future directed causal curve if at each 
pe A the tangent is a future directed timelike or null vector. 

Using this definition one can define the chronological future (past) and causal future 
(past) of a spacetime event or a set of spacetime events. The chronological future of 
peM, denoted by /+(p), is defined as 

/+(p) = {q e M\^ future directed timelil<e curve A(t) s.t. A(0) = p and A(l) = q} 

(A.1) 

For any subset Z c we thus define 

/+(Z) = U /+(p) (A.2) 
per 

In analogy to /+(p) and l^{T.) we can also define the chronological past l~{p) and 



/"(Z) by simply replacing "future" by "past" in (|A]T 
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In analogy to the chronological future we can also define the causal future of a 
point peM, denoted by J^{p), by replacing "future directed timelil<e curve" by "future 
directed causal curve" in ( |A.1[ ). The definition of J^{T.), J'{p) and J {T.) then follow 
accordingly. 



Very important for the formulation of the entropy bound stated in Sec. |3]2] is the 
concept of the domain of dependence of a spacetime region. 

Let T be any spacelil<e hypersurface of M. We define the future domain of depen- 
dence of r, denoted by D+{T), by 

D+(r) = {p e M |every past inextendible causal curve through p intersects T} . 

(A.3) 

Analogously, one can define the past domain of dependence of T, denoted by D~{T), 
by simply replacing "past" by "future" in ( |A.3| ). 



Using the definitions above one can then define the future and past Cauchy horizon. 
Let T be any spaceli[<e hypersurface of M. The future Cauctiy fiorizon of T, denoted 
by is defined as 



H+(Z) = D+(Z)\r(D+(r)), (A.4) 



where is the closure of D+{T). The past Cauctiy tiorizon, H~{T), can be defined 
in analogy to ( |A.4| . The different concepts defined in this subsection are illustrated in 
Fig.[3Jl 



A.2 1 Volume of the causal region between two events 

In this appendix we calculate the volume of the intersection of the causal future of an 
event p with the causal past of another event qeJ^{p) in IVIinl^owsl^i spacetime M'', i.e. 
the Alexandrov region J+(p) f] -/ (g). This volume depends only on the proper time r 
between p and q and hence we set Vold(r) = Vol(J+(p) H-/ (g))- 

Note that the volume of a c/-dimensional ball with radius r, Sd{r), is given by 
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Using this we obtain 



T 

Volrf(r) = 2 c/tVol(Sd_i(t)) 
Jo 



t'=:D,t', With Drf = -£^. (A.6) 



Specifically we get Vol2(r) = ^r^, Vol3(r) = ^r^ and VoUlr) = ^r^. 



A.3 I Asymptotics of the generalized hypergeometric function 

In this appendix we obtain the asymptotic expansion of the generalized hypergeometric 
function as given in ( |3.29| ) to first order. 

We use the general expression for the asymptotic expansion of the generalized 
hypergeometric function pFq with p=q which can be derived from [21 1| 

M^i a,; 6, b,. z) = ^) + ^ (^) ) + 



where 



7 



^a/,-6/, and a/, ^ a^ V k=^j. (A.8) 



In the special case of ( |3.29| ) we cannot straightforwardly apply ( |A.7| |, since 3p^i — 3p. 
However, since we are only interested in the first order behavior and ai < aj with ai 7^ a; 
for all 2^jXp, we can obtain the asymptotic expansion of ( |3.29| ) to first order, yielding 



/2 4 d d + l d + 3 2d -I , 



C/ — 1 2 d TT 



-^N--^ + 0{N-^, N^oo. (A.9) 
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Annexes: Causal dynamical triangulations 



Lorentzian angles and simplicial building blocks 

In this appendix, a brief summary of results on Lorentzian angles is presented, where 
we follow the treatment and conventions of [21 2| . 



Since in CDT one considers simplicial manifolds consisting of IVIinl^owsl^ian trian- 
gles, Lorentzian angles or "boosts" naturally appear in the Regge action as rotations 
around vertices. Recall from Section [5J] that the definition of the Gaussian curvature 
at a vertex v is given by ( |5.1| ), 

/<. = {f. (B.I) 
"'1/ 

where = 27r - the deficit angle at a vertex v and is the dual volume of 

the vertex v. Recall that the space-lil<e deficit angle can be positive or negative as 
illustrated in FigurelSTl Furthermore, if the deficit angle is time-like, as shown in Figure 



B.I , it will be complex. The time-lil^e deficit angles are still additive, but contribute to 



the curvature ( |B.1| with the opposite sign. Hence, both space-lil<e defect and time- 



like excess increase the curvature, whereas space-like excess and time-like defect 
decrease it. 

The complex nature of the time-like deficit angles can be seen explicitly by noting 
that the angles Q-, between two edges a, and b, (as vectors in IVIinkowski space) are 
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Figure B.I : Illustration of a space-like (a) and a time-like (b) Lorentzian deficit angle 
at a vertex v. 



calculated using 

cost'/ = 1 ^ ^ 1 , sin^; = ^ ^ ^ ^ , (B.2) 

(a/, a/) 2 (b/, b/) 2 (a,-, a,-) 2(6/, 6/) 2 

where (-, ■) denotes the flat IVIinkowskian scalar product and by definition, the square 
roots of negative arguments are positive imaginary. 

Having given a concrete meaning to Lorentzian angles, we can now use ( |B.2[ ) to cal- 
culate the volume of IVIinkowskian triangles which we will then use to explicitly compute 
the volume terms of the Regge action. 

The triangulations we are considering consist of IVIinkowskian triangles with one 
space-like edge of length squared 1^ = a^ and two time-like edges of length squared 
/f = -aa^ with a > 0. The general argument a > is used to give a mathematically 
precise prescription of the Wick rotation, but it can be set to a = 1 after the Wick 
rotation has been performed. With the use of ( |B.2| we can calculate the volume of 
such a IVIinkowskian triangle, yielding 

a2 



Vol(triangle) = — V4a + 1. (B.3) 



Now one can define the Wick rotation W as the analytic continuation of a -a 
through the lower-half plane. One then sees that for a > | under this prescription 
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/Vol(triangle) ^ -Vol(triangle) (up to a C(l) constant which can be absorbed in the 
corresponding coupling constant in the action). This ensures that 

yy ■ e'^R^ss^^'^'") e"^R<=sse(T"")^ a > ^. (B.4) 

In the following we set a = 1 again. Generalizations of this treatment to dimension 
d = 3A can be found in | [T36| . 



Analyzing Hamiltonians for causal quantum gravity 



In this appendix we analyze the properties and spectrum of Hamiltonians which one 



usually encounters in causal quantum gravity such as discussed in Sec.|6|5| Generally 
the Hamiltonian appearing in this context are of the form 



d 



(B.5) 



where ci and C2 are real constants. In the following we investigate this generalized 
class of Hamiltonians. 

It is immediately checked that the Hamiltonian ( |B.5[ ) is selfadjoint on the Hilbert 
space H = C^{R+, dfi{L)), where the measure dfi{L) is given by 



dfx{L) = Lf'dL, fx 



Cl 



(B.6) 



Further, for the boundary components of the partial integration \n {Htp \ cj)) = {tp \ Hcf)) 
to vanish, the wave functions on the Hilbert space must satisfy 







(B.7) 



To have a Hilbert space with flat measure, one can pull the measure into the wave 
function by substituting Mf{L) = L'2iIj{L), where ip{L) is the wave function corresponding 
to ( |B.5[ ). Commuting /. 2 through the Hamiltonian gives 

d (ci - C2Y 1 

I' 



d 
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4ci 



dn{L) = dL 



(B.8) 



This pulling in and out of the measure is similar to what one does when introducing 
marl<ings on the boundary loops. 
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We now consider the eigenvalue problem of the Hamiltonian (|B.8[), 



Let us perform a change of variables and wave functions, 

L=^^'. cD(^) = y^vl; (^^^ , (B.10) 

where the latter guarantees a flat measure dfi{ip) = dip for ^'(v?)- The eigenvalue prob- 
lem then reads 

where we set uj = \/oJ\ and A = 1-4/1^ g (-00, 1). This is nothing but the eigenvalue 
problem corresponding to the Hamiltonian of the one-dimensional Calogero model, 

Note that the parameter range A = l-Af/ e (-00, 1) is the maximum range for which 
the Calogero Hamiltonian is selfadjoint. Two-dimensional Lorentzian quantum grav- 
ity with open boundary conditions corresponds to the value A = 1, whereas for cir- 
cular boundary conditions it corresponds to A = -3. The first connection between 
two-dimensional CDT and the Calogero model was established in [127J, where a gen- 
eralized CDT model was introduced which continuously covered the parameter range 
O^A^l. 

The spectrum of the Calogero Hamiltonian is well-known, further it can be related to 
the spectrum of the radial solution of the three-dimensional Schrodinger equation with 
potential U{r) = Di/r^ + D2r'^ [21 3). Explicitly, the solution of the eigenvalue problem 
( |B.1 1| | can be obtained by doing a variable transformation L = ujip^, yielding 

I^ + ii-il+^ + ^V(I) = 0. (B.13) 

5/.2 25/. 4 2uj 16/.; ^ ' ^ ' 

Note that L is just a scaled version of L, so we could have directly arrived here from 
( |B.9| ), using the wave function transformation from ^{L) to <^(L). From (jBTSj one can 



see that asymptotically for /. ^ 00 the solution should be proportional to exp(-/./2) 
and for small I ^ it behaves like IV4+a*/2 ^q^q ^j^^^ these are the only asymptotics 
compatible with the requirement ( |B.7| on the wave functions. Hence, we make the 
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following ansatz for the wave functions <^{L) oc exp(-I/2)I^/'^+''/2?7(I). Substituting this 
into ([BJ3) yields 

WiD + (1 + /X - + (^A - + 1)^ r^{L) = (B.14) 

This equation is known as Kummer's equation |214], whose solutions are the confluent 
hypergeometric functions, 

r]CL)= iFi(-n;l + (B.15) 

where n= - |(1 + /x) must be a nonnegative integer. In this case the power series of 
the confluent hypergeometric function is truncated to a polynomial of degree n, namely, 
the generalized Laguerre polynomials [57], 

with 

- r(n + /i + i) z' 



yielding the final expression, 

.r,(-. . . i: z) ^ (;;) (^,,)(^;,,,^,,) ^^ (B.1 9) 

Since n = + must be a nonnegative integer, the spectrum of the Calogero 

Hamiltonian H^^'^sero ^g^^g 

E„ = w(2n + /i + 1), n = 0,1,2,... (B.20) 
The corresponding eigenfunctions are given by 

= C„e-¥'^'^l+M,Fi(-n;l + ^;cuv^2)^ d^i^^) = d^. (B.21) 

Further, since the Laguerre polynomials are orthogonal functions, the wave functions 
form an orthonormal basis of the Hilbert space, where the normalization constant is 
given by 

^"-Vr(n + i)r(;. + i)- (^-"'^ 
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where we used the following orthogonal relation for the generalized Laguerre polyno- 
mials [57] . 

dLe-^L'^l'tminL) = { , , ^ (B.23) 

which is valid for /i^O . 

Let us now come back to the original problem, analyzing the Hamiltonian 

) = -Ci/.— -C2— + A/., dii[L) = L^dL, /"=|-1- (B.24) 

The spectrum is obviously the same as for the Calogero Hamiltonian, hence from 
we get 

= y^(2n + /i + 1), n = 0,1,2,... (B.25) 



Further, the eigenfunctions of H[L, di) can be obtained from ( |B.21| by a simultaneous 
variable and wave function transformation, yielding 

i,n{L) = Ae"^SFi(-n; 1 + ^^■ 2^^/-). df^{L) = L^dL (B.26) 
where the normalization factor is given by 



-tf±i /4A\ " / r(n + /i + l) 



Having obtained the eigenfunctions for the class of Hamiltonians H{L,dL), we are 
now able to give an explicit expression for the finite time propagator or loop-loop corre- 
lator, 

GA(/-i,/-2;7-) = (/.2|e-^^|/.i). (B.28) 

Since the eigenfunctions form an orthonormal basis of the Hilbert space, we can insert 
the unit operator / = Xln^o !")("!' yielding 

oo 

Ga(/.i,/.2;7-) = J2{L2\n)e-^^''{n\U) 

oo 

= 5^e-^^"^:(/.2)^„(/-i) (B.29) 



n=0 
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Upon inserting the energy eigenvalues | |B.25[ ) and eigenfunctions ( |B.26| ) into ( |B.29| ) one 
obtains 



Ga(/-i,/-2;7-) 



n=0 



- 



,„+i£+l r(n + /i + l) 



Ci / 



r(n + i)r(/i + i)2 



XiFi(-n; 1 + /i; 2^^/-!) iFi(-n; 1 + /x; 2y^/-2), 



(B.30) 



where we used the notation z = e^^^^^. To evaluate the above summation we make 
use of the following quadratic relation satisfied by the confluent hypergeometric function 
(571, 



f)=0 



r(n + i)r(/i + i)2 



1 



z(x+r) 



^— e 1- (xyz) ' 



/xyz 



(B.31) 



where /^(x) denotes the modified Bessel function of the /x-th kind. Gathering all terms 
together leads to the final expression for the finite time propagator 



A ^ g-y^(Z-i + Z-2)coth(v^T) 



sinh(v^CiA7") 



ci sinh(v^Q7\7) 



(B.32) 



For /i = ci = 1, | |B.32[ ) corresponds to the "unmarked" propagator for the two-dimensional 
CDT model with circular boundary conditions, as obtained in | |6.33| . 
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